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ABSTRACT 


A procedure is presented that allows an engineer to coordinate 
protective devices consisting of fuses, circuit breakers, and over- 
eurrent relays based upon load and fault currents. The resulting 
digital-computer program makes many of the decisions requiring very 
little prior in-depth knowledge of protective devices and power systems 
fPemercs utilization. 

The coordination procedure consists of six separate programs 
that must run in sequence, as the output from one program forms the 
input for the next program. This reduces in size of internal computer 
storage required. 

The theory is presented with an explanation of each algorithm 
followed by a basic flow chart. An example computer run is presented 


merche end. 
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eee LN PRODUCTION 


A prime objective of electrical engineers working on the design, 
operation, or maintenance of electrical power systems is to provide 
reliable, efficient power to their customers. The distribution system 
has always been a weak link in the power chain. Often, the catas-— 
trophic results of a fault are easily seen, but the cause may be 
extremely difficult to find. A distribution system, either in the 
utility or the user's facilities that is not adequately protected from 
electrical faults is both dangerous and costly. An unsatisfactorily 
protected system can cause considerable inconvenience and damage to 
equipment. Adequate protection is not a set definition. Each load 
system is different with its own set of boundary values. 

In a distribution system, the flexibility of design reduces 
near the generator or the load. Usually the generator or load protec- 
tion is designed as an integral part of the equipment by the manufacturer 
and any changes should only be undertaken with their counsel. Moving 
away from the generator or load,the priority shifts from load pro- 
tection to ee eer protection. It must never be forgotten that 
the load, distribution, and generation systems must be protected in 
Such a way) that a device that operates satisfactorily in one system 
cannot have a detrimental effect on another system. In this thesis 


the distribution section is attached from an infinite-bus power 


source to the load-protection device. The programs provide the 





jeepers setting of a chosen device in order to provide good protec- 

tive Me nation. Proper coordination is achieved when a device 
@losest to the fault clears the fault or overload condition before 

the rest of the system is disturbed. In order to achieve this end by 
manual means, extensive mathematical calculations and experience are 
required. As a direct result,proper coordination is often not achieved. 
The use of digital computers to solve the coordination problem makes 

a more exacting solution. However, computers with large internal 
memory are required to solve these problems. Several approaches 


to computer solutions have been published [1-4]. 


Problem Statement 


The purpose of this thesis is to develop the procedure and an 
interactive computer program for a generalized case of device coordina- 
tion. This coordination should include fuses, relays, and molded case 


Circuit breakers. The phases of research are as follows. 


1. Develop a computer procedure to compute the one-line program 


per unit values. 
2. Develop a computer procedure to perform a load-flow analysis. 


Se DVevelop a computer procedure to compute line-to-line ground 


and three-phase faults. 


4. Develop and apply procedures for mathematical representation 
pmeliness/R data, fuse, molded-case circuit breaker, and 


relay-curve time-current characteristics. 





Pee Wevelop a computer Precedure to coordinate the devices 


from phase 4 with the data derived from phases 2 and 3. 


Although the coordination programs are written for the PDP1O 
computer, it requires no more than 25K bytes of internal memory and 
a disk unit to solve any system of 30 buses and 100 elements. Almost 
any large system can be reduced to its smaller parts composed of 30 
buses or less. As each segment is coordinated, it can be represented 
by a load and its protective device. In this manner a 300-bus system 
could be divided into smaller segments and in 11 runs would be com- 
pletely coordinated. The step-by-step approach used here allows for 
insertion or extraction of data before proceeding to the next Step. 
fmeomee 1-] is a basic flow chart of the coordination program. Here- 
after, it will be referred to by its individual parts. 

When the coordination process has been completed, a load flow 
Study and a fault analysis has been performed as well. The major 
assumptions that have been made are listed here. Other minor ones 
are included in the corresponding chapters. These assumptions have 


been made in an effort to keep the solution as general as possible. 


1. The incoming line is assumed to be an infinite bus. 


2. Assymetrical currents are neglected. Devices are coordi- 


nated on the basis of symmetrical fault currents. 


3. Only synchronous-motor contributions to the fault 


currents are recognized. If a large induction motor 


Ser ee => 


Se ee oe - 
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Pieure 1.1. Flow chart of the coordination procedure. 





Figure 1.1. 
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Yes 


is present, it may be added as a fictitious synchronous 


motor. 
4. Line-charging shunt admittance is ignored. 


5. Prefault current and voltages are neglected in fault- 


current computations. 


6. Different sizes of fuses have the same time-current curve 


shape. 
PeeeeOoverioggd conditions are 125% of full load. 


8. The system has a balanced three-phase load. 


In this thesis each phase of research in the problem statement 
has a chapter devoted to it. After describing the computational 
methods used in curve smoothing, the coordination problem begins with 
data input. Next the fault currents are solved followed by the 
actual coordination program and device response graphing. Finally, 


a sample distribution system coordination problem is described. 





2. CURVE SMOOTHING AND DIGITIZING 


General 


In this chapter, the methods for using a data digitizing 'numonics" 
device and the curve smoothing program of Appendix A are explained. 
The various device equations that were found that adequately represented 
the particular device's time current curve are also described. 

A method was described in Wagner [1] for performing this by 
using the PDP10 computer. This least-square curve-fitting program [1] 
was expanded to eliminate any external subroutines that may not be 
available. The coordination program makes extensive use of polynomial 
functions whose coefficients were derived by using the least-squares 


curve-smoothing technique. 


Curve Smoothing 


Electrical Cables. Coefficients for the mathematical modeling 
of resistance and reactance were taken from data for mine power 
cable [5]. This is used as only one type of cable or wire. Any other 
type may be modeled by the same procedure. For instance, either actual 
cable or line data could be used directly [6],or the theoretical resis- 


tance and reactance for each size could be calculated [5-8]. 





Fuses. Two types of fuses are used in the coordination pro- 
meee rhnese are current-limiting and solid-material boric-acid power 
fuses. In each case,the curve shape is assumed to be the same for 
different sizes. The low-current asymptotic value is called the 


"ofset.'' Each fuse has a minimum and a typical operating time. In 


fuse coordination, both values must be used in describing the operating 
Pmevelope Of a particular fuse. The higher typical operating time 
Maeve 1S referred to as "Ofset 1." By this means, any fuse can be 


completely described by four equations, where: 


n 
- - alae 
Log a Ce eee Current -Ofset) 2 G11) 
j=o 
i k 
Poculee— ("2 4)(Leo Current —Ctsetl) — 2 OE) 
yp 'x=0 * 


Fuses are selected also by using mathematical models. To 
select a fuse, the program uses the overload-current value as the 
variable in the polynomial function. The result is the device number. 
In curve smoothing, variables can assume no value less or greater than 
those originally used as data for the curve-smoothing program, or large 
Perors can result. 

Molded-Case Circuit Breakers. Like the fuses discussed above, 
each molded-case circuit breaker consists of an operating envelope 


bounded by a minimum value and a maximum value. Unlike fuses, however, 
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a table look-up approach is used in selection of a device and its 


setting because of the smaller number of sizes and magnetic-trip settings. 

Relays. Wagner [1] discusses the representation of the family 

Smerelay Curves for each type with one polynomial equation comprised 
a th ee 

Sietour coefficients. The 0 power coefficient was left out as a 

variable for coordination. Unfortunately, this method proves unsatis- 

factory when working with more than one relay type. It was found that 

: nee th 

each relay could be represented by five coefficients. The 0 power 

coefficient is necessary in that it provides the time displacement 

Mecessary for minimum time dial setting. The curves were derived 

using the procedure below. 

Westinghouse Company type relays were used. An average curve 
was taken as time-dial 4. This curve was reduced to a polynomial 
function. The time difference between time-dial 4 and time-dial 1/2 

: ; : th 

taken at 20 times the pickup value is subtracted from the 0 coef- 
ficient. Twenty times pickup was a good value for computing the time 
difference as the time value changes very little for increases in 
current beyond that point. The resultant curve was the average shape 
of a time-dial setting of 4 now shifted to minimum setting. This allows 

; : th 
for time values computed during coordination to be summed with the 0 


coefficient. The curves are then represented as follows: 


nh . 
Log Time = ( a, (log(Current Value)?) + 
j=0 


Log (Time Setting) 





ey 


This allows the use of more than one type of relay in each 


coordination as it separates time setting value from relay type. 


Curve Digitizing. The device's time versus current curves 
were digitized by using a NUMONICS electronic graphics calculator 
(EGC). The NUMONICS EGC is designed to translate graphic information 
on any medium into digital information. Complete information on its 
usage may be found in -NUMONICS Electronic Graphics Calculator [9]. 
The NUMONICS EGC then provided a convenient way to enter X-Y coordinate 
data representing any curve onto disk. Since this device works in 
length in inches, the data has to be normalized to its proper units. 

A simple program DIGIT.F4 was written to normalize and store 
onto disk the NUMONICS EGC information. The program flow chart is 


shown in Figure 2.1. 


Curve Smoothing. The curve-smoothing program contained in 
Wagner [1] was used as a basis for the curve smoothing program DEV.F4. 
Unnecessary parts were removed and a double precision matrix inver- 
sion routine was added. This prevented the necessity of using 
subroutines contained in the computer's SSP routine files. The 
inversion routine was converted from The Pennsylvania State University 
Computer Center's matrix inversion routine CMINV which inverts a 
complex square matrix. This routine in turn was converted from the 
SSP routine MINV [10]. 

Since it has its own inversion routine, this curve smoothing 
program can be easily expanded to any number of points and any 
dimension of fit by changing the dimension statements. The flow 
chart and procedure remains the same as in [1]. A complete program 


fisting may be found in Appendix A. 
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Srere 2.1. Curve-digitizing fiow chart. 





mopure 2.1. 


Continued. 





13 





| - 


» Oper CAs? 





Se 


14 


summary 


Polynomial functions describe the various device-response 
m@gmves. line X/R data, and fuse-selection criteria. The coefficients 
for these polynomial functions are found from data entered using a 
digital-data device or manual input. This data is reduced by a curve- 
smoothing program to the desired function. 

Micestart or Ene COOndination problem begins in the next chapter 


with line and load data entry. 
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eee NE OF DISTRIBULTION-SYSTEM DATA 


General 


Doubly sthe first step in performing any calculations on a 
distribution system is to reduce all system values to per unit. The 
methods used to do this in the program will be described in this 
chapter. Errors can be introduced easily when converting to the per 
unit system using manual calculations. Program INPUT.F4 was written 
to perform this function as well as computing line resistance and 
reactance values and eeceee all system data on disk for use in sub- 


sequent programs. 


Per-Unit Reduction 


The program uses the following formulas to convert to the per- 


unit system [8]. (See Table 3.1 for variable definitions used in this 
fnesis. ) 

Z = (KV )* /MVA 

Base Base Base 

ee nase 

I = KVA /KV 


Base Base Base 





maple 3.1. 


Zo) 


ry 
Pq 


Pq,rs 
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Variapler detiniwtions. 


element impedance 


impedance in per unit 


Suberancient reactance 


real power bus p 


reactive power bus p 


complex voltage bus p 


complex current bus p 


Y bus element bus p to q 


admittance value from primitive admittance matrix 


Z bus element bus p to q 


impedance value from primitive impedance matrix 


fault impedance 


fault admittance 








Ly 


KV 2 KVA 
7 = 7 ( Base EVEN) ( Base new ) 
P.U.new P.U.given KV KVA 

base Mew Base given 


Utilizing the above formulas it is possible then to represent all 
elements in the distribution system on a per-unit basis. Loads need 
to be included also in the one-line representation. 

The program allows loads to be represented either by induction 
motor horsepower, synchronous motor horsepower, kVA, or current. 

Loads are entered in the following format: Bus number/horse- 
power, kVA, current/power factor. The program assumes 100% efficiency 


so one must use: 


Cera = HP/Efficiency Ge) 


Tf information on efficiency and power factors are unknown, the follow- 


ing equations may be used instead of 3.1 [9]. 


For induction motors, use: 


HE eas eae x lL O27 
Power Factor = .86 b322)) 
For synchronous motors use: 
Unity Power Factor HP = 1.139 x Horsepower (32) 


.8 Power Factor HP = 1.18 x Horsepower (336) 





ie/ 


KV : 2 KVA 
Z =, ( Base pWvien, ( Base new ) 
P.U.new Pe use iven KV KVA ; 
Base new Base given 


Utilizing the above formulas it is possible then to represent all 
elements in the distribution system on a per-unit basis. Loads need 
to be included also in the one-line representation. 

The program allows loads to be represented either by induction 
motor horsepower, synchronous motor horsepower, KVA, or current. 

Loads are entered in the following format: Bus number/horse- 
power, kVA, current/power factor. The program assumes 100% efficiency 


so one must use: 


Jaan = HP/Efficiency (orale) 


If information on efficiency and power factors are unknown, the follow- 


ing equations may be used instead of 3.1 [9]. 


For induction motors, use: 


HP entered =a XM ode Oz 7 
Power Factor = .86 (322) 
For synchronous motors use: 
Unity Power Factor HP = 1.139 x Horsepower (35) 


.8 Power Factor HP = 1.18 x Horsepower G36) 
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Peemalac for Line Resistance and Reactance 


Element data is entered in the following format: line bus/ 
load bus/cable length in feet/wire size/transformer P.U. impedance/ 
size in MVA/line voltage in KV/load voltage in KV. 

Wire sizes that are entered in American wire gauge (AWG) are 


converted to MCM wire gauge by using the following formula [1]: 


2 


(36 AMG), 


MEM —C. OOa) Cl. 229) x 1000 C7) 


Once the wire size is in MCM, the line resistance and reactance are 
computed using the polynomial coefficients that were found by the 


meecedure in Chapter 2. 


n 
Resistance Or reactance = 2. 


a. [log(McM) ]/ (3.8) 
j=0 7 


miesprogram should not be limited to just one wire type. The program 
does this by reading wire coefficients Z above from disk so that each 
time the program is run it can compute element data using different wire 
types until the entire system has been completely entered. 

Transformers become very important in that they can become the 
highest source of resistance and reactance. By their connections and 
Peters impedance, they determine the magnitude of zero-sequence current 
that will flow in the particular element. The program allows for this 
Bamiiemande writes this imtormation in the particular disk files for use 


ater on in fault analysis. 
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Transformer inrush and withstand values must also be calcu- 
lated. These values must be used in coordinating protective devices. 
These devices must pass inrush currents but open under conditions of 
exceeded transformer withstand currents. The algorithm in the program 
is based on average values. If test data is available this should be 
entered directly by changing the transformer data disk file directly 
after the program is run and before proceeding with coordinations. 


Inrush and withstand values are computed as follows [11]: 


Inrush = 12x(Transformer Base Current) (3.9) 
Withstand @2 Sec= 25x (Transformer Base Current) (3510) 
Withstand @5 Sec = 14.5x (Transformer Base Current) Cale) 


Transformer withstand is some value between the two second value to the 
five second value and depends upon the transformer impedance. It is 
fairly easy, however, to have the device protect the complete range 
of withstand values and the program is written to accomplish this task. 


Synchronous subtransient reactance is computed as follows [11]: 


(PersunitenvA) x (.17) = 4X 


iX, (Ge 12 
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One of the problems encountered in computing transmission line 
impedance is a reasonable calculation of the zero-sequence impedance. 
The computer program assumes that this value js two times 
the positive-sequence impedance [8] if nothing is entered when the 
computer asks for the zero-sequence impedance multiple. If a different 
value is entered, the computer multiplies the positive sequence reactive 
mame by this multiple to arrive at the zero-sequence value. 

Peipdcic flow chart follows in Figure 3.1 with a complete program 


isting in Appendix B. 


Summary 


The first step toward solving the coordination problem is com- 
pleted. All element data has been converted to the per-unit 
representation. The next logical step is to perform a load-flow 


analysis to determine the full load currents in each element. 
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Figure 3.1. Per-unit calculations flow chart. 
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4. LOAD-FLOW ANALYSIS 


General 


In order to establish overload conditions on each element, 
the load currents must be determined. This is done by using a load 
flow analysis. The Newton-Raphson Method using Y-bus [12] was chosen 
because it required the least amount of iterations to obtain a solu- 


tion to the load-flow problem. 


Load Flow 


In this method Y-bus must first be formed from the primitive- 
impedance network. When mutual coupling is ignored, Y~-bus can be 
formed easily by algorithm. Any system can be completely represented 


by a set of equations [8,12]. 


bus bus bus 
i = (Y )(E ) (4.1) 
bus bus ; : 
ES and I are column matrices with one entry for eacn bus 
: bus . : 
im the system. Y is a square matrix of order equal to the «umber 


of buses in the system. Therefore, once all bus voltages are known and 


Petts 1S formed, matrix multiplication gives the desired load currents. 
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However, there is no direct approach to obtain the bus voltages. 
The only approach to use is that of iteration. The Newton-Raphson 
method accelerates the iteration process so that normally an accurate 
answer is obtained after three iterations. 


The power at bus p is (from 4.1): 
eon 
p> *p Pp Pp 


~ 


wn 
Piers YG Ce?) 
eee Sy) Parc 
Be eet Y = 6 =4R (4.3) 
P =p” Pp pq pq” pq 


So 


n 
P = £ fe (e G +f B )+(£ G -e B_ )} (4.4) 
Bre Pq pq q pq q 


al PG qd Pq 


n 
= LoS 
Ss = one +f£ B CEC =e. RB )} ( ) 


)-e 
1 Pd 4 Pq pq pq -q pq 


The Newton-Raphson method uses the Jacobian matrix of partial- 


differential equations as follows: 


AP 
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Ae 


(4.6) 


The elements of the Jacobian matrix are found from the following 


equations: 


oP 
EPL > 
Je 


eG +f B 
P PP Pp PP Pp Ppp 


=ieb 


n 
ee cen Geter B 


q=l 
q=p 


Sige oe 








a 
Fh 
do 

il 
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Ph 
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OQ 
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eee et GC F 
P PQ Pp pg ee 


Th 
=e B +Z.G -e B +2 (£ G -e B SY) (A958) 
P pp fp pp Pp pp aa eee 
q#p 
=e B +f G_ aq#p 
P PG Pp PY 
Th 
=f CG -£ G 4+Z B -23 (£ G -eB ) (4.9) 
P PP P PP E&p PP q=l Gq Pd q PY 
q#p 
=e G +f B qip 
P PG P PG 
Th 


eG +Z._B -eG +2 (eG +f B_) (4.10) 
Pp pp fp pp Pp pp gen “| PS ees 


p#q 


I 
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Bus #1 or the swing-bus voltage is known. The swing bus is 
the bus that will provide whatever power the system requires. The 
program begins by forming Y-bus by using the primitive-impedance matrix 
generated by INPUT.F4 of Chapter 3. Since mutual coupling is ignored, 
the y-admittance matrix can be found by inverting each member of the 
Meee aryix found in disk file TEsT.paT. Y°"® matrix is formed 
using the algorithm below. 

Diagonal elements are formed as the sum of all admittances 
connected to that bus or node. Off-diagonal elements are the negative 
of the admittance connected between the buses or nodes. 

If more than 10 iterations are performed, the program stops and 
informs the user that either an error has been made in design or data 
entry. Upon successful completion of the program, a load-flow analysis 
may be obtained by printing file DISK.DAT. The program puts the load 
currents in the proper format and file for later use in the coordination 
program. A flow chart follows in Figure 4.1, and a complete listing 


is given in Appendix B. 


summary 


This load-flow analysis obtained can be very useful in analyzing 
the distribution-system design. It can be used to point out locations 
for power factor correction as well as potential problems due to under- 
seged wiring. 

The load currents found here are stored on disk to be used in 


the coordination program. The next step before beginning the 
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PeemdinalLiOn program is to solve for all fault currents. Zbus must 
be formed to compute these currents. The methods used to do this 


Pollows in the next chapter. 
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Figure 4.1. Newton-Raphson load-flow chart. 
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>. FORMATION OF IMPEDANCE-BUS NETWORK 


General 


The formation of the impedance bus network Cates) by 


feeorithm is the first step toward solving for fault currents. In 
this chapter this algorithm will be discussed as well as other problems 


encountered in matrix compression. 


Formation of Z-Bus 


_ Once the one line diagram has been completed, the primitive 
impedance matrix that was found in INPUT.F4 is a partially filled 
matrix. However, it is not in the proper format for the formation of 
Z-bus as it contains no coupling. The primitive impedance matrix must 
change from one of order equal to the number of buses to one of order 
equal to the number of elements. This matrix becomes a very large 
Sparse matrix. The memory required to store this matrix becomes too 
Harege to fit in a small computer. A method was derived to only store 
myeenon Zero values and their location in the matrix. Since the matrix 
fmestrianeular only the top half is stored (Figure 5.1). 

Diagonal elements found in INPUT.F4 are stored in "ZD." ZD 
is a three-dimensional matrix, ZD(200,2) where (X,1) is the element 


number and (X,2) is the diagonal value. 





| * | ZD(200, 2) 
| ! | ZU (2001) 
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Elements X,Y,Z 


(A) lr = 
ea / Wy 
/ 





gui 200, 1) 


[IM | r200, 29 





Number locations matrix (A) = (# elenentey— 2 
Number of locations of alternate method = 8 X (# elements) 


Figure 5.1. Transformation of Square sparse primitive impedance. 
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The off diagonal mutual coupling Zm is a three dimensional 
matrix also where (X,1) is the location value At+tjB. Element A is 
coupled to element B. Location (X,2) is the coupling value. Fortun- 
ately, partitioning makes it necessary to assemble a Square matrix for 
inversion of only the coupled elements. Since coupling only exists in 
zero sequence, this applies only to the zero-sequence Z-bus matrix. 

finethe formation of the Z—bus matrix, Bue 1 is the reference 
bus for all calculations. Z-bus is formed by starting with element 
one and proceeding with each element in turn. An element is either a 
fink Or branch. A link means that both buses have been used earlier; 
a branch means they have not been used earlier. 


Like Y-bus in Chapter 4 any system can be represented by 


Ebus? <4 = Zbuse* Thus) * 


Formation of abe can be formed either by the incidence, 


network-matrices method [12] or by algorithm. Formation by algorithm 
is much simpler and is used here. 
After forming the primitive impedance network, it is a straight- 


fOrward application of the equations of Table 5.1. 
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fibre 5.1. Equations for the formation of ZBUS. 


Aa 
Add p-q Mutual Coupling No Mutual Coupling 
Branch eee (ZZ. ) eae age 
qi Piae pag tS ori. si qi py 
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eee 7 
7 = 7 + = Paes Pg tS 7 =7 42 
qq Pq y__,ts qq PC Pq>Pq 
Pq 
vie 2.) 
Link Dae = ea ee pg rs Ti 51 oe = Za 
*pqsPq 
= he Ze aati i #1 
es ec clacal 5, = Se 
Jak pl gi y 11 pay a 
Pq,Pq a2 
P9>»PY 
a. : th 
Elimination of 1 Node 
Fein Dit: 
re sot i : Pl 
Ae uedintcaj=— Zee (bperore Elimination) — —— 
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For example, a primitive impedance matrix looks like this 


noeertive elements: 


1(1-2) 2(2=-3) 3(3-4) 4(1-4) 5 (2-4) 


H elt. 5) 75 
2 ae 2145.5 
Z2° = 3 wie dO 34 vee [2°}-* 
4 eo 144.5 48 
5 44 48 es 


where the off diagonal values represent mutual coupling between ele- 
Memrs. Using this primitive matrix and the equations in Table 5.1 
yields Zbus. 

ZBUS.F4 is one of the longest programs in the series. The 
program asks for the zero-sequence impedance for all loads in the 
system and the mutual-coupling impedances. If the zero-sequence 
impedance of a load is not known, the program assigns it the same value 
as the positive sequence impedance. 

All loads that are added are links, and only loads that are 
added are those of synchronous motors. The subtransient reactance 
is used instead of its actual impedance since the Z-bus is going to be 
used for fault calculations. Any link from a bus to the reference bus 
will cause a fault current. Prefault voltages and currents will be 
ignored as will all contributions from loads other than synchronous 
mentors. A flow chart follows (Figure 5.2) and a complete program list- 


mee in Appendix B. 
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Beeetre 5.2. Continued. 
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Write 
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Mepure 5.2. Continued. 
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Summary 


le : : 
eee is formed by algorithm allowing for mutual coupling 
between elements. The next step in the coordination procedure is to 


peers ror fault currents. This procedure follows in the next chapter. 
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Cree WiE CALCULATIONS 


eeneral 


The final program required before executing the coordination 
peeeram 1s FALT.F4. This program computes all the possible line-to- 
ground and three~phase—fault currents and prepares all other load 
@eeaeftor entry into the coordination program. The procedure to 


accomplish this will be discussed in this chapter. 


Fault Currents. The equations describing the fault currents 


meed to be derived [12]. 


Three—Phase Fault. For a fault at bus P: 


Zz 


ol = 012, 012,-1 012 


O12 
(U+Z Y 


*F pp r ? PCO) on 


where U is unity matrix, a is prefault voltage, and Yp is the 


fault admittance matrix for a three-phase fault. 


0 
ae ae ome 0 (6.2) 
0 








Substituting (6.2) into (6.1) 
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ail 
feo NG 142, 0 
012 
I = 0 1+2Z ¥ 3 
P(F) Vy pp £ 2 
i 0. 1+Z 0 
ppé 
oy 0 70 ~ 0 
142 
2 -y Io 1 9 L / 3 
P(F) F if 
LZ y 
pp -F il 
eal 5 0 
1+Z 
pp °F 
0 0 0 0 
012 I eS 
1 (F) oe [ee 7 0 ¥ 3 (G3) 
EZ, y 
Ppa 
A F iL 0 
2 
1+Z 
pp °F 
1 2 
Since Zpp = Zpp. and Zz. = 1/y, 
2 y 3 meron is = 6 (6.4) 
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Using the same method the following equations are found [12]. 


/ 
rl = : “FF Po 
Herr 
and 
(1) 
_ Lee F 
Base 3 (Es (9) x Bape cere 
Z Ate - 
f “PP 
O..2 
i#P E, 0 E Dee 0 (6.6) 
A p(f) 
These equations are identical to those found for a three-phase 
grounded fault [12]. 
. 3 ot ae : O12 
The line-to-ground fault admittance matrix is given by Ya 
below. Substituting aoe into (6.1) the equations below result [12]. 
1 al i 
you2 Fy yy 
F 5 
if al i 
VS 
mole po) (6.7) 
: t O30 


Ze 2 
P PP 
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ee 0) 
a pp 
Y3E ° 
eas P (Ding Biay. (6.8) 
Pp z (0) 57 (1) 3, pp pp F 
PP pp F a (1) 
i) 
0 Z (0) 
vo ae 
O12 - are ma p (0) (1) 6 
mi 60 10) | Pl > @,, ©,,. | we 1°” 
0 pp pp ae . (5 
ifp : ip 


The value of 24, is entered from the keyboard when asked by the program. 
All of the equations are available now to solve for all fault currents. 
All post-fault bus voltages are solved. Using these bus voltages 

each element fault current is determined. These currents are stored 


onto disk so that they may be used in both the coordination and graph- 


ing programs. 


Summary 


The programs thus far have solved for all load and fault cur- 
mits DUL NO Coordination has been discussed. Now the coordination 
problem is ready to be addressed. In the chapter that follows, the 


protective devices for the distribution system will be coordinated. 
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meeure 6.1]. Flow chart for FALT.F4. 
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(ae LOU CETVE-DEVECE COORDENATION 


General 


In this chapter,the protective-device-coordination problem 
will be discussed. A computational method is put forward that allows 


a digital computer to solve the problem. 


Device Coordination 


One of the best methods to approach coordination is a step by 
step device selection and setting working from the load toward the 
source. A load is the best place to start because there is not much 
meexibility for device setting at a load. The load flow analysis of 
Chapter 4 provides the overload values for the various loads. This 
enables a definite overload setting for each device protecting a load. 
Since overload currents are a primary concern at these buses, all fault 
settings can be placed at six times the full load currents. This should 
allow for normal starting currents. Care must be taken so that an 
instantaneous tripping element will not open under normal starting 
conditions. From the loads it is a straightforward process to place 
a line side device curve next to its load side curve to insure that 
proper coordination is achieved. By taking this approach, each device 


depends only on one other device. 
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This program 1s written with the option of entering two devices 
fomeach e€lement—-one located at each bus. From this approach, 
Obviously two possibilities can occur. Either the device is on the 
line side of the element or on the load side of the element (Figure 7.1). 
If it is the line-side device, its setting depends only on the one 
load~side device. It must be coordinated for the fault current at 
the load bus and not the line bus [1]. Devices must be coordinated 
for a minimum possible fault current which will be discussed later. 

If the device is on the load side and it does not protect a 
Meade then it must protect a following bus. This device must then be 
meerdinated with all line-side devices connected to the bus it pro-— 
tects. This load-side device must be coordinated for fault currents 
at the protected bus. 

Devices should be set to (1) open under overload conditions, 
(2) open under minimum fault conditions with proper coordination, and (3) 
open under maximum fault conditions with proper coordination. Often, 
however, all three conditions cannot be met simultaneously unless the 
proper devices are chosen. Criterion (1) is not as important as 
meteerion (2) or (3) when in the distribution system not at a load bus 
since an overload condition below fault conditions should not do damage 
to the wiring. For this reason wire sizes must be checked to insure 
that upon final coordination they can withstand the overload and fault 
Semditions that can occur [5,13]. At this point it may be necessary 
to begin the entire coordination procedure again if an element cannot 
withstand the short-circuit conditions. Allowable short-circuit cur- 
rents can be checked by using the final coordination graph and the fault 


@iterentS Obtained in FALT.F4. 
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Bus 1 Bus 2 


(a) Coordinating line side device of an element. 





(b) Coordinating load side device of an element feeding a bus with 
load elements connected W, Z. Y. 


meoure /.1. Two possibilities for coordination element Z. In 
Figure (7.la) coordination device A according to curve 
jdeviceosbanpmuamnault bus 2, In Figure (/.lb) coordina- 
tion device B depending on load flow current element Z. 
Ugorndimaasedewern device curves C, D, E, up to fault at 
pus’. 








47 


This coordination routine utilizes fault currents and load 
Currents in achieving protective-device coordination. The load cur- 
rents determine the overload or pickup setting in each device. The 
coordination problem is then in two parts. 

Devices must be coordinated for overload conditions. This 
provides the most important coordination step for fuses in that they 
contain no variable fault-protection settings as do relays or molded- 
case circuit breakers. Also, another reason is that based upon only 
fault currents and load-side device settings, it is impossible to 
satisfactorily achieve coordination. 

This inability to achieve coordination enters when setting a 
device that is to protect a bus with other elements attached to that 
bus as loads. Tae laying a line curve next to a load curve with the 
highest time value at any current value will provide an erroneous 
answer. This device must be set by determining its overload current 
value first since this current is the sum of all currents in the load 
elements that are attached to that bus. Once this point has been 
established, each load device curve must be polled in turn to insure 
that the device setting responds to a fault in a load element slower 
than the load element line side device. The device needs to be 
coordinated up to the maximum fault current that the device will see. 
It is possible for this device to see a smaller fault current than that 
in a load element. 

Coordinating a device located at a line-side bus is much simpler 
Since the pickup or overload setting is the same as the load-side 


device. The line-side device must be coordinated up to the maximum 
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M@eeescurrent that will occur at the bus where its load-side device 
ms located. A higher fault current will occur the closer one gets to 
the line-side device; however, this higher current value will only cause 
the line-side device protecting the element to respond quicker. The 
method of coordination will be to first find a load point based ona 
load flow, and second coordinate the device at any current value up 
to some maximum fault current such that the line-side device has a 
fenger response time than its load-side device. 

pbanee the program coordinates three different types of devices, 
Mine different possibilities can occur. Each device may be coordinated 
with one of its own type or one of the other two. 

The program structure lends itself well to the use of subrou- 
tines. Two of these are SELECT and SETDEV. SELECT ore a device 


and sets its own load point, and SETDEV coordinates it. 


perce selection and Coordination 


Due to the large number of fuse sizes available, look-up tables 
are not as efficient as using polynomial functions to select devices. 
Once a load current is known, this value is substituted into one func- 
tion providing the proper device number. This device number is next 
E@eetituted into two other functions providing the OFSET and OFSET l 
values referred to in Chapter 2. Since there are two types of fuses 
in the program, current-limiting and solid-material boric-acid power 
fuses, there are six functions used in selecting fuses. These all have 


the following form: 
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Y= f a x0 pal) 


X will be the logarithm of overload current or device number. Table 7.1 
ests the various coefficients aa 

iece tie fuse 15 Selected, 1f must be coordinated. The program 
Meesethe values of OFSET and OFSET 1 with the applicable fuse curve to 
solve for a time value for a set-current value. The program searches 
for the load-side device to find the highest time value to that current 
value. If the time value is larger for the size, it is properly coor- 
feted. If not, it increases the fuse size by one, solves for OFSET 
and OFSET 1, and tries again continuing until it has found the proper 


rose size. 


Relay Selection and Coordination 


Relays provide the most difficult coordination problem. First 
aeeurns ratio for the current transformer (CT) must be chosen. The 
computer selects the proper tap setting to prevent saturation. It 
does this by always selecting the highest turns possible while 
selecting the tap setting on the relay as well [12]. The computation 
uses a 600/5 CT [13] with tap settings and CO relays [14] with tap 
settings. This combination of tap setting for the CT and the relay 


combined with the load current establish the pickup point: 


Pickup =Log [(Turns Ratio CT x Turns Ratio Relay) | (Faz) 





mage 7.1. Fuse-selection coefficients. 
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pemcestelay Curves are in units of multiples of pickup this 
pickup value will be used often. Only the relay pickup has been set, 
its time dial must be coordinated. Table 7.2 lists the relay-curve 
Secrricients. 

Due to the inverse time shape of a relay curve, many points on 
the curve must be checked with its load-side device to insure that 
proper coordination exists. The manner in which the family of relay 
curves are represented simplifies the problem somewhat. A relay 
Merarine time of 0.4 seconds is assumed throughout. At each current 
value, the load-side-device operating time is summed with the relay 
Operating time. If this time is greater than the line-device time, the 
difference is the relay time-delay setting. Thus, when all currents 
up to the fault currents are checked, all points on the line-device 
curve are at least 0.4-seconds higher than all points on its load-side 


device. 


Molded Case Circuit Breaker 


The molded-case circuit breaker (MCCB) with its different 
characteristics provides a different set of problems. The MCCB curve 
cannot be represented by one equation. This device has one set of 
curves for the thermal element plus another for its magnetic element. 
The MCCB used in the program is a General Electric type, K-1200. 
Although there is some difference in the curve shape for the current 
ratings, 300-1200 A, they are similar enough to be represented as one 
type. Each MCCB has one curve corresponding to minimum total clearing 


time and a maximum total clearing time. 
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The magnetic portion of the time-current curve cannot be 
represented adequately by a polynomial function, but is represented 


by: 
ie Ctmrenk = Magnetic Setting 
Minimum Time = 0, Maximum Time = .025 seconds. 


As will be seen later, .025-seconds response time is quite fast and 
uSually does not present a problem during later coordination. 

The thermal elements of molded-case circuit breakers are 
susceptible to temperature change so the operating temperature must 
be entered. The temperature merely shifts the curve to the left or 
right so it does not affect the curve shape. This offset is added 
or subtracted from the multiple of current rating when using the thermal 
curve polynomial functions. 

The problem of coordinating the thermal elements is the Same as 
that for fuses. The magnetic element coordination is much simpler in 
that the magnetic element is set separately. Its load-side device 
operating time is checked at a current value. If the operating time 
meeepreater, the magnetic element setting is increased. The end of the 
curve for the magnetic element is determined by the maximum fault cur- 
rent at the bus. The total decoupling of the overload from the fault 
device of a MCCB makes coordination much simpler than the fuse or relay. 
The more or less proportional-time shape of the curve will ensure that 


it will be properly coordinated once its proper current rating has been 
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established. Table 7.3 contains the polynomial coefficients used in 
the program. 
A basic flow chart (Figures 7.2 and 7.3) follows with a complete 


program listing in Appendix B. 


Summary 


In this chapter, the coordination problem was solved to allow 
Pases, molded-case circuit breakers, and relays to function so as to 
Meedeea fault in the area of the fault before interrupting the rest of 
the distribution system. There is enough flexibility in the program 
to allow device settings to be changed and checked using the ee 
PLOTD.F4. An example problem would be beneficial at this point to 
show how to run the programs. In the next chapter,a sample distribu- 


tion system will be coordinated using the coordination programs. 
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Bepure 7.2. Coordination flow chart. 
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Figure 7.3. Device-graphing flow chart. 
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8. OPERATING THE PROGRAMS 


General 


imechis chapter the capabilities, limitations, and use of 
the programs will be described. The example system of Figure 8.1 


will be coordinated using these programs. 


One-Line Diagram 


INPUT.F4&. All data is: entered via the keyboard in the format 
detailed in Chapter 3. Data entry is a straightforward process 
exe ept aoe entering transformer data. The program is written for a 
ey-y  cOnfieuration. If one or both sides are "4" configured, then a 
large number must be entered for the transformer neutral impedance 
"(ZN)" on the bus connected to the A side. After the program is 
Bempleted, the zero sequence value for the transformer element con- 
taining the transformer may be assigned an arbitrary high value if the 
transformer wiring is such that there can be no ground current through 
Maeselement. This is done by changing its value in disk file TEST.DAT 
(see Appendix C for the file format). 

When entering the loads the program asks which type of load 


Bombe entered. The various loads are described as: 








300 hp SO hp 


mieure 8.1. 


One-line diagram of sample system. 
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fe — Synchronous motor 
HP1 - Assynchronous motor 
Veuea—eein KVA 


GURS— in amperes 


S - No more loads 


Load Flow 


INPUT.F4. There are no entries from keyboard for this program. 


Printing DSK.DAT will provide a load-flow analysis in per-unit values. 


Pettit Oon Of Fault Currents 


Zevsoer4. In addition to forming Pea this program allows 


the input for mutual coupling and for zero-sequence impedance for loads. 
Mutual coupling tends to increase fault currents in the mutually coupled 
elements and small neutral-to-ground impedances can create higher line- 


to-ground fault currents than three-phase faults. 


FALT.F4. Besides the solution of fault currents described in 
Mraprer 6, FALT.F4 also allows the input of fault impedances. This can 
provide an interesting study of the system because the user can go 
back and place different fault impedances into the fault analysis to 


Peewitr the protective devices will operate as desired. 





OZ 


Device Coordination and Blopting 


COORD.F4. The program forms a two-dimensional complex matrix 
@atled ADEY that contains all the information required to plot all the 
devices in the system plus all actual settings that are to be made on 
mieedevices. An example row 3 in the matrix is 3,(3,4), (-4,2), 


mimes), (1.03,1.26), (3,667,0), (.67,0), (20,0) in which: 


imener 3, 3 is the row location of the matrix and corresponds 


to the element number. 


PeetoOr (3 4).0906 ise the first column llocation and is the bus 


numbers of the element. 


femeber (452),-4 is a relay type #4; 2 is a current limiting 


fuse type #2. 


eboreie7o.205), 1226 1s the logarithm of the pickup point; 
203 contains the tap information. Two is the tap of the 


Cl wancdO3.1S5 the tap of the relay. 


omer ont.260), 1.03 1sethe logarithm of the current for 


the minimum fuse open; 1.26 is the typical fuse open point. 


rmeommcs700/,0), 3.00/ 1s the magmetic setting of the relay; 


O since the fuse had no magnetic setting. 


i OLatGroy sO) ene oy 81s ashe time-daal setting for the relay. 
The actual time dial setting must be read from the manu- 
facturer's family of relay curves using the formula below 


Ae nOo maximum fault current. 





63 


time-dial-setting found = (found-time curve) - 


(time-dial 1/2 setting) 


Beeeeror (20,0), 20 is the turns ratio of the current trans- 


former. 


Device numbers the programs asks for correspond to the 


devices below: 


- Current limiting fuse 


IL 

2 - Boric fuse high voltage only 
3 - Molded-case circuit breaker 
i 


- Overcurrent relay 


In device 4 the program will ask later what particular type 
@rerelay to use from Table 8.1. 

Upon completion of coordination, a device listing is typed 
on the keyboard, and a printout of the same is available by typing: 
PRaINt*. LPT. 

When performing the plot, having a device listing as well asa 


matt current listing is helpful. 


PLOTD.F4. This program plots the devices that were coordinated 
above. By using a graph, it is possible to visually check coordination. 
The program arrangement allows the user to change any device setting 
@eesize in ADEV.DAT to tailor the coordination to his particular require- 


ments (see Appendix C for the file format). The program uses the lowest 
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Table 8.1. OQOvercurrent-relay types. 


Program Type # 


7 


Relay Type 


CO-2 


CO-3 


CO-4 


wO=> 


CO-6 


CO-11 


Relay type numbers are listed as negative numbers by the program. 
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voltage in the system as a base for plotting currents. The user can 
Select any plot offset. This has the effect of plotting only the 
area of concern. In this way a very detailed graph can be obtained. 

Tables 8.1-8.4 provide a listing of program device type numbers 
to the actual device type and size. 

Table 8.5 details what is contained in the various disk loca- 
tions created by the programs. 

To coordinate the devices in Figure 8.1, the procedure begins 
by entering element data using INPUT.F4. The program asks if data is 
to be read from disk or from the keyboard. The first element must be 
entered using the keyboard. 

hes proper procedure is to group the elements to be entered 
by wire or cable type as the program only handles one type of wire or 
cable at a time. When completed with one type, answer zero for all 
other questions and the computer writes the element data that it has 
computed on disk in the file RELA.DAT and stops. The new wire type 
polynomial coefficients are copied into RES.DAT and the program is 
ready to be run for the new wire type. 

The cable coefficients used in the test problem are for mine 
mower Cable [5]. The coefficients are listed in Table 8.6. 

After selecting keyboard, the base power and volts are entered. 
The computer then informs the user to begin element data entry. The 
aaeaetormat (Chapter 3) is listed to aid the user. For example, the 
element from bus one to two of Figure 8.1 is entered as follows: 
1/2/3000/6. For the element containing the 1 MVA transformer the follow- 


mews input: 2/3/35/6/.001.005/1/15/4.1. The bus base voltages -are 





feple S.2. Current—limiting fuses. 


Program Type # 


mo CON DW NO &F W PDP FF 


PR Pe RB Be 
wm & wb NH FE Oo 


Limiting Fuses 


aE 

7E 
LOE 
LSE 
20E 
Zot 
DUE 
4OE 
5O0OE 
65E 
80E 
LOOE 
2 OE 
I50E 
200E 
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mapre &.3. Boric fuses. 


Program Type # 


16 
ey 
18 
io 
20 
21 
ee 
Z3 
24 
Z) 
26 
Z/ 
Ze 
Zo 
30 
oa 
32 
33 
34 


Size 


NE 
205 
25E 
30E 
40E 
SOE 
65E 
80E 
LOOE 
er 
L50E 
eo 
200E 
2505 
300E 
4O00E 
2-250E 
20005 
2-400E 
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Table 8.4. 


Program # 


34 
35 
36 
5 / 
38 
59 
40 
41 
42 


43 


Molded-case circuit breakers. 


Size 


300 


350 


400 


450 


500 


600 


700 


800 


1000 


1200 


amps 


amps 


amps 


amps 


amps 


amps 


amps 


amps 


amps 


amps 
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meme G.5. Disc files. 


reir. F4: 


DEV. F4: 
INPUT. F4 


LDFLO. F4: 


mauS.F4: 
MALT .F4: 


COORD. F4: 
PLOTD.F4: 


Fea. DAT : 


iol. DAT: 
RELA. DAT: 
XFOR. DAT: 
EGCUR. DAT: 


DEV. DAT: 
DSK. XFR: 
ok. DAT: 


Pes .DAT: 
ZBUS. FOR: 


PeUr.co: 


PALT.DAT: 
eis . DAT: 
eCUR. DAT: 
CUR. DAT: 
ADEV.DAT: 
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Curve-smoothing program 


Numonics input program 


Element data-input program 


Load-flow program 


HOrmation of ZBUS 
Solution of fault-current 


Coordinate device 


O12 
program 


program 


program 


Plot device curves 


Contains 
Contains 
Contains 
Contains 
Contains 
Contains 
Contains 
Contains 
Contains 
Contains 
Contains 
Contains 
Contains 
Contains 
Contains 


Contains 


X/R line coefficients 

element and load data in per unit 
element and bus voltage data 
transformer inrush and withstand values 
load currents 

relay coefficients 

load-flow currents in P.U. 

load-flow analysis data 


ee a 


zBust? before loads are added 
mutual coupling data 

per-unit fault currents 

per-unit fault voltages 

real values of load-flow currents 


beal Vvalnes “of fault currents 


coordinated-device listing 





Table 8.6. Mine power cable polynomial coefficients. 


R-resistance h- reactance 
ay 027500905267 34D+01 0.6731452381407621D-01 
ay = .1108155100578006D+02 - .1408237729509665D-01 
ay 0.6290760668063740D+01 - .1721580927352306D-02 
a. ae Lom 2 1.428156510Z2D+01 -6845126542245339D-03 
O.1604746751660286D+00 Oe 





Volt 


input next with the voltage in volts. Upon completion of entry of all 
bus base voltages, the computer types the computed element value in 

per unit. The load data is entered next. The program types the data 
entry format for the user. The 50 hp motor is entered as: 8/50/.85. 
The only other data to be entered are mutual coupling and zero sequence 
load impedance in ZBUS.F4 and fault impedance in FALT.F4. The sample 
problem uses no coupling and a fault impedance of zero. COORD.F4 is 
run next. The device type desired is entered when asked for by the 


program. A sample data input is given in Appendix D. 








The following 


eurrents , 


Bus number MOad secur rent 
4 84.0396680 
4 84.0396680 
ii 79356613400 
o Wow od 3 55/0 


The following is the current flow 


iomeaeelenGtaO Lr swodd. lLocatiLons 


and 


in the systen 


maaer full load conditions as found from the 


load flow analysis. 


Bus Pp Ou usS Q 


HS MW WN ee 
COUSIN Une to 


Current 2c P 


34.2.2226440 
Boe 22 2b o30 
Bow sy 257.0 
ecnlus 2.27 i) 
So aly ee ey) 
313545357 U0 
51-6573760 


Current ate 0 


34.2226440 
P25 C457 00 
88 .0372570 
See 22 1) 
316570129 506 
SS 6357 00 
Sle 650 37 Oe 


The following is the computed fault currents for a 


feprase fault and a line to ground fault. 


Flt Bs P Element 
2 ] 
Z 2 
2 3 
Z. 4 
Z 5 
2 6 
Z yi 
3 ] 
3 Z 
3 3 


Bee hase 


OS OUew 2 92 O00 
7 Gules oo 5 7 OU 
O77 ~8lyO200 
280-.0413400 
Zr oe CO 

222 Ono WOO 
383./966800 

2 eee one OG) 

eZ Gy a OO 
846.6523 100 


Ln toe Gne 


7 Opn 27 80010 
teow? 685400 
Deo 2 | 200 
eg Te aleve eNO) 8) 
ore 5 5 7 00 

1496.4606000 
27 O93 ob OOO 
wee O20 7 00 
826.2558200 
re 32.00 





Flt Bs P 


PPh HWW W Ww 


—_ 


Come, Co Co Co ee CoS SS tO Ss OO Gs Sn nt Sr rt rn rr i 


Element 


SOME WNH BSE NSDUH WN ENDO PEP WH BANDON WD BND LW Ph HR ~AIKU 


3 Phase 
S49 o72 695 00 
Sadie OO 1 UO 

273059159000 
479.3659400 
iO 5 JOO 
PON Sa40000 
4003.09338000 
JAS Sie SPA 850 BONG, 
295923352800 
Pee all> O70 00 
391.0159300 
Some 773500 
6831.8899900 
593 .6094700 
B805- 2500000 
553. 1902700 
2009 «5142000 
463.991 2400 
791.2624200 
PAs Seen) 
Do @0o00 OOO 
3416.9009000 
S4G.o 7 I9SC00 
2559 5040000 
aye 4 4500 
Soe oO ly oll 
BOO 2 4000 
Pas OOO) 
151375939000 
Sass. 359 50100 


14273.8630000 


Zag or) Onl ae 
eo Gree O11 96:00 
Por 082 1.90.0 
ON es SONONS O10 
67642499600 
G7 O24 O55 00 
Ole 7S 00) 
Oy Cee,» 10010 


Lo £0 Crd 
206604350 00 
20iGea2 07 aod 
1648.7648000 
3075662200 
469.6438400 
Olio. 4 90 0 
27 NO. 2o0y 0.00 
121.9024600 
M2507 2000 
972./7699300 
La OL e0'0 
306653295700 
SUG a2 17.00 
WAG Os Waele 8) ONE, 
Wage 17 25000 
LOS. 2970700 
DU lez soo) © 
183.8482000 
148.9287400 
14854773900 
125.4932500 
668.5896000 
569.4760500 
517 .3090900 
le7 ee ell OO 
WO3 nel 31-4010 
103. 1008700 
70-4041400 
445.7008900 
443.97/70900 
4235.3970000 
G4 26075200 
Zoe 7 Oa 
23«5401980 
pray OO 0 
POZe 2265000 
100. 13874100 
84.07806/70 
Oso 3405000 
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Using the device listing of ADEV.DAT and the plotted curves, 
tt 1s Peecibie to obtain all device settings. 

Daevieeenumpensone 1c a molded=case circuit breaker. Its size 
is 300 A and its magnetic trip is on minimum setting. 

Device number 10 is a current-limiting fuse number 10. 

Device number 12 is a solid-material high-voltage boric fuse 
number 24. 


The relay data follows: 


Device ipe Time 

Number Number Cre Tap Dial Magnetic 
3 i. Sack 0) iy 2 -—= 
4 I 3 8 8 6 

3 3 3 14 5 6 

6 4 3 3 4 6 

i 2 3 4 6 6 

8 5 3 5 6 6 

0 6 3 6 EO eye 
2 1 5 13 » Dios 
ia a 3 4 10 6 
3 7 3 iu 2 6 

14 Z 3 8 ez 6 


The magnetic trip of device number 9 had to be reduced after 
coordination to allow the device 12 to set properly to provide pro- 
tection for the 1 MVA transformer. The graphs that follow een the 
devices plotted to the lowest voltage in the system, 440 V. All the 
graphs were computer generated. The first graph shows how a one-line 


computer diagram would be placed to aid in understanding the graphs. 
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Summary 


The sample distribution system of Figure 8.1 was coordinated 
using the coordination programs. Using the device-plotting program 
PLOTD.F4 this coordination was checked to see if the computer generated 
coordination was satisfactory to the user. 

During the research for this thesis,many observations were made 
and many features could be added to these programs in the future. This 


will be discussed in the chapter that follows. 
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CONG LUSLONS 


In this thesis the author set out to represent a one-line 
distribution system in a per-unit system and to solve for three- 
phase and line-to-ground faults and to perform a load-flow analysis. 
The various protective devices-response curves needed to be represented 
by polynomial functions. These requirements were fulfilled using the 
methods detailed in Chapters 2-7. 

The last task was to develop a computer procedure to coor- 
dinate a distribution system's protective devices. Fuses, molded- 
case circuit breakers, and relays were coordinated and their coordi- 
Hation graphs were plotted. There are many extensions for further 
work that could be done. These extensions will be discussed later. 

While working with the various programs, several observations 
have proved useful when performing coordination. Coordination should 
not be an afterthought in designing a distribution system. it is 
easy to design a system that is impossible to coordinate correctly. 

A radial system is easier to coordinate. Moving from a load toward 
the source, the time-current curves shift to higher current-—time 
responses so that when finally arriving at a transformer, providing 
proper transformer protection and coordination simultaneously may be 
impossible. 

Coordination of the thermal-tripping elements of molded case 
Circuit breakers is very difficult due to the difference between the 


minimum and maximum tripping time curves. Coordination of these 
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devices is best done manually with some overlap of the thermal time 
curves. The magnetic elements could then be adjusted to provide 
proper fault coordination. Those manually adjusted devices can be 
added to ADEV.DAT and be checked using PLOT.F4. 

In systems with different voltages, the transformers should be 
of different sizes so that the smaller sizes are on the load side. 
This will allow proper coordination and protection. Relays are 
preferable to other types of protection on the high-voltage parts 


of the system in that their response can be controlled. 


Recommendations for Further Work 


1. Using an analog computer and the programs presented here, 
an entire distribution and generation system could be 
looked at in detail. A transient-stability study would 


show the effect of various device settings. 


othe INPUT E4 program could be enlarged to automatically 


account for different transformer configurations. 


3. Additional devices could be added easily to the program. 
The program allows up to twenty different relay types. The 
additional coefficient data is merely entered in disk 
DEV.DAT. Additional molded-case circuit breakers would be 
Sirgittlveanore ditticult im that if would require some program 


mMmoditicativon. 
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The plotting program could be enlarged to include drawing 
a one-line diagram in the space provided and adding other 


data directly on the graph. 


The program could be expanded for the case of having one 


device in an element. 
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APPENDIX A 


CURVE-SMOOTHING PROGRAMS 








C 
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iGgeewbel tera CREATES FILES FROM DATA 
USING NUMONICS ARM 


KKEAKAKKKHKKEKKKKKKKKKAKKKEKEKRKKRKKKKREKARKRKKKEKKAEKKKKKRKRKEKSE 


1 
3 


5 


10 


15 


20 


30 
3 


100 


Lone 


200 


210 


AY 


TYPE 5 

ROR ALGoeENPUT DEVICE NR...” ) 

ACCEPT 8,IDEV 

FORMAT(1) 

OPENCUNIT=IDEV) 

mye 10 

FORMAT(’ INPUT X-AXIS DIMENSION, Y-AXIS 
Bete SHON ff * s6, 7) 

MO@G@EPE 15% Ax. VAX 

FORMAT(2F) 

TYPE 20 

Por wean LNPIT  DENGTH PER LOGIO:X;Y..”) 

PGGHPs 354 XL4 VL 

FORMAT( 2F) 

KAX=1N0*XAX/XL 

VAx= LOOSYAX/ YL 

Cheer nNcgO. DEVICE’ , NUM’ ,“MODE’ ,’ IMAGE’, 
bene Grios — SEOIN ,° BUFFERS’ ,4) 

GCMMieN URES C30, 4) 

GC ie bock( 1, )) 

Mere OG. END=200) IDATA 

eae iy  bOooOO00 ye AND.'37777 
EV pee AND. "37777 

Pemacs—G@1 pATA.AND.''740000)/'"40000 

X=1X/XKAX 

yy AX 

Heir oGipey, 150)1LFLAGS, X,Y 
PORMAT( 1. 2F ) 

CO TO” 100 

CALL NUMCLO 

CMne mweoCk (1,1) 

Cros Cunt 1 DEV) 

PYPE 210 

DORM An = CONTINUE?..,° S$) 

OCP i 220.1 C 

FORMATCL) 

Mere ehong)CcO TO 1 

STOP 

END 








_— = ee 
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Kx KK KKK KKK KKK KKEKKEKEKK KKK KKEKEKKKKKKKKRKKKK KK 
C CURVE FITTING PROGRAM:DEV.F4 
BKK KKK KKK KEKE KEKKKEEKEKKEKK KER KKKKREKKEKKEKEKEEKKEKH 
INTEGER L(7),MM(7) 
aie © CLO), YS( 100) 
eres oeeClom ),C(7.7),U(7),A(0/6),6G(7,7),ERROR(100) 
TYPE 100 
100 FORMAT(1X,’CHOOSE DIMENSION OF FIT ’$) 
RECEPT 110,ND 
110 FORMAT(1) 
Tybee 20 
120 BORMAPCIX, CHOOSE NO. OF POINTS oe) 
ACCEPT 110,M 
TYPE 130 
0 MORIMATC/° CHOOSE INPUT FILES’/*% ..., oe 
MmeCEPT 110,1UNIT} 
DO 140 I=1,M 
MeAmiCrIuNtrl, 135)J,xX(1),YS8(1) 


3D FORMAT ( 1, 2F ) 
ioe y= ALOGLOCYS( 1) ) 
140 CONTINUE 


DO 10 K=1,M 
bO 10 I=1,ND 
WII )=UC1)+YS(K)*X(K)**(CI-1) 
DO 10 J=1,ND 7 
Come =e Te )4+ x OK) ** (1+3=2) 
10 CONTINUE 
CALL DMINV(G,7,ND) 
DO 20 I=0,ND-1 
DO 20 J=1,ND 
NOt AC 1)+G (141 ,3)*UCJ) 
20 CONTINUE 
Hows0 l=), M 
DO 30 J=0,ND-1 
YMCL)=YMCI)+AC J) *X(1) ** I 


30 CONTINUE 
DO 40 I=1,M 
40 SQUERR=SQUERR+(YS(1)-YM(1))**2 
TYPE 160,SQUERR 
160 FORMAT(’ THE INTERGRAL-SQUARED ERROR IS: a Es) 


OPENCUNIT=40, FILE=’A.DAT’ ) 
WRITE(40,170)A 
CLOSE(UNIT=40) 


170 FORMAT(D) 
Terre 99 -1yN rT 
99 POV nGae lit 1) 
Cue = 1M 
50 ERROR(I)=DABS((YMCL)-YS(1))/YSC1)) 


OPEN(UNIT=41, FILE=’YINOUT.DAT’) 
Wome o yC% Cdl ),YS( 1), YNCI),ERROR(1),1=1,™) 
180 FORKMAT(2F, 2D) 





c 


nO oe euai f= ]\) 
Soe 
END 


ewer ys SUBROULTIRNE TO INVERT A REAL SQUARE 


MATRIX DOUBLE 


C 


a 


10 
eS 


70 


D0 
5 


38 


40 


45 


50 
a) 


60 
6 2 
65 


70 
nD 


PRECISION 

SUBROUTINE DMINV(A, MDIM,N) 
REAL*8 A(MDIM,N) 

MILEGEREL« 100),MC 100) 


COVER EDeEROM SSP ROUTNE MINV 


Howe he PRECISION BIGA, HOLD 
poecoe—-1, N 

men) —% 

MC K)=K 

BIGA=AC(K,K) 

Demo. Jk. Ny 

Dom 0 1=k,N 

Mae nSCBIGA)=DABSCACI,J))) 15,20,20 
BIGA=A(1,J) 

eC) = 1 

Hood 

CONTINUE 

ee 1c) 

mnCd—-) 35,35,25 
Demo. r=, i 
HOLD=-A(K,1I) 

hs )=ACIJ,1) 

hol) =HOLD 

I=M(XK) 

reac) 4545, 38 
Doro sJ=1.N 
HOLD=-A(J,K) 

Glen) —70Gd, 2) 
Gee) = 110.1) 

DOmoo) t-l,N 

Photek) 50,55. 50 
aie) =o) / (— B81 GA) 
CONTINUE 

Domes Lol N 
HOLD=A(1,K) 

MO) 1S ST 

IF(I-K) 60,65,60 

Mtv Kk) G26 9 6 2 
NGsero— HOm peak At I yd) 
CONTINUE 

Dome 5 = 12.3 

IF(J-K) 70,75,70 

ACK, J)=ACK,J)/BIGA 
CONTINUE 
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100 


105 


108 


iO 


wae 


125 


150 


Me v=— (hen 0000 )/ BIGA 
CONTINUE 

K=N 

K=K-] 

Mack ei sO, 150.105 
P= kK) 

Pcl) 120,120,108 
DO 110 J=1,N 
HOLD=ACJ,K) 

med J==AC J, 1) 
A(J,1)=HOLD 

J=MCK) 

Pci= kor 100,100,125 
Por asOel=1,N 
HOLD=ACK,1) 

ACK, I)=-ACJ,1) 

Sel t)=HOLD 

GO TO 100 

RETURN 

END 
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ie ee OO" ei 
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APPENDIX B 


COORDINATION AND PLOTTING PROGRAMS 








a1 


PRIOR KK RAK KKKKKKKRREKKKEKRKKKKKKKRKKREKKKEKKRKKRKKERKKKKEKKEKE 


C INPUT.IT INPUTS DATA FROM A TERMINAL, PLACES IN 
C THE PROPER FORMAT FOR OTHER PROGRAMS 
KHKEKEKKEKKKKKKKKKKKKRKEKEKKEKKRKRRKRKKKKRKKKKKKKKEKRKKKKKKKKEEEKER 
C 
C 
DIMENSION 
meet, 2), PU(30,30,1),RES(30,30,2),E(30),XL(30, 30, 2) 
DIMENSION 


mee), ZER( 30, 30,2),RC(7),XC(7),2Z6(30),xXDC(30) 
PoMroexeetl., TPU,PU,S,S1,PU1,PZ,CUR,ZER,ZG 
REAL*§ A(7),AAC7) 


REAL L 
KEEKKKKAKKKEKKEKKRKKKKRKE KK RRR UEKK KKK KK RE KE KKK KEKE RKEKKEKEKEEKE 
C ieee er hOGRAM HAS ALREADY BEEN RUN ONCE ALL DATA 
C (eee NOT BE ENTERED AGAIN 
KAKKKEKKKKK KKK KK KKK RK KK KK ARK KEKE KKK ARK KAKA KKK RE KEKAEKKKKKEKERE 

ies 3 
B FORMAT(” DO YOU WANT TO READ DATA FROM A FILE OR 


FROM KYBRD?°/ 
occ > EORSINPUT FROM TTY,! FOR READ... ) 
NCCMeET 1), 1 


WiGmeNE.0)GO TO 400 
SHRM KRKRKKEKKR KKK KKK KARE KEKRK KKK KEKKEREKRERKKKRK KKK KKEEMKEEK 


C INPUT OF DATA ON LINES 
Deer or KKK KKK KKK KKK KRKKREKKKRKKKRKRKRKKKRKEKKRRKEKKRKEK KKK 
TYPE 5 
5 FORMAT(2X,’ENTER THE PROBLEM BASE IN MVA‘’) 
ACCEPT 10,ZBASE 
TYPE 6 
6 FORMAT(2X,’ENTER THE PROBLEM BASE IN KV’) 
ACCEPT 10,EBASE 
10 FORMAT(F) 
fel FORMAT(I) 
COmUG) 14 
ee KKK KKKKKKKKKKAEKRKEKKEKEKEKKKKEKAEKEKKE KKK KEKSE 
C PROGRAM RETURNS HERE AFTER READING DATA AT 400 
Meee RK KKK KKK KRAKKKEKEKKKKRR KKK KKK EKEEKAKRKEKKEEHRK 
a ier 13 
3 FORMAT(’ OLD DATA IS ENTERED,READY FOR NEW LINE 
ATA...’ ) 
Nein 1 Ld 
14 TYPE 15 
US FORMAT(/,2X,’ENTER ELEMENT DATA AS 


FOLLOWS: BUS( )TO( ) (LENGTH) 
ese smh / , 2X, ° CTRANSFORMER 
ME DENCECP.U.))(TRANSFORMER SIZE 
Dien eeop 2k. (TRANSFORHER WIGH VOLTAGE IN KV) 
SP‘ ’ / ’ 
Semester sboOl VOUTACE IN KV]«.-.,’/) 
KRKEKKKKRKEKREKKER KKK RK RK KKK KKK KKEKK KKK KKKKEKKKKEKHK 





C 
C 


Gite teel eerie o> CONTAIN COEFFICIENT 


PROGRAM: CURV.F4 


KREEKCKKRAEKKKKKKEKKKKKKKUKKEKKKEKEKREKKEKREKKKKKEKKEKERKEKKKKKEKKEKE 


Z 1 


26 


a7 


31 


32 


33 


34 


KKEAKKAKEKKKKEKKKKAKKKEKEKRKEKRKKEKKEKEKEKKEKEKKKEKKEAKEKKEKREKEAKEKE 


C 


RHEEEKKKEKKEKKKEKEKKEKEKKKKEKKEKRKEKKKKKEKKEKKEKRKEKKKEAKEKEKEKKRKARKRKEKK 


5) 


KRKEKEKEKEKKKEKEKEKEKKKEK KKK KEKE KKK KKEEKKK KKK KKEKEKEKKKKRKEKEKE 


Ok veer Ok bE SRESTSEANCE AND REACTATCE 
KK KKK KKK KKK A KAR RK KKK KK KKK KEK KKK KK KEK KKK EK KEK KEKK KKK KKK K 


C 


OPEN(UNIT=40, FILE=’ RES.DAT’ ) 
OPEN(UNIT=35, FILE=’ XFOR. DAT’ ) 
READ(40,290)(AC(1),1=1,7) 
READ(40, 20)(AAC(I),I=1,7) 
FORMAT(D) 

Moy 21 t=). 7 

RCC 1I)=ACT1) 

XC(1)=AACL) 

CONTINUE 

CLOSE(UNIT=40) 

CLOSE(UNIT=41) 

Me@eet 24-120, %Li.WIR.11,761.TVH, TVL 
FORMAT(21,7F) 

MEGIENEQ.0) CO TO 110 

TYPE 26 

JJ=1 


FORMAT( ’ ewe i Ns IeenO™ Bievcey 5 5) 


MOCEPT 27, J 

FORMAT(L) 
IF(JJ.EQ.0)JJ=1 
Me@mreieet ds i= LL 
iaevERehO.140)GO TO 31 
IF(WIR.~EQ.130)GO TO 32 
BGR. hOsl20)C0 TO 33 
Per IReEO.110)GO TO 34 
SA=WIR 

SB=WIR 

ieGstebn~37)GO TO 35 
GO TO 40 

SA=-3 

Core 35 

SA=-2 

COmEO: 35 

SA=-] 

O20. 35 

SA=0 

Cor TO. 35 


CONVERSION FROM AWG TO MC? 


§C=1.1229**(36-SA) 
SD=.005*SC 

oo Dest 2 
SB=1000*SE 


DAT AVEO G OP BELOW 
ity Sh eSCkhEATED USING THE CURVE FITTING 


ee eee eee eee OL," Oe eee, ee, eee 


a a 








os 


40 
meet, 3, JJ )=2C( 1 )+RC(2)* (ALOG10(SB) )+RC(3)*(ALOGIO(SB 
ee )+ 


memes )<(ALOGIOCSB) **3)+RCC5)*(CALOG]IOCSB)**4)+RC(6)*(A 
moe loc SB) 
2**5)+RC(7)*(ALOG1O(SB)**6) 


Serr, J )=XC(C 1 )+XC(2)*(ALOCIOCSB) )+XC(3)*(ALOGIOCSB) 
Bare? ) + 


meaty )*(ALOGIO(CSB) **3 )+XC(5)*(ALOGIO(CSB)**4 )+RC(C6)*C(A 
moelocssB) 
2**5)4+XC(7)* (CALOG1IO(SB)**6) 
eer ot ie) ) 4 SA, SB, XLC1L,J,JJ) 
DO Oa Ce ee NE RES. ~,3F,/,° LINE REACT...’ ,F) 
imeernSG(t1).6b9.0)GO TO 85 
mel so —11*(ZGbASE/TS1)*(C(TVH/ EBASE) **2) 
mye 5.5, L407 
55 BORMATC  TNPUWT TRANSFORMER ZN IN P.U. REFER TO BUS 
Paepe , (BUS 
ele ) OV (BUS ~,13,°)..,°) 
MeCEPT 60,S1,PU1 
60 FORMAT(4F) 
TV=TVH*1E3 
ZG(L)=S1*(ZBASE/TS1)*(CEBASE/TVH) **2) 
Ae =PUl-(CZBASE/ TS1)*((CEBASE/TVL)**2) 
Cm B=TSI*lEG/TV 
TINRSH=12*CURB 
TWSTI=25* CURB 
Desi 2=14 74+ CuRB 
Hie Cos, Ol. !.1INRSH, TWST1,TWST2,TS1 


70 FORMAT(21,4F) 
85 MiOntANsGR.1) CO TO 90 
IMAX=I 
90 Maenii<sered) CO TO 00 
JMAX=J 
100 Ae ee = CMPEX(0.0,XL(1,J3,JJ)) 
iho? 2 ER(L,J,05),S9S1 
102 FORMAT(3F) 
Gomero. 22 
m0 OPENCUNIT=31,FILE=’ RELA.DAT’ ) 


Ore GUNIT=32,F1LE= TEST«DAT’ ) 
OPEN(UNIT=33,FILE=’LOCUR.DAT’ ) 


IB=INAX 

COMEO 11S 
lS IB=JMAX 
KAKA KKK KKK KARA KK KE KKK KAKA KKK AKIRA KK AKA KKK KKK KAKKAKK 
C HU OLiAGEo ARE Liput 


RAEKKAKKKKEKKKEEKRAEKKKKEKKKEKKAEKRRAEKRKEEKEAKRR KK RAR KKEKRKEKAEKEKREKAKKEKS 


ee 7 








120 
2 i 
2 2 
22 


126 


128 


30) 


94 


ioe Oat — len 11 

TYPE 119,1 

noe aminrut BUS .13, VOLTAGE..:,°,$) 

eC EP a ThO. EZ 

WCh2e30.0)CO TO 121 

Ext )=E2 

CONTINUE 

TS WAY 

orate DATA 1S ACCEPTED AND A LIST FOLLOWS’ ) 
Deb? 5.2845 5 

POMC er ePROBLEM BASE IN MVA..,° ,F) 
WRITE(31,126)ZBASE, EBASE, 1B 

FORMAT(2F,1) 

E1=E(1)/CEBASE*1E3) 

Petes? .l28)518, 51 

FORMAT(1, 2F) 

WOn135 1=#1,30 

pOm133 J=15.30 

DO 33: JH] 4:2 

PaGnGimaieid ib Oe.) CO TO.133 

NN=NN+1 

Ea eC eee tie ex (il, J, JJ)))*C.0014L(1,3,JJ)) 
PZZ=(ZBASE*1E6)/CEC1)**2) 

PZ=PZ*PZZ | 

ie Pt Pete.) J) 

Ao An Roni id )*P2ZZ*(.001*L(1,J,JJ))+T1 
iat, 130e Ni. bClsJ.IJ),7T1 

Home coe b END 9 9503.5%, LENGTH’ ,F,/, 
11X,’TRANSFORMERPER UNIT IMPEDENCE PLUS LINE 


IMPEDENCE’ , 2F) 


meets) ) iste) .bCl,J,JJ)SRES(L,J,JJ),PUCIL,J,JJ), 
fect, J, JJ) 


131 
132 


3 
3D 


| She 


37 


13S 


3 9 


140 
145 


MeZG( Ly .260)) 

FORMAT( 31, 9F) 

a Cone i win Tl, ZERCL,J,JJ),Z6°1),Z6( 3) 
FORMAT(41, 8F) 

CONTINUE 

CONTINUE 

WRITE ( 312136) 

WRITE( 32,136) 

DO aemoGmomoro 0.0 0.0 0.0 0.0 0.0 0.0 0.0’ ) 
Wem s or. eae 

Vivier i 3 lena bs bi 1) 

FORMAT(I1, 2F) 

CONTINUE 

WWE ( 31,139) 

BORAT ( ” 2020/0 02:0 ) 

Ta) VAs 

Poti Clementine OF BUSSES IS...,°’,1) 

TYPE 150 





7D 


KEEKKKRKEKKKEKKEKRKKKEKRKK KKK KKRKKKKKEKRKREKKKRRKKKREKEKKRKKKKKKKKKEK 


15 0 Bonds toy eet OWOENTER THE LOAD DATA IN THE FOLLOWING 
PPANNER: ... 
lee eeChvomm ChOAD IN KVA)CPWR FACTOR)...OR..-< 
er eeeebinice™) (LOAD DN HP)CPWR FACTOR)....OR....°/ 
3° (BUS )CLOAD IN AMPS)CPWR FACTORY’ ) 
DO 300 N=1,1B 
TYPE 160 
160 FORMAT(” WHAT METHOD WILL YOU USE 
feeVA,HP,OR AMPS??°%) 
ACCEPT 162,MD 


NG 2 FORMAT(A4) 
TYPE 164,MD 
164 RORMAT(2X0AG.° ENTER...’ ) 


PaeubebOn HP1° ) GO TO 170 
PECiMeenOn He.) GO. TO. 170 
PEG Oe AMPS ) GO TO 180 
DEG. 5° ) GO TO 305 
MD=’VA‘ 
ACG dt V6 5... J 4/A5 14 PR AC 

65 FORMAT(1, 2F) 
SZ2= hol) CZBASE* 1E3) 
SZR=SZTZ*PFAC 
PCOS=ACOS(PFAC) 
SAI=SCIZNS INC PCOS) 
Si -CMPEXCSZR, SZL) 
pis * 1h3/ CECI) <SORT( 3.0) ) 
WRITE(33,168)J,AI 


168 FORMAT(I,F) 

ST=MD 

Co. TO, 300 
70 ACCEPT 175,1,PRH,PFAC 
175 FORMAT(1,3F) 


SZTZ=PRH*746 

SZR=SZTZ*PFAC 

PCOS=ACOS(PFAC) 

Pal oA ees PCOS ) 

Sal) enteeaco ZR. SZ1)/ (ZBASE*1E6) 
Puen anO. He )XDPCI)=.17/CABSCS(1)) 
Ml o7 E77 eSORl(3.0)*E(1)) 

rem os, 178)1, Al 


178 FORMAT(I, F) 
ST=MD 
GO TO 300 
180 N@Cnnd soon lL, CURT, PRAC 
185 FORMAT( 1, 3F) 
VREnNCa sels 7 )1, CURT 
187 FORMAT(1, F) 


CURR=CURT*PFAC 
PeOo=NCOSCPr AC) 
CViki=CURT*SINCPCOS ) 








300 


B05 


510 


Bae 2 


Bil 


B20 


50 


55 | 


2) 34 


BD 


400 


411 


450 
456 


470 
472 


457 
458 


900 


96 


He CuPLXC CURR, CURL) 
CUR=CUR* (SQRT(3.0) *EBASE/ZBASE) 
Sol JE BASE=~CONJG(CUR)/SORT(3.0) 
ST=MD 
CONTINUE 
Wom o20.N=1,18 
NN=N 
wee, oO. Nes N)sMD 
MoOmtAtTC BUS ~,13,2F,A4,2F ) 
WRITEC31,312)NN, SCNN) 
PORMATCI,2¥F) 
helo 2s OU SCN). XDCNN) 
HORMAT( L, 3F) 
CONTINUE 
Wielte (31), 330) 
WEEE (32, 3 3:1) 
FORMAT(’ 0 0.0 OQ. 
BORMATC” O 0.0 
eet 33,332) 
mOmiAr(- 0 0.0 ) 
Perte( 35,335) 
MemiAnC 20) 0 0.0 0.07 0.0 0.0 7) 
eo. lo 900 
Coe CUuNTT=36.,FTLE= RELA.~DAT’ ) 
oso ll) 2BASE, EBASE, LB 
IMAX=I1B 
RORMATCZF,1) 
Mee wie oreo Ea lu kL, RES I, PUL; XL2,T1,S1 
PORMAT( 31, 9F) 
Para ho.) CO TO 470 
oer. IJ )=xXL3 
Meo. JI y=RE oO] 
BUC 3, 1 )=Pu) 
MeGin FoI )=XL2 
An led sb =CMPLX(0.0,X1L2 ) 
ZGC1L)=T1 
Ze Ca) = Si 
DiGlaNE.o) GO. LO 450 
Mewes ow472)0.,E¢1) 
BORMAIC £2 2F ) 
Paget sas GO TO “470 
REA 36,456) 1,5 | 
FORMAT(1,2F) 
S(1)=S1 
Pe@eibend) “GO TO 457 
KEWIND 36 
J J=0 
CLOSECUNIT=36) 
Cia Loe | 2 
CLOSE CUNT T=3 1) 
CLOSE CUNT T=3 2) 





oli 


Chose. UNTT=3 5) 


Sel P 
ERD 
PKK IKK KKK KK KER KEKE KEE RIK KEK KEKE RR KK KKK KKK KKK KK KKK KEK KKK 
C porno wbaskA LOAD FLOW PROGRAM 
KKK KIRK KR KK KK KKK KEK KKK KEKE KR KKK IKK KKK KKK KEKE KKK KRKKEKEKEKES 
POMENS LON 


mummeee 3, 33),PS5CH(33) ,QSCH( 33), PCAL(33), 


umes) TP ULEEF( 33) ,ODIFF(33),10(033),8(33),JACC68, 68) 
> 
mew eC33),EFDV(64) ,O@QVEC(33), 
omen gee ey oo 33) CRELO(C33,33),UC68),VC68), 
4GEDIFF(33) 
MPRTEGER P,Q,M,U,V,PP 
Coie. TOC Y.Z,YS,Y1,PFLO,PLOS,YSH,EDIFF, 
oe CRP LO, PFLO1,Z61,ZN1,ZNJ 
BOoUlLVALENCE (CRFELO, Y1) 
eats JAC 
DIMENSION EIDV(64) 


Mee UNLT=30,DEVICE=" DSK’ ,FILE=’ TEST.DAT’ ,DISPOSE=’SA 
WE ) 
Mahou er=-3 5. DEVICK="DSK’,FILE=" DSK.XFR’ ) 


See) VALUES OF IMPEDANCES YSHUNT ADMITTANCES 


C S0658O***xXXTF SHUNT ADMITTANCES ARENOT IGNORED 
REMOVE 
men THE C IN THE LINE YSCQ, P) RRR RAR ARKKEKRKKERAAKKRKRK KE 
CHEK AKKKKKKEKEK KKK KKE RE KKK KEK KE KK EKA RE KKK KRKKK KAR KR KEKE 
READ( 30,2) 1B,El 
Z FORMAT(C1, F) 
Per =CuriEx(E1l,O0.0) 
DO 10 L=1, 1000 
Pees OmogNNG PP. O,JJ,Z,261,2NL,ZNJ 
5 FORMAT( 41, 8F) 
Peers bOo.0)eGO TO 15 
VsChO)=VSH 
C WoO, PHY oH 
Meer et 1 1Cr, 0) V.NE-0)2=1.0/(Y1(P,9Q)4+1/Z) 
alr CO =the 0/2 
eGo re = YC r sO ) 
BeOS Os Y1(P,9),YS(P, Q) 
6 FORMAT(21,4F) 
10 CONTINUE 
KKK KKK KKK KEK IKK KKK KKK KKK KKK KR KA KHER K KK IKE KE KEKE KE KKK KK KK 
BeeouPTLE TERMS FOR YBUS 
CHRKKKAKKK KKK KKK KEK KKK KEKE KK KKK KK EEK KARKKERKKAKKEKKEKEKAEK 
eS DOO P=1, 18 
Y Gore) =0 








93 


WO 40 O=1, 1b 
PiCrebowewoyeGo TO 40 
Meer omer y+ VY 1(P,O)+FYS(P,Q) 
Ge = — CP, 1) ) 
40 Cony ENUE 
50 CONTINUE 
CHAEKRKKAKKKKKKKKKRKEKKKKKKKREKEKKKKKKKKKKKKKKKKKKKEKEKK 
Meer INT OUT VALUES OF YBUS 
CHRKKEKKKKKEKEKKKKKKKKKKREKKKKKKKKEKRKKKKKKEKRKKEKEKKKKEKEKKK 
ERD T <5 | 
51 Borne OF , 10K, BUS IMPEDANCE MATRIX, YBUS‘ ) 
Dee ao, l=] ,18 
ioetieo4.( Y(L,J),J=1,18) 
54 FORMAT( “0° , 1OOF) 
53 CONTINUE 
CAKKKKKAKKRKKKKKERARKKRKKKKREKAKKAKKKEKKAKKEKEKRKKKKKEKKE 
fea VALUES OF POWER AT BUSSES 
CRREKAKRKKKAKKKKKEKKRAKKKEKKKEKKKKEKKKKKKEKEAKEKKKKKEKKEEEKE 
POeo0 P=1,18 
mers), 59) 1,PS,90S 
mel. BO.0)CGO TO 62 
ioe EO.,1 CoO TO 60 
Bscncl)}=—-Ps 
QSCH(1)=-QS 
5Y BORMAT( 1, 2F) 
60 CONTINUE 
CHEK KRKKKKAKKKKKKKKKKRKKERAKKKKKKKKKKKKAKAKRKEKKKKKKEKKE 


BeeeooUhE VALUES OF VOLTAGE AT BUSSES 
CHAK K RAKE KKK ERK KKK KK KK KEKE KR KEK KR KKK KRKERKKE KR AKKEKKKEEK 


62 DO 70 P=2,1B 
me@ny= (1.0.0.0) 
70 CONTINUE 
K=0 


DO 80 I=1,IB 
Tee Seley) 


75 FORMAT(1, 2F) 

80 CONTINUE 

122 PRINT 123,k 

iD 3 men te OM 85%. TPTERATION = ° ,12) 


CRHKEKKKKEKERKKKRKKKEKKEKRKEKRKKEKKKKKEKKKKKRKKKKKAKEKKKEKE 
C CALCULATES REAL REACTIVE BUSS POWERS 
CRAKKARKKRRKKKEKKKKKKKKRKKKKERKARKRKKKEKKKKARAAKEKAKKKEEKEH 

eso) P= 2018 

PSUM=0 

QSUN=0 

me) WAS Ce ie 


P1=REALCE(P))*(REALCE(Q)) *REALCY(P, Q) )-AIMAG(E(Q))*AI 
MAGC(Y 


1(P,Q)))+AIMAG(ECP))*CAIMAGCE(Q)) *REALCY( P,Q) )+RCALCE 





oo 


co) )* 
PAMAGCY'CP.0) ) ) 


mrtAG ( Y 


meee) ))-REALCECP))*( AILMAC(CE(Q) )*REALCY( P,Q) )+REALCE( 
ap) * Al 
BIIAGCY CP, Q))) 
PSUM=PSUM+P1 
QSUM=QSUM+Q1 
140 CONTINUE 
PCAL(P)=PSUM 
QCAL(P)=QSUM 
Peer b45,PCAL(P),OCALC(CP) 
a5 FORMAT(2F) 
150 CONTINUE 
CHKEKKKKKKKKAEKKKKEKKKKKRKKEKAEKKKKEKEKKRAKRKKRKKEKKKKEEKRKKKKKKE 
C PRINT OUT VALUES OF CALCULATED BUSS POWERS 
CRAEKARRKAKKKKAKRKKKKAKRKKAKAKKARKRKKEKRK KKK KAKRKKEKAKERKEKKAR 
BRENT 155,K 
eo 5 momar ( O 4 REAL AND REACTIVE BUSS POWERS K=’ ,I2) 
WOmlss L=2,18 
peeiiel 57 ,PCALCL),OCALCL). 


157 FORMAT (’0’,2F) 
158 CONTINUE 
M=0 
DO 160 P=2,IB 
EP=0.01 


PDIFF(C P)=PSCH(P)-PCAL(P) 
QDIFF(P)=QSCH(P)-QCAL(P) 
TDIFF=ABS(PDIFF(P))+ABS(QDIFFCP)) 
PACGIhr.lt.«EP) GO TO 160 
M=] 
160 CONTINUE 
CRAREKKKKEKKKRKRKEKKRKEKERKEKKKAKKEREKEKEKKKEAEKRKKKERKRKKKKK 
MeeereNT OUL VALUES OF DIFFERENTIAL BUSS POWERS 
CHERKEAKRKKKKKAKRKKAKKKKKEKKEKEKRKEEKRKEKKERKKERKKRKKKERKEKK 
Priv 1655 
165 FORMAT(’°0O’,”’ DIFFERENTIAL PQ = oe yl 2) 


POWbo oc) b=2,.15 
ria werlors. PDEEFCI) SODIFFCL) 
eG 7 Orem i O° 4 2F) 
ros GONTINUE 
GoGlimnGwcjmGoslO 400 
DO aO) ©r =2., 1B 
Nomen LX CPCAL(I),QOCALCL)) 
170 CONTINUE 
CRAKE KAKI KRAKKERKRKERAKAKKKAKKKKEKKAKEKKEAK RK KKK K 
Mec ALOCULATES BUSS CURRENTS ECP) EXCEPT AT SLACK BUS 






100 


CRAKKKKKKKKKKKKKKKAKKKKRRAKRKKKKKKKKKKEKKEKKEKEKE 

Mem > l=2518 

ince) CONE (SCL ))/CONIGCECL)) 

Boe 1 32 1C( 1) 
173 FORMAT(F) 
ee 5 CONTINUE 
CxXeKKKKAKKKKKKKKKRKAKRKEKKKKRKKRKREKKKKKKKKKKKKEKK 
MeereciMLATE JACOBLAN USING PARTITIONS ’°J1,J2,33,J4° 
CREKKKKEKKEKKKKAAKRKEKKKEKEKKKKKKKKEKERKKKEKKEKKEEEKS 
MEACOBLAN ’J1° 
CXKKKKKKKKKKKEKKK 

Pepg=2*1p—2 

TBJ=1B-] 

poel9O M=1,18J 

mel so N=1,1BJ 

P=M+] 

Q=N+1 


JAC(M,N)=(REAL(E(P))*REALCY(P,Q)))+(AIMAG(E(P))*AIMAG 
mar, Q))) 
MEGPRNESO) GO. TO 177 
JAC(M,N)=JAC(M,N)+REAL(IC(P)) 


77 CONTINUE 

c PRINT 178,JAC(M,N),E(P),Y(P,Q) 
1 8 FORMAT(5F) 

180 CONTINUE 

190 CONTINUE 


RKEKKEKKR KKK KEKE KKKKKKKKRKREKRRRKEKKRKKRKRKKEKREKKKRK KKB KK KKEKEKKKREKAKEEKSE 
BEJACOBIAN ’J2’ 
wee KKK RK KKK KKRKKEKRKKKKKKKRKEKKKKKRKKKEKKRKRKKRKRKRKEEKK 
DO 210 M=1,1BJ 
DO 200 N=IB,IBDJ 
P=M+] 
Q=N-IB+2 


JAC(M,N)=-REAL(E(P)) *ALMAG(Y(P,Q) )+AIMAG(E(P))*REALCY 
e250) ) 

iMiGraNHeo) CO TO 195 

JAC(M,N)=JAC(M, N)+AIMAG(IC(P)) 


B95 CONTINUE 
200 CONTINUE 
Z10 CONTINUE 


Mee JACOBLAN ° J3’ 
RHEREKEKKKKKK KKK RHR KKK KKK KKRKEKREKEKKRKKRKE 
Dome t—=1 5. 1 BDJ 
Gwe =), 1 8) 
P=M-13+2 
Q=N+1 


JAC(N,N)=-REALCLECP))*ALMAGCYCP,Q) )FAINAGC ECP) )*REALCY 





LOL 


oe, )) 
ieCPetEe@) CO TO 215 
JAC(M,N)=JAC(M,N)-AIMAGCIC(P) ) 
215 CONTINUE 
220 CONTINUE 
230 Cour tue 


thx keke ee KKK KKK KKK KKK KKK KKRKKEKKKRKEKREEKREKKEKKRKEKS 
SEPJACOBIAN ’J4’ 
Kk wa KKAK KKK KKK KKK RK KEKE KKRKEKKRR RK KKK KEKE KKKEKEEES 

DO 250 M=IB,IBDJ 

DO 240 N=1IB,IBDJ 

P=M=-1B+2 

Q=N-1IB+2 


JAC(M,N)=-AIMAG(E(P) )*AIMAG(Y( P,Q) )-REALCE(P) )*REALCY 
oe Q)) 
MEOPLNE.QO) GO TO 235 
JAC(M,N)=JAC(M,N)+REALCIC(P)) 


235 SON LINUE 
240 CONTINUE 
250 COur TINUE 
bO 334 M=1,1BDJ 
C Pimp 343, 0 JAC(M,N),N=1,18DJ) 
B33. FORMATC1S8F) 
334 CONTINUE 


RK KKERKKKKEK KKK KKK KEK KKK KK KER KKK KKK KKK KKKKKER KKK KKK KK 
KHKKAKKAKKKKREKKEK KEE KK KKK KEK K KK KKK KKK KKK RK KR KK KKK KKKKEKKKS 
PeeecOUPUTE INVERSE OF JACOBIAN 


CALL MINV(JAC,68,IBDJ) 
Skee KKK KKK KKK CREEK RKC KKK KER KKKKEKARKEAKKKKS 


C PRINT 344,(JAC(M,N),N=1,IBDJ) 
344 FORMAT(13F) 
345 CONTINUE 


RHEEKARKEKKKEKRKEKKKKKKKRKEKRKARKRAKKKEKKRKEKRKRKEKKERKEKREKKRKEKKREEKKEEKSE 
Memeo rure DIFERENCE IN VOLTAGE USING JACOBIAN INVERSE ° JAC’ 
KRAURKKKREKEKKKKEKKKAKKEKEKKKKEKEKAKRKEKRAKEKREKA KEKE KRKEKEEREKKREKERK 

pO 346 1=2,18 

ae 

Bowe CCl =PpLrr( 1) 

CON ea COM) =ODIFECI ) 
346 CONTINUE 

Jb =) 

WOM TMS Ted 

EF DVD=0 

EFDQD=0 

EILDD=0 

ELQV=0 

DOB 49 i=l... J 

EFDVD=EFDVD+JACC(CM,N)*PQVEC(N) 
349 CONTINUE 

Doms Oo r=—a EB IBDS 





Giz 


NQ=N-J 
EFDQD=EFDQD+JAC(M,N)*QQVEC(NQ) 
350 CONTINUE 
EFDV(M)=EFDQD+EFDVD 
35 | CONTINUE 
DO 360 M=IB,IBDJ 
ELDD=0 
ELQV=0 
DO 358 N=1,J 
EILDD=EIDD+JAC(M,N)*PQVEC(N) 
B58 CONTINUE 
DO 359 N=I1B,IBDJ 
NQ=N-J 
EILQV=EIQV+JAC(M, N)*QQVEC(NQ) 
359 CONTINUE 
ELDV(M)=EIQV+EIDD 
360 CONTINUE 
DENRA A KAKA KKK AKKKKKKARKKKEKKEKKKKEKRKEKKARKEKREKKKEKEKKRKEKKKRERKKE 
C COMPUTE NEW BUS VOLTAGES USING ’EFDV’ EFDIFFERENTIAL 
VECTOR 
REKKKKEKAEKKEKEKKEKKKRKKKKKKKKKRK KKK KKK KK KEKE KRKEKRKKEKEKKEKKKRKK KK 
Domo 0) P=2,18 
M=p-1 
NeH+J 
LBZ=2*(1B-1) 
MeCNeGr.I1bBZ) GO TO 370 
EDIFF(P)=CMPLX(EFDV(M),EIDV(N)) 
PRINT 368,P,EDIFF(P) 


368 FORMAT(I1, 2F) 
E(P)=E(P)+EDIFF(P) 
370 CONTINUE 
PGK eHOslO)CO TO 390 
K=K+] 
Como. 122 
390 TYPE 395 
395 FORMAT(’ AN ERROR IS MADE IN DATA, LOAD FLOW DID NOT 
CONVERGE’ ) 
400 COM NUE 
500 CLOSE(UNIT=30) 
MEmicn ea K aR RK KK KKK KREKKKAKRKKKAEKKKKEKEKKEKKRKKEKKRKKKKEKKKKKKEKE 
€ COMPUTE POWER FLOWS BUS P TO Q 
C COMPUTE CURRENT FLOWS AND POWER LOSSES 
C COMPUTE SLACK BUS POWER 


KRREEKKKKRKKKRKKEKKKKKKRKRKRR KKK KKK KK KKEKEKEKEKKS 
OPEN(UNIT=31, FILFE=’DSK’ ) 
S(1)=0 
503 FORMAT(L) 
DO 520 P=1,IB 
DO 515 Q=1,1B 
MEGPHEO.0) GO TO 515 
Paeccreonsnn.(0.0,9.0)) GO TO 515 








ro 


PP=P 
Perera KAKA KKK KKAKKKKKKEAKAKKKKHRKKEKKKEARKKEKKKRKKKK 
C IF YOU ARE NOT NEGLECT SHUNT ADMITT REMOVE THE NEXT 


c 
Br HK HAM K UB AKKKKRAKAKRKEKAKKKKREEKEKKRKARKKAKKKEKAKAEKKEKRRKREKRAEREKSR 
momo co NGO CP) (CEC P)—£(0))*Y1CP,Q)) 


C 1+(CONJIG(E(P))*E(P)*YS(P,Q)) 
Pinot Cone ECG))<*((ECO)—E(P))*Y1(P,Q)) 
C Precio iGCcnh(O))=ECQ)*YS(P,Q) ) 


Chino P.O )=PFLO/ CONIGCECP)) 

PLOS=PFLO+PFLO1 
Mie a CK A KK KKK KKRKKEKKEKKKHKHEKEKKKKEKKERKEKRKEKKEKRKRKKRKKRKEEEKE 
C TAKE THE CONJUGATE OF THE POWER FLOWS FOR FORMAT 
P-JQ 
Pcs eae KK RK KKK KKKKEKK KKK KKK KKEKKKKEEEKEKREKEKRKRKE KKK 

Bees CONIC (PLOS ) 

FEo-CONJG(PFLO) 

Dee@eeNE.l) GO TO 505 

S(1)=PFLO+S(1) 


505 WRITE(31,509)P,Q, PFLO, CRFLO(P,Q),PLOS 
WRITE(31,510)PP,E(PP) 
509 FORMAT( 6X,’ POWER FLOW 


eet? , TO’ ,12,°1S’,2F,2X,’ CURRENT 
ManGerS em POWER LOSS 1S ,2F) 


510 Memint( So. BUS ,12,° VOLTAGE 1S’ ,2F) 
5 CONTINUE 
520 CONTINUE 
WRG sl .540) S(1) 
540 FORMAT(6%,’ POWER AT SLACK BUS 1 IS’,2F) 


DO 550 M=1,I1B 
DO 545 W=1,1B 
MeeeARS(CCRFLO(MN,N)).GT.0.000001) 
WRITE(33,543)M,N, 
1CRFLO(M,N) 


543 DORMATC ZT, 2F) 
545 CONTINUE 
550 CONTINUE 
WRITE( 33,555) 
55 5 hOeer ae © ©.0 0.0’) 


Chose (UNIT T=3 3) 

DO 600 M=1,1B 

DO 580 N=1,1B 

See -E( uM) CONIG(CRFLO(M,N) ) 
PFAD=CABS(S(2)) 
Gums) GO. TO 580 
PFAC=-(REAL(S(2)))/(PFAD) 
WRITE(31,570)M,R,PFAC 


57 0 omni ~euRoiieR FACTOR BUS’ ,13,°T0’,13,°1S’,F) 
580 COME IUE 
G00 CONTINUE 


CEOSECUNIT=31) 





oO 


mo 
sy 


29 


30 


5 


Bo 


40 


45 


50 
35 


60 
62 
5 


70 


n5 


80 


100 


SOP 
END 
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Pen steve OUD ROUTINE FO INVEKT A REAL SQUARE MATRIX 


SUBROUTINE MINVCA,MDIM,N) 
DIMENSION ACMDIM,N) 
INTEGER L(100),M(100) 


CONVERimhD FROM SSP ROUTNE MINV 


DiowaO k=), N 
mer) k 

MC K)=K 

pee A=ACK. K) 

DO) 20 J=K,N 

Dom2O L=K,N 

mers BIGA)—-ARS(CACL,J))) 15,290,290 
Been ( 1 Ii ) 

LC K)=1 

Meloy= J 

CONTINUE 

Gk ) 

ieG@i{(C) 354135; 25 

Boe sO l=, 
HMoiED=—ACK., 1) 

mG )=ACi, 1) 

Gg) OND 

I=M(K) 

WEG Ko 5, 38 
DOO = 1. N 
HOLD=-A(J,K) 
ACJI,K)=ACJ,1) 

howe) HOLD 

DOs eel— le 

itt) 50.55, 50 

Nie y= AC 1, K)/(—BIGA) 
CONTINUE 

DOmoon l= yo 

FOUL, BAU C I 1): 

Ome omer 

Ter K ) 900165 , 60 
DiGi) 62555, 62 
Noles — Omp~ ACK, J)+AC1,J) 
CONTINUE 

DOmmeo = lor 

Peer) 70,/75,70 

Gio = AC Ke), BLGA 
CONTINUE 

Ne wne)=( 1) .0090000)/ BIGA 
CONMINUE 

eh 

K=K—-] 





BOD 


iota. lo), 105 


oO5 =k) 
he Glsierel2@. 120.108 
POS Doerr J=1.N 


HOLD=A(J,K) 
ACJ,K)=-ACJ,1) 


110 AC J,1)=HOLD 
120 J=M(K) 

mej) 100,100,125 
5 DO 130 I=1,N 


HOLD=A(K, 1) 
AC T)==A( J, 1) 


Po 0) del = 10LD 

GO TO 100 
150 RETURN 

END 
KKK KKKKRK KKK RAK KEKE KKK KKK KEK KKK AR KKARKKKARARKKKKARKEKER 
C Memes Woeta ACCHRRTS DATA FROM TEST.DAT AND DSK.XFR 
FROM 
C ROU OWePROGKAM LDFLO.F4 AND FORMS ZBUS BY 
C Siemens tiles NILE BE USED TO COMPUTE ALL 
PAU LT 


C CUM to. OAD Fh GOW COMPUTED MAXIMUM LOAD 


CURRENTS. 
KKEKKKKEKEKEKKKKKKKKAKHEAEKRERKE RK KERR KKKAEKEKRKKEKKKEEKEKKEKKAKES 


RKEKEAKAREKKKKKEKKKEKRKEKEKKKEKKREREKERKEKEAREKKEKKRKAKKKKEKKKKKERAKEKEEE 


C Reo DATA FROM TEST.DAT 
KHKKKKKAKKKK EKA KKK KEE RARER KKK ARK KEK EK KER KK KR KKKKAKRKEKAKEKR KE 
DIMENSION 


mee), 3), 210300, 2),RID(30, 30, 3),ZG(30),ZN(30), 


mime O00 ),R1(30),Z28(30,30),2B0(30,30),2MU(5,5),1DD(5), 
meC30), 
PUaBItG@sO), Z60(30),2L0030) ,EB(30) 
COMIC ZN.2C,TI.ZN,ZN1,ZN1,ZNJ,Z2G1,Z0, 
RaReaABOMZi 2hB,2L,20b,ZBR,ZC0,ZL0,ZLLO,EB 
NNB=] 
EQUIVALENCE(ZBR, ZLO, EB) 
OPEN(UNIT=31,FILE=’TEST.DAT’ ) 


6 OPEN(UNIT=32, FILE=’DSK.-XFR’) 
ReApcaiy lO) 1B. EB! 
10 FORMAT(I1, F) 
20 Pees iesOo Mb. Lal, JJ,T1,261,Z2Ni, ZNI 
30 FORMAT(41, 8F) 
PiGweEO.d) GO TO 40 
DN=NB 
NN=NB 


Cet = CMPLXC DN; 0.0) 
ZDCNN,2)=T1 

ZW, 3)=Z26C14+3*ZNC1) 
ZWCI)=ZNI 
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ZNCJ)=HZNI 
RICCI yy EOE 
PaCvbeNeaO)GO TO 26 


40 DO 45 I=1,IB 
C READ(32,43)T1 
43 FORMAT( 2F) 
45 CONTINUE 
DO 48 I=1,IB 
C Rea 46 ol, Tl 
46 FORMAT(1, 2F) 
EB(J)=T1l 
48 CONTINUE 
49 RemNbes1,50)1,T), XDI 
50 FORMAT(1, 3F) 


ier. 4.0) GO TO 60 
MeCGABS(T1).EQ.0) GO TO 49 

Cc Mi GCONIGET 1) 

fc ZGOl )y=3*(CABSCEBC(I))**2)/T1 
Mmme@oD lane .0)2GCL)=C4UPLXCO.0,XD1) 
BeEocl )=ZECL ) 
Mui. Nis) ) GO TO 49 


RaARKEKKKKEKEKRKKKRKKKKKKREKKEARKREKERKEKKREAEKKREKRARAKEKRRKEKKRKKK 


C ENTER THE ZERO SEQ IMP BUS TO REF(GROUND) 

Matec Aki k eee KRKAEKKARKRKKERKKKHRRERKEKKKREKRKK KAR K 

60 TYPE 62 

62 FORMAT(’ ENTER THE SYNC ZERO SEQUENCE IMPEDENCE FOR 
THE LOADS 


Pee Aon oe BUS. 3°2N,4MUST BE ADDED, IMPED IN PU..,’) 
DO 65 I=1,I1B 

MGGEPT G4, 11,T1 

ail ~b0.0)GO TO 66 


Ameo iy )=T) 
64 FORMAT(1, 2F) 
65 CONTINUE 
66 Dome srl, is 
PieGAnSe7C0 (1) ).hOQ.0)2G0(1)=ZG(1) 
68 CONTINUE 
PERERA a uw AK aKKKKKKEAKAKKAKKEKKEKRKKKKAEKEKAKRKAEARAREKRKEKKRERKAKRKNSE 
ce ENTER THE MUTUALLY COUPLED ELEMENTS 


RARAEKKKEKKKKEKKKAKEKEKEKEKKEKKREKEKKEKKKKEKRAKEKEKEKKRKEREKAEEE 


CEOSECUMLT=31 ) 


iets 70 
70 Oe went eRe nee MUTUALLY COUPLED ELEMENTS IN PER 
rar’ /’ 
womb ovw ie rORMAT. TEST. DAT. HAS THE ELEMENT 
a, /- 
Peewee l oe TOCEEREMENET )SP(O SEQ VALUE) ‘%) 
75 AGGERT. 30 .,EL,ELL,Z0 
&0 FORMATC4F) 


MGue aoe) CO. TO 100 
ete 
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ZMC(II,1)=CMPLXC(EL, ELL) 

Bact), 2)=Z0 

eoro 75 
KKeRKKK KR KK KKK KKKKEKEKRKKEKKEKKKEKKKKKKEKKRKKEKEKKARKKRKEKEE ESE 
C ALL ELEMENT DATA PREPARED PERFORM POSITIVE AND 
NEGATIVE 


C SEQUENCE COMPUTATIONS TO FORM ZBUS 
KKKKKKRKEKKEKKKEK EK KKK KKK KK ERR KR IK KKK ERK KKRKKKKKK KE RKKKKKAK KK 


RRR KKEKRKKKKKRKKAECKKRKKKERKRAEKKREKKAEKRKKKRKRRKKKKKKRKEKKKKKEKRRKKKK 
100 DO 500 IEL=1,NN 
BRA=0 
ID(IEL)#=1 
nO 120 .L= 1. 30 
DO 120 J=1,30 
Ome OOeei=1, 3 
Peer Ch iJ J).NEeTEL) CO. TO 


120 
Behr) sbO. 0) GO TO 
LO 
BGR) .ee 0) Com rO 
aS 
BRA=].0 
PeP=MiNOC 1, J) 
IQQ=MAXO(1I,J) 
COMEO 130 
110 PEGhelned a Oa) ) GO et © 
fel 5 
RI(1I)=1.0 
IPP=J 
I1QQ=I 
Goyer O a 30 
lS Tee) ) = lene) 
RIC1I)=1.0 
PPP=4ENOC LSJ) 
ITQQ=NAXOCI, J) 
GOwtO 30 
120 CONTINUE 
KKEKKKEKKKEKKKRKKAKKKEKKKERRAKREKEEKAKREKRKEKKKKEKEKKEARKRAKEKKEKKE HK 
C Peed Oke NK WAS BEEN DETERMINED BRA=0: BRANCH 
KRHEKEKKKEKKKKKEKKEKKKKKKEKKEKKKEKRKKEKKEKKKKEKREKKKEKKKKEKRKKKEKKKKKE 
130 DOrwo> IL=),18 
Ze CMPIEX(0.0,0. 0) 
181 PPCBRA.WE.2O) GO TO 186 
Peano. nO.lt) 'GO LO 190 
Znivell) =2b( TPP, 1L) 
Coyi1o 190 . 
186 een) =Z2rtebere tie LEC LH), 1L) 
190 CONTINUE 
95 CONTINUE 


Zee MP Ly C060, 0.0) 
IFC BRA-NE.O) GO TO 250 
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HOr2oo Th=1. 1B 
AGUeO ml) =2 BR CLL) 
AGEL TOO) =Z8R( TL) 


200 CONTINUE 

240 Zero =AnCiP Pe LOO) ZDCLEL, 2) 
Dopo EX CO. 1, ().-0') 
COP lO. 3.00 

250 ZA Obra beroaO}arZoCIEL, 2) 


DBR =seMPEXCOeOn 0. 0) 
DO 290 IP=1,1B 
DO 290 IQ=1,1B 
PGureer ae) coO) 10.290 


mmenP,1Q)=ZBC(IP,1Q)-(ZL(IP)*Z 


mc Q) 
i 11.) 
Pee n Owe) SCOrTO 290 
Zh ir )=Z2BC LP. 10 ) 
290 CONTINUE 
SOO D Omer Ome" 5 
DOms50 J=1,5 
Aen) = CiuP EX (0.0.0.0) 
meal = 
350 CONTINUE 
Di eO Oe eee 
Donec C10 sn — alas 
Pies sou..!. J,ZBC1L, J) 
380 FORMAT( 21, 2F) 
400 CONTINUE 
NNB=1 
500 CONTINUE 
Peirce sce oo oe Oe AK KK RAR KKK KKK KKK KKK KRKRKKKEKREKEKE HK 
GC ZBUS POS/NEG IS FORMED 


AReaeaRKAKKAEKAKKAARKKKEKRARKEKKEEKRKKKKKEKRKKKEREKKRKEKKKKEKAKRKARERKARER 


DO 505 1=1,NN 


IDC 1)=0 
505 CONTINUE 
perslo T=). 18 
RI(1)=0.0 
10 CONTINUE 
DO 1000 IEL=1,NN 
BRA=0 
ID(IEL)=1 
NNB=1 


bDO 620 1=1,30 
DO 620 J=1,30 
DO 620 JJ=1,3 
mre wen@iues. JJ).NE.LEL) GO. TO 


Pearl Leb .0).CO. TO 
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ECR TG Is bs O 160-10 


ot 5 
BRA=1]1.0 
Pera TNO CL. i) 
POC=MANO CI J) 
CO) 10 625 
610 PGR IC J) .£0.0)) CO. TG 
65 
Riel =). 6 
IPP=J 
1QQ=I1 
GOMLO7 i625 
or 5 RI(1)=1.0 
RIC J)=1.0 
IPP=I 
HOO = 4 
GOVT Oro 25 
620 CONTINUE 
MPM nA KAKA KKKKEKRKARRKEEAEKKKEKKKKEARKREAKRKEKRKEKKKKREKEKRER KS 
C Seach FOR MUTUAL COUPLING IN ZERO SEQUENCE 
KRKEKKKKKEKKKKARKARKREKRKKKKKKAEKRAKRERKARAKRKAKKARKKKAKKKAAKAKRER KK 
625 IDD(1)=1EL 
DO 650 I=1,NN 
Pie 
TE CREAIRCZMC1, 1) )2.EO.D1)MUT=) 
LPCATMAG(ZMC1,1)).590.DMN)NUT=1 
PeCuviene.1)GO 10.650 
Mu —O 
Pile th LX REAL ZMG1, 1)) >) 
fee Wei inGie ZMC1,1))) 
triewghert |) .20.0)Co) TO 650 
Bvem@reimi2).60.0)GO TO 650 
PGmmierrOsaleL)GO TO 630 
mis = Pr] 
I1l=112 
L12=113 
630 EOC y=1 11 


NNB=NNB+1 

Tee (NNB)=112 
Dil=E12Z 

DOTGA > J l=) ,NN 


ieee EALCZH(CJ1,1))-EQ-DM)MUT=1 


MmereATMAGCZMN(J1,1))-<EQ.DM)MUT= 

] 
ete. |) GO TO 645 
MUT=0 
meee ek ALC ZC J1,1))) 
ieee I PIAG{ ZMNCJ],1))) 
Paoli. sy. O0)GO TO 645 
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Pro 0ll2),60.00G0 TO. 64 5 
PaChhliskoe DM) GO ro 643 
ers —ir 
eb =i 12 
it 227 1) 
643 een roetli) CoO TO 645 
NNB=NNB+4+1 
IDD(NNB)=I1112 
645 CONTINUE 
650 CONTINUE 
DO 660 I=1,NNB 
ZMU(1,1)=ZD(IDD(I),3) 
DO 660 J=1,NNB 
DD1=IDD(1) 
DD2=IDD(J) 
DO 660 KK=1,NN 


LF(ZMCKK,1).EQ.CMPLX( 
Die DD2) ) 
IMUT=1 


meee RK, 1)-£EQ.CMPLX(DD2, DD1) 


) MUT=1 
IF(MUT.NE-1)GO TO 660 
MUT=0 
ZC Te J YaZ OKK , 2) 
ZMU(J,1)=ZNCKK, 2) 
660 CONTINUE 
Pere FO HK AKAK KK KK KKK KRKKKAAKKRREKKRKEKRKKKAKHKKAKKEK 
C ZMU IS FORMED 


Perec KARR KKK KKK KKK KKK KKK KARR 
DO 665 L=1,NNB 
DO 665 J=1,NNB 
PRINT 663,1,J,ZMU(I, J) 
663 FORMAT(21, 2F) 
665 CONTINUE 
Meneses te Hr Ks a AK KARE KK KKK KKK KKK KKRKEEKKEEKEKEH 
CALL CMXINV(ZMU,5,NNB) 
DEM RK KK RRR ARKH KKKKEKRKKEKKKRKKKEKKKKKKKKEKEKRKEEKKSE 
DO 668 I=1,NNB 
PR o67,(Z2NUC(1,J);J=1,NNB) 
667 FORMAT(10F) 
668 CONTINUE 
DO 695 IL=1,I1B 
ROGUNE EQ.) CO TO 621 
AS B= GMPL NM  Oa0), 0.0) 
DO 680 J4=1,NNB 
J=J4+1 
I=IDD(J) 
PGE orOon GO TO 630 
DO 670 Tl=1, 30 





8 a 


DOUG 0s eb 2 tees © 
DOM] 0 


LF(RI 
meni ,12,13 
Deiat ycomuo 670 


IP=11 
ee 1 2 
Ae «iy Mite) *(Z2ZBOCIP,LL)-ZBOCLOQ,IL))+Z BB 
670 CONTINUE 
630 CONC LNUE 
bo | PGB AS WE. 0) GO TO 686 
EEO sro. lL) =GO TO 690 
ZERCTLY=ZBOCILPP,LL)+ZBE 
GOO. 690 
686 Ani j=- Soba CmePr TL )—ZBOCTOO, IL) 
690 CONTINUE 
o> CONTINUE 


ZBB=CMPLX(0.0,0.0) 
Me CBNeNE. 0) CO TO 750 
DO 700 IL=1,IB 
ZBO(I1QQ,IL)=ZBR(IL) 
Zr Lao l=Z8R CL) 
700 CONTINUE 
MeGinibebowl) CO TO 740 
DO 730 J=2,NNB 
I=IDD(J) 
DO8720 Ii=1.30 
pO 720 I2=1,30 
DO 720 13=1,3 
GaermpGi le lo «19).NE.I) CO TO 720 


IP=I1 
IQ=12 
20 CONTINUE 
me—=(ZNUC1,J)*(ZBOCIP,1Q0)-ZBOCIQ,1I QQ 
y/ 
LAC 11) + ZEB 
730 CONTINUE 
740 Mere Oo )=7B0C1PPs1OO)+ZBB+1.0/ZNU(1,1) 
fee enh in(O.0,0.0) 
Gor hOmu0 G 
70 PeCNNBoEQ.1) GO TO /81 
Dery s0 J=2,NNB 
I=IDD(J) 


DOw/7 7 O11 —1., 30 
DO 7 OF el2=1 4.30 
DiGmn 0 eL3—)] ,3 





ie? 


CR OGL Ie, 
malo ele. 1)GO TO 770 


eal 
1Q=12 
770 CONTINUE 
mee ZHUCL,J)*(ZLC(1LP)-ZLC(1Q))/ZMUC1,1 
) )+ZBB 
780 CONTINUE 
781 Aub weber yale LOO)+2BB+C1.0/ZMUC1,1)) 


Zoi = Cre LX 00.0, 0.0 ) 
DO 790 IP=1,1B 
DO°7 90 LO=l, 18 
mele ecil. 1O)CO 10 790 


ZBOC(IP,1Q)=ZBO(IP,1Q)-(ZLC(IP) 


ar 1Q) 
ey ZL) 
(PE GGeeaOn LO) GO TO 790 
Aomori y=2 60 CLP, 10) 
790 CONTINUE 
800 DOr es 50s l=) 5 
Horton = 1,5 
Ze) = OME Xx (O.0,0.0) 
Loc lO 
850 CONTINUE 
Borovgon t=1.,15 
Vero s=1. 15 
Bites oO. a. 2b0C1, J) 
330 FORMAT(21,2F) 
900 CONTINUE 
NNB=1 
1000 CONTINUE 


Dee AAA KRACK KKK KKHKEKKKKKERKKHKKEEKKKEKKRAKKKEAKRKEKER 
nok a oe Kx 
C ZBUS ZERO IS FORMED 
Dts oO KORO KK RAR KKK KKK KAR KKK RK KKEKRKKEKKRKKKREKKHRAKHK 
OPENCUNILT=33, FILE=’ZBUS. FOR’ ) 
Dolo TO 1.18 
DO 1010 J=1,1B 
Pam Gomes) Ieee ZB0CL. J),ZBC1,.J3),ZB(1,J3) 


1005 FORMAT C21, 6F ) 
Pet 10107 
OO 7 i Ce enstommelOeumbuS  ,3%, ZERO SEQ “ ,5X,° POS 
SQ ) 
Poeun J. ZBOCT,J),ZBC1L, J) 
1008 FORMAT(C216,4F) 
ol © CONTINUE 


CEOSE CUNT T=3 2 ) 
KRARKKKKKKKKKKKEKREKKEKKKKE KEE KKK KEKE RAKKKAKKEKKEKKEKKERAEK 


C vO wtb LOADS. ALL LOADS WILL BE LINKS WITH 





eS 


NO 
€ MUTUAL COUPLING. 
KEE KKKKKEKKKKKKKREKEKREKREKKKEKKKKREKKKKKKKERKRKKRRRKKRKKRKKK 
pemtz00 T=), 18 
iiGe wo C7G GMOs 0) GO TO 1200 
DOM LOO Mutts 
Punta) == 2 BC 1. 1 1.) 
AiO y=—-2 R001, LL) 
1100 CONTINUE 
ZLL=-ZLCI)+ZGC1) 
ZLLO=-ZLO(1)+ZG0(1) 
DO 1150 J=1,I1B 
DO 1150 K=1,IB 
ieee Gor To. 1150 


mee K)=ZBCJ,K)-ZL(J)*ZLCK)/Z 
LL 
ZB = 2B (J, kK) 


me, kK) =Z2B0(J,K)-ZLO(J)*ZLOC 
roy, ZLLO 
ZBO(K,J)=ZB0(J,K) 
ba 50 CONTINUE 
1200 CONTINUE 
Come CUNm@=-354 hILE= ZBUS.DAT’ ) 
DO AO) 1 = ae: 
mor 1300 J=1,18 
Mivmencouel2 50)1,J,Z28(1,J3),ZB0C1, J) 
250 FORMAT(21,4F) 
1300 CONTINUE 
CLOSE(UNIT=34) 
moon CUNT T=3 3) 
GeG@uUnNli=36,FILE=’MNUT.CO’ ) 
Domes oo  1=—1, NN 
Mim PiAnCZMCL, 2 )) 
GCE nESEQ.0.0)GO TO 1322 
Di —AMYAG(ZMC1.,1)) 
IELI=IFIX(EL1) 
PERQ=1FIX(CEL2) 
pemmla70)Lil=1,18 
Hoes 20° Jl=1,18 
Doee1320 Ki=1,18 
Peele Kl).kO.ELL)GO TO 1310 
MEGmnmG@btn sd) ,K1).5fO0.ELZ2)GO TO 1315 
ComLoul 320 
310 ieee | 
Poy) | 
COecro 11320 
315 ii =1 1 
mee 1 
20 CONTINUE 


— 





mo 2 
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POO. )GO TO 1400 


Prema t6., 1350)EL1,EL2,IBIL,1BIJ,IBJL,IBJJ,ZM( 


ilies 2 ) 
1350 
1400 


1450 


C 


eo) 
i 


aG 


25 


30 
5D 
BS 
40 
45 


=) 
.D 


60 
oz 


POMmiIAt(2F 413 .2F ) 
CONTINUE 
Wight pes Gols 5 ) 
HQ Gro. 0.070 0 0 0 0.0 0.0’ ) 
CLOSEOUT =3 6) 
STOP 
END 
TITLE: CMXINV-SUBROUTINE FOR INVERT COMPLEX MATRIX 
SUBROUTINE CMXINV(CA,MDIM,N) 
COMMER wevChD IM. N) > BIGA, HOLD 
eG 10.0). MOC1lLO0O) 


COM. hE DO RKOM SSP ROUTINE MINV 


Dome 0) k=l. N 
LCK)=K 

MiG = % 

Bee A= ACK, kK) 

DO 20 J=K,N 

bO 20 I=K,N 

Meee Sen GA )—CARBSCACL,J))) 15,29, 20 
BlGASAC 1. J) 

L(K)=1 

Gros 

CONTINUE 

SLO) 

ee 35,35 ,25 

Dome oe t- 1 N 
HOLD=-A(K,1) 
mR. Te Cp) 
ele) HO 1 1 

I=M(K) 

DEM!) 745.45, 38 
Demet — lh 
HOLD=-A(J, XK) 
Nein y= AC), 1) 

Gia )—HoL) 

Domo l=) N 

@r—-1 ) 50.55, 50 
ACI,K)=ACL,K)/(-BIGA) 
CONTINUE 

DOmoo l=) ,N 
Hen A Ci, kK) 

Die oe t= 1, 1 

Piem=3o) = 60.65, 60 
Pete) 62.6 5,62 
Glo WOen ACK. J)FACI,J) 





Lie 


65 CONTINUE 
DO 75 J=1,N 
IF(J-k) On? 527. 


70 ACK, ID=ACK,I)/BICA 
75 CONTINUE 
atk» K)=(1-0000000,9.0900000)/BrcA 
80 CONTINUE 
K=N 
100 K=K~-] 
EF(K) 150,150,105 
105 I=L(K) 
ee eo) nl 2 io aie 
108 DO 110 J=1,N 


HOLD=A(J,K) 
ACJ,K)=-A(J,1) 


110 A(J,1)=HOLD 
0 J=M(K) 

at) 008 100. 125 
125 DO 130 I=1,n 


HOLD=A(K,1) 
ACK, 1)=-A(J,1) 


130 ACJ,1)=HOLp 
GO TO 100 
150 RETURN 
END 
GC nap pttttAsstnnansinsunssnacersnaunsrnvanenac.. 
C TITLE: FALT.F4; TAKES DATA FROM ZBUS.DAT AND COMPUTES 
G ALL FAULT CURRENTS AT ALL BUSES 


| OIHERRNAASAHSASUNSSAdanERANAsEED ESE EE, BUSES 
DIMENSION 
Be 120039, 30), PcuR(300,3,2),era17(30, 5. 9) 
~2 BC 30), IDELE( 300) ,E1B( 305 
COMPLEX 
pags Ak ™S®*0E21,201, 2801, 281,201,201, 293, 
pee e* 2). MDELE,FCURI, FCURD alccs 
eee a OP ME “Pest. par’ ) 
ie een Zbus pat) 
READ(34,5) IBLE] 
5 FORMAT(I,F) 
IBB=IB¥IRB 
DO 7 I=1,1B 
pe CUPIX (CEN, 0.9) 


7 CONTINUE 
9 eee tO) NT Tey. Kk 
10 FORMAT(4T) 


FENESEQ.0)GO To 15 
BG ee. 2)) =: 
LE@i@eHOMms)) j=— 4 

RI=I 

RJ=J 
'DELE(NI)=CMPLX(RI, RJ) 





(eS 


IPCRE~NE.9)GO TO 9 


i) Hoe soe i 1B 
Rego 25) 1, J3,2B.2B01 
2G FORMAT(21, 4F) 


ZBI, J)=Z281 
ZBL, J )=Z B01 
30 CONTINUE 
CLOSECUNIT=35) 
PoaeuNbi=37,FTLs= MNUT.CO-) 
CVG i—-so.rliLE= FALE. DAT’ ) 
OPENCUNIT=39, FILE=’EFLT.DAT’ ) 
OPENCUNIT=40, FILE="RELA.DAT’ ) 
OPENCUNIT=41,FILE=°FCUR.DAT’ } 
OPENCUNIT=42, FILE=’ DSK.XFR’ ) 
SPE NCUNDT=43,FLLE= LOURsDAT’ ) 
READ(40,32) RMVAB,RKVBA,IRR 
a7 FORMAT(2F,1) 
RKVBA=RKVBA*1E3 
FLPU=RMVAB*1E6/(CSQRT(3.0)* RKVBA) 





33 READ(40, 34) IRR 
34 FORMAT(L) 
IF(IRR.NE.O)GO TO 33 
36 READ(40,37)J,E 
37 FORMAT(I1,F) 
BCI <fO.0)CO TO 38 
ELB(J)=E 
GO TO 36 
38 TYPE 35 
35 FORMAT(’ ENTER THE FAULT 


IMPEDENCE-REAL(SP)IMAGINARY..,’ ) 
ACCEPT 40,ZF 
40 FORMAT(2F) 
DO 500 IFLT=1,I1B 
MiGenkocZhGlrin sIFlT)).EQ.0)GO TO 500 
DOmaOl T—leei8 
[MGmerO.1FERT)GO TO 45 


mayen (L,2,1)=SORT(3.0)*(EB(1)-((ZBC(1, 
LFLT)*ER 
MGMT) C7 R+2n(1FLT, IFLT)))) 


Meee CL, )],2)=-(SQORT(C3.0)* EBCIFLT)*ZBO 
omer LT) / 
ea OG hel tidy oor 2*2BCLFLT, IF LT) ))) 


feet Cl, 2,2)=SORT(3.0)*EBCI)-CCSQRTC3 
-O0)*EB 


eee) 2D DPT) )/ CZBOCIFLT, IFLT)+(2*(Z8BCLFLT,IFLT) 
De o*ZF))) 





ey 


Meee ( 1,3,2)=EFALT( oe ZyVeCSOorRT (3.0 )* 


Bee 1 ) ) 
Commo 50 
45 
Meee IFLT,2,1)=SORT(3.0)*(ZF*EBCIFLT 
\/CZF+ 


Pao PRET. LEE.) ) ) 


PRALTCIFLT,1!,2)=-SQRT(3.0)*EBCIFLT)*Z 
Boel E LT, 
LionDy earn OCnneme TELT)+(2*Z8(LFLT,IFLT) )+3*ZF) 


Pia@eee( TF LT, 2,2)=SQRT(3.0)*EB(IFLT)*(Z 
mom LELT, 


Mmeemer + ZBCIFLT, 1FLT)+3*ZF)/CZBOCLFLT,IFLT)+(2*ZBCIFLT 
meet) )+3%* 
DEP) 


EFALTC(IFLT, 3, 2)=-SQRT(3-0)*EB(IFLT)*Z 
Perr LT, LFL 
Mey (2BOCIFLG. TP ET)+(2*ZB( 1ELT,IFLT))+3*ZF) 
50 CONTINUE 
REWIND 34 
REWIND 37 


meme 39,53)CIFLT,1,EFALT(1,2,1),EFALT(1, 1,2) 


> 


Mee Geo. 2) BEFALTCT,3,2),1=1,18) 


53 FORMAT(213, 8F) 
BBADC 34.55) L1,5 
BD FORMAT(I1,F) 


DO 90 K=1,1BB 

READ(34,60)NN,1,J,KK,ZR1,ZN1,ZNI,ZNJ 
60 FORMAT(41, 8F) 

ie Cuneo. co Too | 


en Oe CmrALT¢ 12. ])-BFALT(J,2,1))}/ZR1 


meme NN, 1,2)=CEFALT(1,1,2)-EFALT(J,1,2))/(ZN1 
+3%* 
LC ZNI+ZNJ)) 


meee 2 )=CbhrALT( 1,2,2)—-EFALT(J,2,2))/ZR1 
FCURCNN,3,2)=C(EFALT(L,3,2)-EFALT(J,3,2))/ZR1 


DO 85 12=1,1BB 


fe 37,65)RN1,RN2,1B11,1BIJ, 
BBII,IBIJ 
1,ZN2 
65 FORMAT(2F,41, 2F) 
NNI=IFIX(RND) 
NN2=IFIX(RN2) 


— ; 
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GER es ieeiealv CO Oey O 
Peewee cored ehG On ak@ 74 
PECNNIT.EOQ.0)GO TO 36 


€07 10 65 
70 emer) — (CR PAD EC TBJT, 1,2)-EFALTCIBJJ,1,2))/ 
1(ZN2))+FCURCNN,1,2) 
Co tO 65 
74 eviMeuieme2z ) —( (EPALTCIBII,1,2)-EFALTCIBIJ,1,2))/ 
1(ZN2))+FCURCNN, 1,2) 
85 CONTINUE 
86 REWIND 37 
90 CONTINUE 
91 REWIND 34 


DO 300 NN=1,I1BB 
IF(CABS(FCUR(NN,2,1))-EQ-0)GO TO 300 


WRITE(38, 250) LFLT,NN, FCURC(NN, 2,1), FCU 
mid, 1,2) 
PRG URCNN, 2, 2),ECURCNN, 3,2) 
0) FORMAT(213,8F) 
300 CONTINUE 
DO 400 NN=1,18B 
MeCewpoch elicit 2/1)).EO.0)CGO TO 400 


fee T= CRKEVBA/EILBCIFIXCABSC(REALCIDELE( 
NN) 
ino) ) +? GURCNN, 2,1) 





FCUR2=(RKVBA/ELBCIFIX(ABS(REAL(IDELE( 
IN) 
Dy) )*(FCURCNN, 1,2)+FCURCNN, 2,2)+FCURC(NN, 3,2)) 
FCURL=FCURI*FIPU 
FCUR2=FCUR2*FIPU 
IRI=(LFIX(REALCIDELE(NN)))) 
IRJ=(LFIX(AIMAGCIDELFCNN)))) 
CUR1=CABS(FCUR1L) 
CUR2=CABS(FCUR2) 
aCe row oul) CO TO 400 
WMH Gale 350 )IFLT, NN, CURI,CUR2 
350 FORMAT(21, 2F) 
400 CONTINUE 
500 CONTINUE 
DO 600 I=1,I1BB 
erenes? 525 lL, JL, ALCUR 
MZGuEeEOr Oo GO TO G10 
MGC NOUR) .LE.0.0)GO TO 600 
BLCUR=CABS(ALCUR) 
CLCUR=FIPU*BLCUR*RKVBA/EIB(JL) 
BLCUR=FIPU* BLCUR*RKVBA/ELBCLL) 
UMPCs 3525 Led, RLOUR, CLCUR 
B25 FORMAT( 21, 2F) 





600 
610 


602 


CONTINUE 
Wirt CAle O02 ) 
WRITE(43,602) 
DORtMAre 90 0 0 0 0 
CLOSE(UNIT=36) 
CLOSE(UNIT=37) 
Cost CU Nir = 38) 
CLOSE(UNIT=39) 
CLOSE(UNIT=42) 
CLOSE(UNIT=43) 
STOP 

END 


7 


jee, 


— 





20) 


KREERKKAKKKKKKEKKKKKERKAKEKRAKKRKKEKRKEKKKEKKKKEKKEKEKKEKEKK 


C Deron asUoloe DATA FROM INPUT.~.F4,FALT.~F4 TO 
COORDINATE’ 

6 ieee oneel LON DEVICES FOR THE WHOLE DISTRIBUTION 
exo lkEM. 


KAR K KKK KKK KK KKK KRKEKAHKRREEKKRKKRKEKEKREaAERRAEKKSE 
OPEN( UNIT=30, FILE=’ RELA.DAT’ ) 
OPEN(UNIT=31, FILE=’LOCUR.DAT’ ) 
OPEN( UNIT=32,FILE=’LCUR.DAT’ ) 
OPEN(CUNIT=33, FILE=’FCUR.DAT’ ) 
OPEN(UNIT=34, FILE=’TEST.DAT’ ) 
OPEN(UNIT=35, FILE=’XFOR.DAT’) 
OPEN( UNIT=36, FILE=’DEV.DAT’ ) 
OPEN(UNIT=37, FILE=’ ADEV.DAT’) 
DIMENSION EB(30),CUREL(2900),CURLO(30),Z(100) 
Pe enema ALDEN (100, 7),xXFR(200,2),11D(3),R1D(100 


PDE y C2000) 4) 
Cour E 7 SAnDEY sel )c2.71,22.R1D,CUREL, CU1,CURLO,XFR 
Kee KK KRACK KKK KK RAKE KKK KRRAKRREKRKRKRKKKKRKAKRKKKaEKH 
G READ IN LOAD AND BUS VOLTAGE DATA 
wesw KKK eR RK RRR RRR KKHRKRKEKKRKEKKKKKKEKEKKAAKRERKRKKS 
READ(30,5)R 


5 FORMAT(F) 
READ(34,10)1B,El 
10 FORMAT(I, F) 
hi READ(30,12)1 
2 FORMAT(1) 
eels Ngeaonco TO 11 
13 Rape sO, 4), El 
14 FORMAT(I, F) 
MGeehO.0 eo TO 15 
EB(1)=El 
Como 13 
15 Dem20 l=), 1B 
REGS 1. 1G iJ, CU2 
16 FORMAT(I, F) 


PGiskOs0)GO TO 921 
CienoC I —CMPLXCCu2*1.25,0.0) 
20 CONTINUE 
21 pee24) 1120 
Mew Ome DEV( 1, J),J=0,4) 
Pec hO.0 Go TO) 23 


D2 FORMAT(1,5F) 

24 CONTINUE 

KEE KEKKEKKKKEREKRKEKEKAKKEKEE KEKE KEAKAEKEKKAKEEREKEK SE 
C NOW TO SET ALL LOAD DEVICES 


KAKA KE AK KKK KKK KKK KA KKK KKK KAKA K KK AKA KARKKESR 
23 NaraEe «2:5 








Sa 


deze 


ZS FORMAT(’ ALL LOAD CURRENTS ARE READY,USE DEVICE 
NUMBER LIST,’) 

DO 100 I=1,1B 

MCOnbSCGthbO( TyyebO.0)GO TO 100 


ID=1 
30 TYPE 35,I1D 
35 FORMAT(’ SELECT THE DEVICE FOR THE LOAD AT 


Pee, 13,’.-.,°S) 
ACCEPT 40,IDEV 
Prev eviEC lL, L)=FLOATCIDEV) 
40 FORMAT(1) 
CUbl—-k EAC CURLOCE)) 
CALL 
een CURL, EB(1),AIDEVL(I,1),AIDEVL(1,2),AIDEVL(CI, 3) 
BeALDEVILCL, 4) > 


100 CONTINUE 
Aw KAKA KAKA KEK KK KKK KR KKK AK KKKAKAKKKKKKEKE 
C Ove LO Theil "CUChDINATING ELEMENTS IN THE 


DISTRIBUTION 
wee eR KA KEKE KKKKEKREKRKKRKK KKK KKEKKRKEKRKREKEKKRKEEKKKK 
eO READ(S e120) NNT, J,K,R, xX 
120 FORMAT(41, 2F) 
MAGN. EO.0)GO° TO 130 
NE=NN 
CIl=I1 
CJ=J 
RIDCNN)=CMPLX( CI, CJ) 
Z(NN)=CMPLX(R, xX) 
Gor TO. 110 
130 DO 140 IE=1,NE 
Reo 32 ol. Cul 
35 PORMAL( 271 .2F) 
Pome. OyGO LTO 140 
Veo 
CoE SHRCHCNE, RID, I1E,1I,J,1A,11D) 
Z1=0.0*Z1 
Z2=0.0*Z2 
NO 136 IEL=1,I1A 
Ziel oTZ( WilCTEL) )) 
QR= WPM 
136 CONTINUE 
Omer Eib= 14.1 A 
Cina Glennie) = 22) *(i/ZCLIDCIEL)))*CUl 


138 CONTINUE 

MeGrhA.Giel) TYPE 101,.J 
11 PO wOAULPON THERE IS MORE THAN OWE 
SOURCE 


Lee OuER TO ih BUS’ 14) 
140 Cont ENUE 


—— 
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HOw sO l= ie NE 
ees eee i, TENRSH, TUST1,TWST2,TS1 
145 FORMAT(21,4F) 
act iaboso Go TO 151 
Wie see MEAR ID, 1,.1.J,J,1A,11D) 


io 
147 XFRCIID(IA),1)=CMPLX(TS1,TINRSH) 
XFR(LIDC(IA) ,2)=CMPLX(TWST1, TWST2) 
TA=0 
150 CONTINUE 
wWamKRKKKRKKRKKKRKKAEKKEKKKKKRKAKKERREKKEKKRKKKRKKRKKKRKHSE 
é NOW TO DETERMINE WHAT ELEMENTS HAVE LOAD BUSSES 
CONNECTED 
pee KK RH KKKKHARHREKAKRERKERKRRKKKKEKRKEKEKKEKEEKKHREKE 
eel DO 190 I=1,NE 
IR=LFIX(REAL(RID(1))) 
IM=IFIX(AIMAG(RID(I1))) 
IF(AIDEVL(IR,1).NE.0) GO TO 155 
IF(AIDEVL(IMN,1).NE.0)GO TO 156 
EO “16190 
155 LI=IFIX(REAL(RID(I))) 
JJ=IFIX(AIMAG(RID(1))) 
Oe TO. on, 
156 P= TELXCAINAGCRID(1))) 
JJ=IFIX(REAL(RID(L))) 
r57 GnT=11 
aie 


Panny (Tel y=CMPLX(CJJ,CI1L) 
Aipevicie2)=CMPLX(0.0,AIDEVL(11,1)) 
Pen Glew ySeMPLXCALDEVL(II,2),AIDEVL(II1, 3)) 
TPCALDEVL(I1,1)-.GE.35)AIDEV(1,5)=CMPLX(0.0, 
Ieee vinchm. 4) ) 
Per vevwnCrn. 1) ste) ALDEV(L, 7)=CNPLX(0.0,AIDEVL(II,4 


)) 
159 Pare SO al 1 
158 FOKMAT(’ SELECT DEVICE FOR 


Mier, 13, PROTECTING ELEMENT 
PenienGay 8tOAD AT BUS’ ,14,°’..,°$) 
ACCEPT 160,IDEV 
160 FORMAT(C1L) 
Prope v lb bOATCIDEY ) 
DURR=AIMAGCAIDEV(1,4)) 
GUiek—-REALCAIDEV(1L,4)) 
Poo imc mnhy (la )).Gl.0)CURR=10**(DURR) 
CURR=REAL(CUREL(CI)) 
eons oc BC lly, ALDEVL,ALDEV2,AILDEV3,CT1) 
Cail 
aa VC CUnm@mal JJ) LL,ALDEV1,AILDEV2,ALDEV3, 
Peep ores xERCL, 2), RE, EBCJJ),EBCII1),RID,CT1,DEV 











re) 


200 


205 


210 
ELEMENT’ 


Z30 


EV ) 


Zz 20 


730.0) 


EV ) 


500 


400 


410 


sas 


Gil =0.. 0 
ieGeketEs09Ge TO 159 
CONTINUE 
IFEL=0 
NEL=NE-IEL 
IF(REAL(AIDEV(NEL,1))-NE-0)GO TO 300 
INDEX=1 
I=IFIX(REAL(RID(NEL) )) 
J=IFIX(AIMAG(RID(NEL))) 
Prune ALDEV(NEL,2)).NE.0)GO TO 250 
i hee ess NEL 
RE=0 
FORMAT(’ ENTER THE DEVICE FOR THE BUS’ ,13, ‘OF 
» 14) 
NG@neT 290, 1DEV 
FORMAT(1) 
PGubEy SEOs 100)GO TO 400 
AIDEVI=FLOAT(IDEV) 
CURR=AIMAG(CUREL(NEL) ) 
SS npeceipen( CURR. EBC ),AIDEV],AIDEV2,AIDLV3,CTI1) 
Piao yAIDEVCNEL, 5)=CMPLX(0.0,.414973) 
Pom ine COURE I NEW, J,AlDDEV!,AIDEV2,AIDEV3, 
LAIDEV,XFR(NEL, 1) ;XFRCNEL,2),RE,EB(I),EB(J),RID,CT1,D 


MiCkE. NE. O)GO TO 205 

ime ee 22001, NEL 

MECEPT 230,1DEV 

AIDEV1=FLOATCIDEV) 

RE=0 

CURR=REAL(CUREL(NEL) ) 

SH CHUNGnCCURR EBC I), ALDEV1,AIDEV2,ALDEV3, CTI) 
IF(AIDEV1.GT.34)AIDEV(NEL, 5)=AIDEV(NEL, 5)+CMPLX( 4149 


Cece Cia NEL, lyg,ALDEVi,ALDEV2,AILDEV3, 
Paavo richie |) .,XFRCNEL,2),RE,ERC1),EB(J),R1ID,CT1,D 


PCB NE.O)CO TO 250 
LEL=IEL+1 

ieananelLt.RE) GO TC 200 
IFC INDEX.EQ.0)GO TO 400 
IEL=0 

INDEX=0 

COmEO! 200 

TYPE 410 

PRINT 410 

MOmeinr(  HEVICE LISTING’ ) 
DO 440 I=1,NE 

Tvemee 2091, (ALDEV(1;3),J=1,7) 
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Dotcom. (ATDEVCL, J), J=1,7) 
iL ie@e 20) 1 CARDEV(1,J),J=1,7) 


420 FORMAT(I1,14F) 
440 CONTINUE 
WRITE(37,442) 
442 Dr we sO OnOee O10 0-0 0.0 0.0 0.0 0.0 0.0 0.0 
0.0 


lO OO O..0° ) 
2000 CLOSECUNIT=30) 

GEOSECUNIT=31) 

CLOSECUNIT=32) 

CLOSECUNIT=33) 

CHOsECUNIT=34) 

CEOs e CUNT T=3 6) 

Chose CUNT T=37 ) 

SLOP 

END | 
Hee He KKK KKK KK RK KKH KKK KKK KKK RK KK KK 
doko eK TK KK KKK RAK AKA AA AK KA RAKE 


C see OUTIINE SU RCHSFINDS THE PLEMENT CORRESPONDING TO 
TWO 
C BUS suiio. 


RARKKK KKK KHER KKK KKKEKKERKE KKK KKH 
SUBROUTINE SERCH(NE,RID,NN,1,J,IA,IID) 


DIMENSION RID(100),11D(3) 
COMPLEX RID, C1, C2 





LA=0 
DO 30 PE —=is 100 
CJ=J 
Cl=1 
Cl CHUPLACC I CI.) 
C2-—CMPT KCCI SCL) 
PACD ln Ly .kO.Cl)y7 CO TO 20 
PGC) ioe hOlne2 GO FO” 20 
GO TO 30 
20 A= 1-2 | 
LIDCLA)=IEL 
30 CONTINUE 
5.1 REeUR 
END 


RRKAEKEKEKKKEKKEKEKEKKEKKEKKEKAKRKKKKKEAARKEAKRKKRAKKAKKRKK KKK KK 


C OUP NE wbeAbT READS FCURSDAT TO OUTPUT FAULT 
SURKENTS 
C Pee Ur getty Sino 


cin KKK aK ARR KAKA KAKA KR KAKKAHEKKR KKK KEKEKAKKKKEKEKEKSE 
SHOUT INE REFALTCIPALT. NEL, FALT, CUREL) 

10 REeiea ail 5) LEF,NN,TFAL, SFAL 

15 FORMAT(21, 2F) 





i25 


iGrerenoROyGO TO 20 

PaCier sibetEALT)CO TO 10 

BE(IN. NEeNEL)GO TO 10 

FALTI=TFAL 

FALT2=SFAL 

FALT=AMINI(FALT1, FALT2) 

FLOCUR=3.0*CUREL 
IF(FALT.LE.FLOCUR)FALT=AMAX1(FALT1, FALT2) 
WALISAMAX IT CFALT1,FALT2) 


20 REWIND 33 
Re PUR N 
END 
POOR RRR KAKA KK KKK KKK KARE KKARAKKKEKKKKKEKKEKKKRKHREKKRKK 
c TCU itt EDDY Sok iS THE DEVICE INSTANTANEOUS 
Peles EN T 
C eNO NGwON THE FAULT CURRENT AT THAT BUS 
BOES 
. COCR DIMATE With OTHER DEVICES 
MTom KAMARA KH KKK KA RRA KKK KKKKEKKKRHKKRKKEKKKKAKRK KK K 
SU SROULTINE 


ey ( CURBL,IEL,I,J,ALDEV1,AIDEV2,AILDEV3,ALIDEV 
eee. KER? RESE BI EB2,RID,;CT,DEV) 
DIMENSION 
mm 100,7) ,CURELC100),DEV(20,0/4),RID( 100) 
Ole cen pe Vv. XERI,KFR2, CUREL,RID 
Poe CeIDEV1) 
3g) embry 1.1.0 )CALL 
Seen, CUREL,IEL,I,J,ALDEV,AIDEVI, 
eee ey LDEVS.xXFRISXPR2,RE,EB1,EB2,RID,CT,DEV) 
Mimenmieylshr.0)GO TO 100 
Pacem VvisekO.99)CO TO 100 
Meeniny Woh e Nl. [DEV GCr.QO)CALL SETFUS( CUREL(CIEL 
PreenrlhiereDEve Div ALDEV2 ,ATDEV3,XFRI,XFR2,RE,EB1 
eb, 
DEV OR ID’) 
Bre visGh «3 5.AND. LDEVI.LT.99)CALL 
See kK BC CURELCIEL), 
Paeeat be eDEV. TNR V I. AILDEV2,AIDEV3,CT,EBI1,EB2) 
100 RETURN 
END 
Mn K KKM K KKK KKK KKK K KKK RAE KARR KEKE 


Osc Oh cete eUs ct iv cde ents ale cle, als. als. als. x85 abe 


C SU etiEeeot ir USS COORDINATES AND SELECTS A FUSE 
KA KKK KK IKK KKK A KKK KKK KAKI KKK KKK KARR KKK KAKA KH 
SUBROUTINE 


Pembis( CUMBMel. J,LEL,AIDEV,IDEV],AIDEV2,AIDEV3 
Peo oe EB!  ER2 DEV, RID) 
DIMENSION 

Pager (100,6), CUREL(C 100) ,DEV(20,0/4),R1NC1]00) 





20 


25 


35 


40 


43 
45 


E26 


Goer eal DEV, XPRIn x FR? , CUREL, RID 

Of= LP ica MAG(AIDEVCLEL, 2))) 

Pelee) CALL ME DEV(1, J, 1EL,ALDEV, COOF) 
Beebe) CALI xhDEV(2,J,1EL,AILDEV, COOF) 
eo Oo 1CURR=ATMAGCCURELCLEL)) 

Pi JieNEDO )CURR=PEAL(CCUREL(IEL)) 

MiG erie ANl.COOF.LI.«cCURR)COOF=CURR 
IF(JJ~NE-0)COOF=(EB2/EB1)*COOF 
IF(JJ~-NE-O-AND.COOF.LT.CURR)COOF=CURR 

PeGnpeEN )2nb. 15 CALL SFUSE(COOF,IDEV1,OFSET,OFSET1) 
Revive isl 5) CALL eSBORPICCOOF, IDEV1!,OFSET,OFSET1) 
Pacisi.2o.0) GO TO 20 

C1l=I 

Cn 

Pepe Vir ?, 1 )=CMPLX( CI, CJ) 

Pie Viele? )=ATDEV CLEL. 2)-+-CNPLXCFLOATCIDEV1),0.90) 
AIDEVCIEL, 3)=CMPLX(OFSET, OFSET1) 

ComroO 25 

Cl=I 

CJ=J 

Pa Ey Creel = CMP LX CCl. Cy) 

AIDEV(IEL, 2)=CMPLX(0.0,FLOATCIDEV1)) 

Nebel nh. 6 )=CNPEx OFFSET, OFSET1) 

Oo TOCSO6 

PeaGCceBSCxPRi).EQ.0)GO TO 50 

AX=ALOG1LO( (CAIMAG(XFR1))/3.0) 

Bauer CN TDEVCIEL,2)).LE~.15)GO To 35 

IFC OFSET.LE.AX )CUR=10**AX 

PGor Sh tebheax) CALL SBORFUC CUR, IDEV,OFSET,OFSET1) 
COr 10 40 

IF(OFSET.LE-AX)CUR=10**AX 

PeG@OrGH ©. DE SAX) CALL SFUSEC CUR, 1IDEV,OFSET,OFSET1) 
ene eC rie 2 = COMPLY (PLOATCIDEV),AIMAG(CALDEVCIEL,2))) 
AIDEVCLEL, 3)=CMPLX(OFSET, OFSET1) 

AR=ALOG1O(REAL( XFR2)) 

Cpr uome2 REALCALDEVCIEL, 2)),AR,OFSET1,TVAL1) 
ivecin= 1 Oe = Py AT | 

BG AleGi=2.0)GO TO 43 

GOe TO. 50 

TYPE 45 

FORMAT(’ UNABLE TO PUT FUSE IN THIS LOCATION 


SupsSTLITUE 


50 


LANOTHER DEVICE’) 

RE=1.0 

Cor 10) 400 

GwhetheAnptCd., LEL.PALT, CUREL(IEL)) 
PAUSE 55 

erMeg hee @) GALL 


ewe were CURE, RID, TEL, 1, J,1TVAL1, DEV, 








DID) 


Z00 


250 


300 


ie 


CVA) 

PPememeeinoALL, SEAL ?(AIDEV,IEL,1,J,TVALI1,DEV, 

eA GVAL) 

PAUSE 60 

ADEV=FLOATCIDEV1) 

Genus ec ll ADEV,CVAL,OFSET, TVAL) 

MaGill. TVALI)GOe TO 200 

ADEV=ADEV+1 

NEGDEV=-IFIX(ADEV) 

IF(ADEV.LE.15)CALL SFUSECCOOF, NEGDEV, OFSET, OFSET1) 
ene Gl ils) CALI SbORFU( COOF ,NEGDEV,OFSET,OFSET1) 
COTO 55 

coe 0) GO TO 250 

MiG Eine? j= GipLXx (0.0, ADEV) 

AIDEVCIEL, 4)=CMPLX(OFSET, OFSET1) 

SO THON SHOo) 

MDE i PEt 2)=CMPExX( ADEV,AIMAGCALDEV(CIEL,2))) 

Piney Gli. 3 )=CMPLX(OFSET,OFSET1) 

RETURN 

END 


PA Cm RRR KKRAKREKAKARKEKKKAREKKAKKKKKKKAEREKKKKEK 


Pea RAR RARER AEKAKKAKKKRRARKKEKAKKARKKKKKEKEEKX 


C 


Cue RENT 


POU ibe ruSE COMPUTES TIME VALUE FOR A GIVEN LOAD 


Mem KKEKKKKKKEKRKAAKRKARKRKKRKRKKAKKRKKERKKKAEREREKEKE 


eo'*3 ) ) 


40 


*3))+ 


BO 


7 


SUBROUTINE FUSE(ID, ADEV,CVAL, OFSET, TVAL) 

TVAL=0.0 

IDEV=IFIX(ADEV) 

[nar ny sere 15)CO TO 50 

X=CVAL-OFSET 

IF(ID.NE-1)GO TO 40 

eee Ot = 36.42 94*%)+(223.7526*(X**2))+(-751.279% 


alee emo eee) aa 222 0.67 *(X**5))4+(440.813*(X**G) ) 
Py TVAL—2.0 

GO TO 100 

Sev AL-ORSET-+. | 

ieee ora S722 se 21 909 eC X**2 ) )+(-213.420*%(X* 


1¢6201.771*(X**4) )+(-96.9836*(X*¥*5))4+(18.6046*%(X**6) ) 
i we=TVAL=2 20 

ie Chala )VAL=—1.95 

GOmtO 100 

X=CVAL-OFSET 

Meith) GO TO 70 

ete 7 oUGt(—4. 26275654X)+C(1~-219340*(X**2)) 

Wie hVAL=2 20 

GOmnC  10G 

X=CVAL-OFSET+.1 





LZS 


ieee nico (= Gag 224*X )4+(3.23181*(X**2))40=.60571 
36* (X**3)) 

iVAL=TVAL—2.0 
100 RETURN 

END 


PEmIO RA KR KRARKKKKEKKKAKKKEKKEKAKRKEKKEKKKKKEKRKEKKKKKEK SE 


KHEAEKKKKKKKKKAKRKKEKEKKKKKKKKKRKEKKKREMKKKAKKKKKKEKKRKKEEKE 


C PP VOUT THE KELAY COMPUTES TIME VALUE FOR A GIVEN 
CURRENT 
C Voter Olen ec re tC RELAY 


KRKEERKKKKKKKEKKKEKKKKKR KKK KEK KEKKKKKKKRKKRKK KAS 
SUBROUTINE RELAY(NDEV,CURR, TVAL, TDSET, DEV, AMST) 
DIMENSION DEV(20,0/4) 

TVAL=0.0 

X=ALOGIO(CURR) 
Mme eh. 1.2)X=1 62 

)O 10 I=0,4 
TVAL=TVAL+DEV(NDEV,1L)*(X**I) 

10 CONTINUE 
IF(NDEV.EQ.7) TVAL=TVAL-2.0 
TVAL=TVAL+TDSET 
MGC eOGT.ANST. AND. AMST-NE-O)TVAL=-2.0 
RETURN 
END 


Kee AMR KAKKK KKK KKKKAKKKKEAK KK KKK KEKRKK KAKI K 


e eon wekhI BER COMPUTES TIME VALUE FOR A GIVEN 
CURRENT 
C VALUE 


Prac kc nam KKK aK KKK AKA KAKKER KAA K EKA KKK KKKEKKKEE EK 

SUBROUTINE CKTBKR(ID,MDEV,DVAL, OFSET, TVAL, AMVAL) 

Sab —=DVAL-OF SET } 

BGhistEsl)GO TO 50 

CVvb=CVAL+. 1122 

IF(AMVAL-EQ.0.0)AMVAL=10 

MACCVAL.Gi.AHVAL)GO TO 90 

eee lot — Oe oo2 oO ( CVAL) )+(23.65419* (CVAL**2) )+ 
34.7922 

ee eee ode fy 2 OCC VAL* <4 ) )+(-21.14536% (CVALRS5 


fe 6.4 
mols oh CVAL**6 ) ) 
TVAL=TVAL-2.0 
GO GO. 100 

a0 CVWwe-CVALH . 2315 


AMV L=AUVALS.07 5 

MEGevALeGreANVAL)GO TO 95 

TVAL=5 -8596722+(-4.394335*CVAL)+(- -50690284* (CVAL**2 ) 
me 10.40 

12739* (CVAL**3) )+(-17.064956* (CVAL**4 ) )+€(13.011543*(C 
VAL**5)) 





—- 





re 


90 


aD 
100 


Zo 


eee O07 24% ( CVAL* 6 ) ) 
ve = by AG 2. 0 

coro hoo 

TVAL=—-2.0 

GO] TO 100 
TVAL=ALOG10(.025) 
RETURN 

END 


RHKEEKEKEKKKKEKKKRKEKKKEKKEKRKKKERRKRKKREREKKKKKRERKEKEKES 
KAEEAKKKAEKKKKEKKKEKKKKKKKKEKKREKEKAKEKKKKAEKKKEKRKKRKKKKKK KK 


6 


oul ENE POETREL;COORDINATES AND SELECTS A RELAY 


KREEKKKRKEKKKRKKAKEKEKRRKRKKRRKKKRKKKRRERKEKEREKKEREKKKERSE 


SUBROUTINE 


feeeeeeCUREL, LEL,1,J,AIDEV,ALDEV1,AIDEV2,AIDEV3, 


oR, XFR2°RE,EBI.EB2,RID, CT, DEV) 
DIMENSION 


melOO, 7), LIDD(30),RID(100),11IDC3) , CURELC 100) 


foe 0° ) 


ibe. 2020/4),FTC30; 2) 

CO lehEX ALDEV, CUREL, XFRI,XFR2,RID 

GOOrF-(0 .0 

CUR R=@ .0 

GhuUR-O - 0 

PvP E> 

Bema os OnSe THis RELAY HAVE AN INST. ELEMENT,IF 


Reece e tT (7 1b 

FORMAT(CL) 

oe i XCAIMAGCATDEV(CIEL,2))) 

MDEV=LFIX(-ALDEV1) 

Petes. O ibe v LELe | )-CHPLX(FPLOAT(L),FLOAT(J)) 
Miri OO par DEVCLEL.2)=CHPLKX(O0.0,AILDEV1) 
Beet ib volt la) —CMPEXCALDEV],O.0)+ALDEVCLEL, 2 


Cee Geo CAbre DEV Cl, J,1EL, ALDEV, COOF) 
Muerieveew) CALL AFDEV(2,J,LTELYALDEV,COOF) 
Petre 60) GURR=AITMAGQGUREL(LEL) ) 

PGi. Neo) CURR=REAL( CURELCLEL) ) 

Pete: Ora. COOH. Lt. CURKR)COOF={CURR 
Duel hee eOOr—(EB2/EBI)* COOF 

CH LOG lo CCOOF ) 

miei OeA ND. COORwLT. CURR)COOF=CURR 
PRGCDUR.GT<ALDEV2)AIDEV 1=4.0 

(imecUR CG hel DEV 2eAN De Js bO.9) CALL 


SELECT(COOF,EB2,AIDEV1, 


LALDEV2,AIDEV3,CT) 
Py GenurMecneAlDEV2eANDsdJ.NE.O)CALL 


Pep CT(COOF,EB1,AIDEV1, 


PAIGE 2 AEN CT) 
CUES ATMAGCCURELCIEL)) 





wy) ) 


130 


Meteor DE VOrh les) —OMPLX(ALNDEV2,AIDEV3) 
Peano it lin) —ClPLX(ALDEV2,ALDEV3) 
Piel toNce Oy CURE-REALCUCUREL(ILEL) ) 
oe ere Gini —CUPLX(CLT,0.0)+AIDEV(CIEL,7) 
cies Alley Giri. 7 y—CMPLX( 0.0, CT) 
ie ier) ) CO BOs 00 
Ceo Noto), Leer ALa),CURL) 
aG@iett.0)GCO TO 32 
PA HAL *( EB1/EB2) 
Pil eb ily (hOes( REALCALDEVCIEL,4)))) 
IN=0 
Dem l=). 30 
Ceoeeor Renata tite), JL, 1A, 11D) 
PGi. kos eo. LTO 10 
DO By LIAS 7. ILS 
IN=IN+1 
epeKN y— Ie Di LA) 
CONTINUE 
CONTINUE 
Doms0 LYN=l, LN 
Lo=11l pir 
PCr wsonlEL) GO TO 30 
IJ=IFIX(AIMAG(RIDCID))) 
CURJ=AIMAG(CURELCID)) 
Ceee aL rcls 1D, FALT, CURJ) 
PvP 2=AbOGlOCrALT) 
MR eAmeAlLDEV( 1D, 2) ) 
ie le GA 11) ) 
Pol wee GE.0)GO TO 16 
Pewee EALTO=FALT/(10**( REALCALDEVCLD, 3 


NDEV=-I1FIX(REAL(AIDEV(LD,2))) 

KEAL=IFIX(10*FALTO) 

DOM to eK=15 KEAL, 10 
RK=K*. 1 
PaGMMOGIORI@rsersl.2)GO TO 16 
Gib 


fee y (NDEV, RK, TVALI,REALCAIDEVCID, 6) ) 


PebiYy Reel ALDEVCLD,5))) 
eae 


MAY (MDEV,RK,TVAL2, TNSET, DEV,0.0) 


15 
16 


HimO-( LOssTVALI)+.4 
TLIMER=TVAL2 
TIMES=ALOGIO(TIMEO) 
iinet l—=TINES-LIMER 
PeG@@ist l1l.abE.0)CGO TO. 15 
Dani + ENOL T 

CONTINUE 

Mi@reip. bas )GO £O 25 





Je '3 


OF SET=AIMAG(AIDEVC(ID,3)) 
FNSET=FALT2-OFSET 
BaGenNokbecieie os )FALTI=OFSET+1.5 
KVAL=LFIX(10*OFSET) 
KFAL=LFIX(10*FALT2) 
DO 20 K=KVAL,KFAL 
RK=K*.1 
[eC Al. 1.35) OMe esh< 2. ALD. RK,OFSET,TVAL1) 
Ma G@AiDiecin. 35) CALI, 


CKTBKR(2,1FIX(ALID), RK, OFSET, TVAL] , REAL( 


)) 


AWG VLE I. 5) ) ) 
CUR-ClUS=RK)/ (1O0**C(REALCAIDEVCLEL, 3) ) 


CALL 


RELAY(MDEV, CUR, TVAL2, TMSET, DEV,0.0) 


20 
Ze) 
30 


Bed 


TIMEO=10**TVALI+.4 
Pe R= Ty AL 2 
Pies -AbOGNOC TIMEO) 
TNSETI1=TIMES-TINER 
Paewish tiene. O)GO TO 20 
TMSET=TMSET+TMSET1 
CONTINUE 
Cor le 30 
CONTINUE 
lon Clemo CuUPLX (0.0, IMSET) 
hick T=0. 0 
Comino 45 
wip A bi AGCATDEY ( LHI;, 2) ) 
FALT2=ALOGIO(FALT1) 
Perl h-Cl.0 GO 10 35 
PAL hiakALEUA(CEBI/ EB) 
Rae wwe ie) FALT=FPALTI/ (104s (REALCALDEV(IEL,4)))) 
NDEV=-LFIX(AIMAG(AILDEVCIEL,2))) 
KFAL=IFIX(FALT#10) 
Diet FAL—15 
ID=IFIX(DIF/50) 
Pen iO) 1D—1 
Domes oe K— oe KEAL, LD 
RhKen- . 1 
Mie RK«~CEBZ/ERI)*(10**(REALCALDEVCIEL,4)))/10** 
1(REALCAIDEVCIEL,3)))) 
Ghee ta 20)GO 10.34 
CALL 


PomaAY (NDEV, RK, TVAL1,AINAGC(AIDEVC(IEL,6)),DEV, 


PraG ATDEYV CIEL, 5 ))) . 
Cierny. RK2,TVAL2.TMSET,DEV,0.0) 
TIMEO=10**TVALI1+.4 
TENE S=ALOGIOCTINEDO) 

TUR LV A lL 2 








5 
34 
oy) 


e532 


TMSETI=TIMES-TIMER 
MGnISEmieLheo)GO TO 33 
fiSer=TNSETICIMSET 
CONTINUE 
AIDEVC(LEL, 6) =AIDEV(LEL, 6)+CMPLX(TMSET, 0.0) 
TMSET=0.0 
MeCATMsLT.0)GO TO 45 
OF SET=AIMAG(AIDEV(IEL,4)) 
FNSET=FALT2-OFSET 
LF(FNSET.GT.1.5)FALT2=OFSET+1.5 
KVAL=IFIX(10*OFSET) 
KFAL=IFIX(10*FALT2) 
DO 40 K=KVAL,KFAL 
RK=K*.1 
mG ess) @AlIemhUGE( 2. AlD.RK,OFSET,TVAL1) 
ME CGAUDeCR. 35) CALL 


CKTBKR(2,1FIX(AID),RK, OFSET, TVAL1, AIMAG( 


40 


45 


PeCEALT1 , 


100 


leamE Vv ( 1 BIL. 5) ) ) 
CUR=RK-REAL(CAIDEV(LEL, 3) ) 
CURTA=10**CUR 
Quran, Epil) ~( 1 028 Re ) 
CuRis—cURy (1022 (REALCATDEVCIEL,3)))) 
MaACeUrn bw. ls5)CO TO 40 
vier ea CMbDEY CURTA. TVAL2, THSET, DEV,0.0) 
TIMEO=10**TVALI1+.4 
PIMER= 1 VAL 2 
ie eoeOcwoc TIMEO) 
TMSETIL=TIMES-TIMER 
MoGhiCek Lhe eh.0)GO TO 40 
TMSET=TMSET+TMSET1 
CONTINUE 
Pepe Cl EiaG =ALNEVYC CIEL, 6)+CMPLKX(CTMSET,0.0) 
TMSEL-0. 0 
Ee. bOsOeAND. FALL? .CT.6)FALT1=6 
Pein Omori J) 0.0) ALDEV( LEL,5)=CMPLX(0.0, FALT1) 
Pett! Ady ClO**REALCAILDEVCIEL,3))) 
Pace Abri.Gl.o)FALT1=6 
Pan how An DaJJ.NE.O)ALDEVCIEL, 5)=ALDEVCIEL, 5)+CMPL 


Ore0)) 
RETURN 
END 


KRAKAKKEKKAKEKAEKKAKEKAKAKREKAKRKRKKEKAKKEKKRAKKERKKAKEKE SE 


C 
BOWCEST 
C 


Sew k sides PROVIDES THE DEVICE WITH THE 


ities beep owi Sr TRok DEV LN ELEMENT. 


KRAKKEKKAKKREKKEKAEKEKEKAKKEKKKAKKKEKRAKRKEKREKRKKKEAKRKKKK 


SOU TENENSEALTI(AIDEV, CUREL,RID,IEL,I1,J,TVAL, DEV, 
1ICVAL) 








oS 


DeeiomonesmEDEYOlO0O,6),CURELC100),RIDC100),11D(3), 
erp coo). DEV (20,074) 
COMPLEX AILDEV,CUREL, RID 
Pore J1=), 30 
Cambie cE RCH! RLOMmea, J,J1,1A,11ID) 
BeG@uA. bio.) GO TO 10 
DO. 8 JA=1,1LA 
IN=LIN+1 
Lipper Ny=l LDC LA) 
8 CONTINUE 
10 CONTINUE 
DO 30 INN=1,I1N 
ID=ILIDDC INN) 
P= 1 Pe CALMAGCRIDCID) ) >) 
Clini. MrAG CCURELCID)) 
Corie eee lejel Dp FALT, CURS ) 
FALT2=ALOGIO(CFALT) 
Siva Cun mu 1D, 2 ) >} 
Peep ir Ex Cal D> 
BED «Cle 32)CG0O TO 100 
Evie Rea Al Dby Clb, 4) )+. 5 
FALTO=10**(FALT2) 
pee sea) FAL T=FAIL1TO/ (10** (REALCALDEV(CLD,3 
yD) ) 
Niev=—IP P< CREALCALDEV(ID,2))) 
Paes ann CALL RELAY (NDEV,FALT,TVALI, REAL 
mew CID soo)  DEV,REALCAILDEV(ID,5))) 
PMR EAGbCALDEV( LEL,4))+.5 
WmGr AL. LI.kKALT2) FALT2=F AL 
iG pael. Oo )CALL 
meee 2,Al1D,FALT2,AIMAG(ALDEV 
Gleb 3)))) pov Ie) 
PG winch. EVALIYGO TO 25 


GO TO 30 
Zo Pieri. ble Oo )CVAL=—FALT2 
CVAL=FAL 
TVAL=TVAL] 
30 CONTINUE 
100 RETURN 
END 


KRKKKAAKKEKAKRKKKKKEKKKEKKKKAKEKRKKRKKEKEKREAKKRAKEEAKE 
KKEKAEKKEKKEKEKKKKKAKKAKREKEKKRKRKEKREKKEKKERKARKRKAKRR KRHA KK HK 


C SO ier AS SIN  1,COMPUTES TIME FOR LOAD 
mene OF 
C BREMEN T AT FALT. 


SeEKKERKKKRKEKKKKKKRKKREKERIKEKRKKKKKREKRKRKEKHKKEKEKSE 
Uo mito Ann aCATDEV,IEL.,1,J3,TVAL1,DEV,FALT, CVAL) 
DIMENSION AIDEV(100,6),DEV(20,0/4) 
COM ALOEY 





ea 


FALT2=ALOGIO(CFALT) 
AID=AIMAGCAIDEVCIEL,2)) 
FAL=REAL(CAIDEVCIEL,3))+-5 

FALT=10**( FAL) 

Woto-P Ali ( bO** CAIMAG(AIDEVCLEL,3)))) 
iG@een. bt. 0 CALI. 


Remy CNPEV, FALT3,TVAL! , ALMAG(AIDEV(IEL 


1,6)),DEV, AIMAG(AIDEV(IEL, 5))) 
MmaGmAbe hb). FALT2)FALT2=F AL 

Mnemneen.d )CALL hUSEC2?, AID, FAL, AIMAG(AIDEV(IEL,4)) 
ieavAL 1) 

CVAL=FALT2 

RETURN 

END 


ARRKKKKAKKKKEKKRKKKKAKKKEKKKKKKEKAKKAKRKKRAKKEKKEKKKKKKEKR 


HRHAREKKKKKKEKKKAEKKKKAKKKEKKAKRKERKRKEKEKKKKKKKEKKEKKEKKK 


30 


a2 


EV ( 


DEV ¢ 


40 


50 


SUBROUTINE XFDEV(NN,I,1IEL,ALIDEV,COOF) 
DIMENSION I1D(30),AIDEV( 100, 6) 
COMPLEX AIDEV 
COOF=0.0 
IX=0 
BEGUN. NEST )GO TO 50 
Doms02 LA=1.,100 
AILD=REAL(CAIDEVCIA,1)) 
mem Ix (AID ).NE.L)GO TO 30 
IX=1X+1 
LIDCIX)=IA 
CONTINUE 
Om oc K=1>51xX 
CONTINUE 
DO 40 JA=),1X 
Cmote— Ab MAGCA LDEV( LLDCJA), 3) ) 
PacnealCALDEVCLID( JA),2)).LT.0)OFSET=REALCAID 


CTA, 3) 
PiGwe wevE VGC pC pA), 2)).GE.35)0FSET=REALCAIL 


rela) = 5) )+0FSET=.05 
PieweNnG@wlUEVC(IEL,2)).LTr.0.AND.REALCAIDEV(II 


ewe Or GlelernieREALCALTDEV(LIDC JA) ,2)).LT.35) 
ZO SE I-ORFSET+. 15 
[PeCCOreGieQFSET)GO TO 40 
COOF-OF SET 
COMELNUE 
GO TO 60 
Pee we ee serie 2) ).6CT.0)COUF=AJNAGCAINDEVCIEL,4)) 
ewe Cet vty 2)) bl. YCOOF=REALCALDEVCLEL, 4)) 
PiemrtencwOmY CLE 2) ).GE.35)COOF=AIMAGCAIDEVCLEL, 5) ) 


(eS OR ee EE SESS SSF OR We 





5 


-.075 
1+COOF 
Pee eweorey lil? oat. oAND.AIMAGCAIDEVCIEL,2)).GT 


eo eG tory lib ps Ll. 35)COOF=COOF+.15 
60 COQOT=COOR+. 025 

DOOF=10**( COOF ) 

COOF=DOOF 

RETURN 

END 
KKKKKKKKKKEKEKKR KERR KERK KR KKK RR KR KK KKKKRKEKKEKKKK 
KKK KKK KERR ERE KKK KKK KKK KKK KEKE KKK KKKKA KEK 


C Ut twom Cm CORUINATES AND SELECTS A THERMAL 
MAGNETIC 
C BREAKER. 


Pca KKK KA KKKKK KKK KAKA KKKKKKKEKKKKEKKEKK 
HRM KRKKKKKKEKKKKKEKKEKRKKKMCKKKKKEKAKRKKRKEKHRKKEKKREEH 
Soro ShdChE CUREL, LoJJ,1TBEL,AIDEV, LDEV,OFSET 
On Sh iieANAG, EBL. EB2 ) 
DIMeENSEON AILDEVC100,7),CONC5) 
COMP LEX -ALDEV, CUREL 
CON(1)=2.6 
COG? y=3.. 5 
CONC3)=5.0 
CON(4)=6.0 
CON(5)=7.0 
Mee CAITMAG( AIDEVCIEL,1))) 
Pee O yOnkL XFDEV(l, J,1KE,AILDEV, COOF) 
Pere. O)CALL XFDEV(2, J, IEL;ALDEV, COOF) 
DiGi O.0) DOOH—-ALOGIO(C COOF)-AIMAGCAIDEVC(IEL, 5)) 
IF(J.EQ.0)COOF=10**DOOF 
imei via )DOOP=ALOGIO( COOF )-AIMAGCAIDEV(CLEL, 5)) 
Paiedeet ie 0) COOF=10**pDOOF 
IF(J.EQ.0)CURR=AIMAG( CUREL) 
Pee N EO McCuRR=REALCCUREL ) 
GG blest R2)GO TO 45 
MaeGaogr abt. CURR) COOF=CURR 
Cire LO.46 
45 RE2=(EB2/EB1)*COOF 
ue 2 Gh CURR) COOF=RE2 
Peeve 2k t.CuRR PCOOF=CURR 
46 Cee obiekCGOoOr. 1DEV3,O0FSET,OFSETI,DMNMY) 
DVO E Ver r.IDEV3)LDEV=IDEV3 
PeGIpEV.GE. LDEV3)ITDEV=LDEV+1 
Cara 
Ca Jk) 
Piva bel y=CriPLXCClL, CJ) 
Paes iii. LEL. 2)=CNIPLX(O0.0;,FLOATCIDEV) ) 
Peer lDey Cli, 2)=CNPLXCFLOATCIDEV),0.0)+AIDEV 
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WeLEL, 2) 
IF(J.EQ-O)AIDEVCLEL, 4 )=CMPLX(OFSET, OFSET1) 
[Monee PATDEVC(LEL, 3)=CNPLX(OFSET, OFSET1) 
IF(J.EQ.0)FLT=4* AIMAG(CUREL) 
LF(J.NE.O)FLT=4*REAL( CUREL) 
SET=10**(OFSET-.12) 
DO 50 K=1,5 
TRP=CON(K)*SET 
Gu Pac womeerNeoO LO) 60 

50 CONTINUE 

60 LF( AMAG. NE.0)AMAG=ALOGLO(CON(K)) 
IF(J-EQ-O)AIDEVC(IEL, 5)=CMPLX(0.0, AMAG) 
LF(J-NE.~O)AIDEVCIEL, 5)=CMPLX(AMAG,0.0)+ALDEVC(IEL, 5) 
RETURN 
END 

Kk kee KH KEKE KKK KKKEKEKEKAKKKKKHREKE KKK 

Swe K KKK KKK KEKE KK KKK KKKKKEKKEKKKEKRERKEKREEKEKKH 

e SUBROUTINE SELECT;SELECTS A PARTICULAR DEVICE 

Peo KKK KKK KRKRKKKKKEARKEKKEKKKKKKEKEKRKKEKKEKEKEKKERKREKEH 
SUBROUTINE SELECT(CURLO,EB,ADEVL1,ADEVL2, ADEVL3, ct) 
IDEVLI=IFIX(ADEVL1) 

54 [PiGRDEVEl nO. oO eo TO 99 
BGmnE yi. bO. 1)GO) TO 60 
[PVGENEViteEO.2)GO TO 70 
PrGrprV Lis bO.4)CO TO 80 
HOG@UmEVbleGr.4)GO TO 53 


GOO EO 58 
sO ver 5:5 
a> Hotwire NOmoUCH DEVICE OR DEVICE ERROR, RE-ENTER 


[DEN LCR NR 1-422 2's) 
MOCHPE 56, IDEVLI 


hg) BORMAT( I ) 
GO TO 54 
53 Be@ebeGh.600-.0) GO TO 53 


Ce Wwieoh RMR CCURIO, TDEVL1,ADEVL2,ADEVL3,ADEVL4) 
ADEVLI=FLOATCIDEVL1) 
mri Vv 14 
COe GO 100 

60 Gwe use ( CURLO, IDE VL], ADEVL2,ADEVL3) 
ADEVLI=FLOAT(CIDEVL1) 
Como 10 0 

70 Di@rurLE.500)GO TO 53 
Cpe GhOREUCCURLOIDEVLI.,ADEVL2, ADEVL3) 
ADEVLI=FLOATCIDEVL1) 
GOO, 100 

80 Creo ele eCURIO. LDEVI],NIAP,NCT, PKUP,CT) 
ADEVLI=IDEVL1 
ADE 2=PKUP 
ADEVL3=NCT*100+NTAP 





iy 


GO TO 100 
2), ADEVL1=99 

ADEVL2=CURLO 

ADEVL3=CURLO 


OD CONTINUE 
RETURN 
END 
BRR KKK KKK RK KER EK KK EKER AKKKARKKERRKKKKKKKREKEKSE 
C Ovo h vot lot wet O SELECT A CURRENT LIMITING 
FUSE 
C Poo ere eee nek HIGH OR LOW VOLTAGE 


WHKKEEKKEKKEKKKEKAKRKKKEKRKKKKKKKKKEKEKEKKAAKKEKERKEKRKEKEKERKRREEEEES 


SUBROUTINE SFUSE(CURLO, NDEV, OFSET, OFSET1) 
MAGNDEY.LT.0)GCO TO 10 
DEV=-11.15221+(23.865*ALOG1O(CURLO) )-(16.487*(ALOG10( 


CURLO 
1)**2))4+(2.9288%*(ALOG1LO(CURLO)**3))+(3.19672*(ALOG1O( 
CU 
Pro yeeayy—-(1.46877*(ALOGIOCCURLO) **5) )+(.1776431* CAL 
OG 


310( CURLO)**6)) 
NDEV=IFIX(DEV)+1 
mew Devecr.15)NDEV=15 
10 DEV=ABS( FLOATCNDEV) ) 
OFSET=.9532-( .N46433*DEV)+(.117885*(DEV**2) )+(-.02935 


2 CDEV 

n= *3) +0 .0033/*(DEV**4 ) )+(—-.0001847*%(DEV**5) )+(.00000 
oo ( 

AW c#* 6.) ) 

Oboe DI—OFSET+. 1 

PeGMDEVs LT.—-15)NDEV=-15 

BeGiDeEy Gl. 15)NDEV=15 

RETURN 

END 
Hew KKK KKK KKK KKK RK RAE K KKK KKK KKK KKK KE 
C Vimeo bORPUUorD TO SELECT A BORIC ACID HIGH 
VOLTAGE 
2 aioe 


REAR KAKKKRKAKEKEKKEKREKKEREKAKHKKEEREKEKAKKEKEKRAKRAK KS 

Sieur tit sSBOkKPUGCCURLO, NDEV,OFSET,OFSET!) 

ALOGC=ALOG10( CURLO) 

LaGwweEy.brks0)GO TO 19 

DEV=88.16434(-215.2235* (ALOGC) )+(188.225* (ALOGC**?2 ) )+ 
(-69. 

13469" ( ALOGE**3))+(7.4695* (ALOGC**4 ) )+(1.65257*(ALOGC 
**KS))+ 

2(-.347786*( ALOGC#X%6 ) ) 

NDEV=iF LX CDEV )+] 

IF(CNDEV.GT.-19)NDEV=19. 





1 
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0) DEV=ABS(FLOATCNDEV) ) 
EPCNDEV.LT.O)DEV=DEV=15 
ao leet bit. 11 34579*DEV+(—.0012654*(DEV**2 ) )+(.00 


m2 13% 
1(DEV**3)) 
Ost |=OFSET+. | 
MONDE V.LE.—19)NDEV=-19 
PeGvoEvVeGr. 19)NDEV=19 
IF(CNDEV.LT.0)NDEV=-NDEV+15 
PMVeNDEV.GT.0O)NDEV=NDEVFI 5 
RETURN 
END 
KKK KEKE KERR KKK KK KRKKRKKRKEKEKEKKS 
C SUBROUTINE SBRKR;USED TO SELECT A THERMAL MAGNETIC 
BREAKER 


Sok KKK KKK KKK KREKKRKEKKREKKKKKKEKRKEEKKKKKKKEEKSE 

Sloe ourtNne —SBRER(CURLO, NDEV,OFSET,OFSET] ,AMAG) 
DIMENSION CON(S) 
CONC 1 )=2.6 
CONG? )=-3.5 
CON 3)=4.6 
Ole ) = 516-6 
CON(5)=7.0 
HYPE. 10 

10 FORMAT(’ ENTER MANUFACTURER (N.A. AT PRESENT USE 

FOR FUTURE 
DEANS LON) ) 


er ke 2.5 
2D POC ENPUT TEMPERATURE FOR AMBIENT SHIFT,1 FOR 
ifieoo LHAN 


Mme meme sr Oh sb0 DEG) C,° //° <CR>FOR NO SHIFT...’ ) 
POC rinm2s. [OFF 

28 FORMAT(1L) 
iPiv@werst Os.) OF S-—. | 
Po@nor.hO.2)OFS=.05 
Cera. NE.0)GO TO 80 
PeGGmriEO.lLE.600.)GO TO 30 
PeecuwLOuwLh.s00.)GO TO 50 
Com non 70 

30 Le GewrRLO.LT. 300. )MDEV=1 
DPiecurrOs. Ll.300.)OFSET=2.47712 
Beem nOnel. 200. 0.AND. CURLO.LEs350.)MDEV=2 
GimbOre DatOO. 0 mA De CURLO.LE.350.)0FSET=2.54407 
PGC OG. 350.0.AND.CURLO.LE.400.)MDEV=3 
Pi eemrnOnG laos. ANDs CURLO.LE.450.)O0FSET=2.60206 
PEeurnOmamtegOiOn. Os AND. CURLO.LE.-450.)NDEV=4 
PE Geto e Teahouse AT. CURLO.LE.450.)0FSET=2.65321 
Peet nonG he 50.0.AND.CURLO.LE.500.)NDEV=5 
Mitnienenom@ec 5). AND. CURLO.LE.500.)OFSET=2.69897 


ote + ee eee ee oe 





a 


IF(CURLO.GT. 500. )NDEV=6 
PiCGURIGeGr ~500-)O0FSET=2.77815 
Gomero) 26 

50 MCCUIEOm IEE. 7001.) MDEV=7 
PaeCcuRinOn LE.700.)OFSET=2.8451 
LF(CURLO.GT.700. )MDEV=8 
mi aGumironel.700.)OFSET=2.90309 
GO TO 80 

70 LF(CURLO.LE.1000.)MDEV=9 
MNCGUREO.LE. 1000. )OFSET=3.0 
IF(CURLO.GT.1000.)MDEV=10 
PPGclmmoreT. 1000.)0FSET=3.3 

8.0 NDEV=MDEV+34 
OF SET=OFSET+OFS+.1122 
OF SET1=OFSET+. 12 
TYPE 90 

90 FORMAT(’ DOES THE DEVICE HAVE A MAGNETIC 

ELEMENT?(0)YES.’ ) 
ACCEPT 95, NOEL 

95 FORMAT(L) 
FCUR=4*CURLO 
SET=10**(OFSET-.12) 
DO 100 1=1,5 

TD=SET*CON(1) 
PGnieecrereiun co. 10.105 


100 CONTINUE 
OS ANAG=ALOGIO(CCON(T1) ) 

LF(CNOEL.NE.O)AMAG=0.0 

RETURN 

END 
KREAEKKKKKKEKKEKRKKKAKKAKKKE KEK KEKE RAK KE KERKEKEKKEKKRK AK 
c Toole, SKELA;USED TO SELECT A RELAY 


RAAEKEKKKKAEKKRKEKKKRKEKAKEAEKEKKKERKAKKKEKAKRAKREARRAKEREKK 


SUBROUTINE SRELA(CURLO, NDEV,NTAP,NCT, PKUP, CT) 


TYPE 10 
10 FORMAT(’ SELECT MANUFACTUER (FUTURE EXPANSION)’ ) 
TYPE 20 
20 FORMAT(’ SELECT THE WESTINGHOUSE CO RELAY ,REFER TO 
TABLE’ ,$) 
ACCEPT 25,NUM 
25 FORMAT(L) 


NDEV=-NUM 
Pee@eurto.LT.30)CT=10 

Maeeiimiloech. 30.AND.CURLO.LT.120)CT=30 
PaGerolo.GE.120)CT=120 
Peete oO. LC PNCT=1 

PaG@eme: oO. 30) NC T=2 

PaGehenO .l20)NCT=3 

Toe cunLo, cr 


ae 
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PCr APeieel.«0)GO TO /] 
Pear ene. le 2)GO TO 72 
Mit ApeLEewt.so)GO LTO 73 
Me hAPSie.<2.-0) CO TO 74 
een Arah nn 2. 5)GO TO 75 
WPCTAP. LE.3-0)GO TO 76 
PelAP bh. 565)GO TO 7/7 
PeGtArPath.4.-0)GO TO 78 
PRCTAPSIE2.5.0)CGO TO 79 
PGi Arehh.-6.0)GO TO 80 
DCIAPeLEe/-.0)GO TO €} 
terre i. )GO LO G2 
PeCrAre bE. tO.0)CO TO 83 
cO TO 84 

7} NTAP=1 
TAP=1.0 
COLO” 90 

2 NTAP=2 
TAP=1 32 
COTO “90 

73 NTAP=3 
TAP=1.5 
COsrO 90 

yay NTAP=4 
TAP=2.0 
CO-1TO 90 

Ss NTAP=5 
TAP=2.5 
Coe 1o 90 

76 NTAP=6 
TAP=3.0 
COMMON 9 0 

af NTAP=/7 
TAP=3.5 
ich lO ee) 

78 NTAP=8 
TAP=4.0 
GO, 10° 90 

79 NTAP=9 
TAP=5.0 
COssEO. 90 

8&0 NTAP=10 
TAP=6.0 
CORT. 9 0 

ot NTAP=1 1 
TAP=7.09 
GO TO 90 

82 NTAP=1 2 
TAP=6.0 
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Go TO 90 
muAP=13 . 

maP=10.0 

GO TO 90 

NTAP=14 

maP=)2.0 
PKUP=ALOGILOC( TAP*CT) 
RETURN 

END 








C 
kk RK KKK 


i) 
20 


AS 


30 


40 
45 
510 
oD 
60 


70 


i) 


80 
400 
410 


420 
440 
445 
450 


Peete wells r 4 PLOTS DEVICES FROM COORD.F4 
KARKKKEKKKKEK KEKE KKE KIKI KEK KKK KK KK ARK KRKKKKKKKKEKKEK 


OPEN(UNLT=30, FLILE=’ ADEV.DAT’ ) 
OPEN(UNIT=31, FILE=’DEV.DAT’ ) 
OPEN(UNIT=32, FILE=’ RELA. DAT’ ) 
OPEN( UNIT=33, FILE=’ FCUR. DAT’ ) 
OPEN(UNIT=34, FILE=‘TEST.DAT’ ) 
EXTERNAL T25HAN 


ee 


Cun selwenesu) NEV 20,0/4),AILDEV(100,7),RIDC100) 


COMPLEX ALDEV,RID 
Rima? , 3 )R 
FORMAT(F) 
Reca 2, 10 ) I 
FORMAT(CI) 
Pele. 0)GO TO 5 
Pees 25 20) 1,51 
FORMATCI, F) 
iene. 0)GO TO 25 
EBC L)=E1 
COMO 15 
HOw40 K=1,100 
Romney, 30), (ALDEVC(C1,J);J=1,7) 
FORHAT(1, 14F) 
Peereee«.0)GO TO 45 
CONTINUE 
RED 345 50)1B,E 1 
FORMAT(1,F) 
Rie S460) NN, 1, J,K,R,X 
FORMAT( 41, 2F) 
IF(CNN.EQ.0) GO TO 70 
NE=NN 
CI=I 
CJ=J 
RIDCNN)=CNPLX(CI, CJ) 
GOset O55 
DOr. 0 <1 =).,:20 
Mewes ley ok, DEV CL, J), J=0, 4) 
FORA T CL, SF) 
MiGeHos0)GO TO 400 
CONTINUE 
Tare 1 0 
POA  DEVECE LISTING’ ) 
bomo40) f=1, NE 
iver 420, b, CAIDEVCI,J),J= 
Peele azo... (ALDEV(L,J),J 
BORAT Gi, 14F ) 
CONTINUE 
Tees 50 


BEST) 
=1,7 


iD 


Perera VlLCEo HAVE BEEN SET LF YOU WANT TO PLOT 
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GRAPH! PAPER 
fee O jelr NOT TYPEC!)” ) 
RCCERT 456-1 P 


456 FORNATCI1) 
iwi.) )GO TO 480 
Hyer be 46.0 
460 DOM pviGGemetiiomw Ube PEROT CURVES BASED TO THE LOWEST 
VOLTAGE 
[eliheernbe SYSTEM.” ) 
EYP he 46 1 
461 HOW wer Ute tHe CURRENT OFFSET FOR THE 


GRAPH, INTEGER..’ ) 
ACCEPT 462,1AOFF 
462 FORMAT(1L) 
Granta Lr 
CALL LOG(IAOFF) 
GALL FRUSH 


430 Tee GAs 5 

485 P@tatr( ©NPUT ELEMENT ..” ) 
peer rT 490, 1EL 

490 FORMATCI) 


Pepi xCREALCRIDCLEL) )) 
J=IFLIX(CAIMAG(RIDCIEL))) 
EB1=EB(1) . 
DOm4o5 > K=1,18 
EB1=AMIN1I(EB1,EB(K)) 
495 CONTINUE 
550 FAC=ALOGIO(EBC(1)/EB1 )-IAOFF 
FACI=ALOGLOCEB(J)/EB1)-LIAOFF 
GUR=REALCALDEV(IEL, 3) ) 
Chel REALCALDEVCIEL, 4)) 
Clio =ALTMAC(ALDEY CIEL, 3)) 
CUR = AIMAGCAIDEV(LEL,4)) 
KERR KKKKARKRKKKEK KKK KAKA KKEKRKREKEKEKEKKEKKKKEEKEKRKKKEKEKS 
C NOW TO PLOT LOAD SIDE OF THE ELEMENT 
KRARAKRKKKAEKKKKKKKEKRKKEKKEKKEKKAEKKEKEKKKKKKEKKAEKRKRKEAEKKEESE 
Clipe ema Cl. LEL,FALT,900.0) 
Geena ory CLE, 2)).GT-0) GO TO 585 
Se- ReEALCALDEVCLEL, 4) ) 
SpA / 107*STL 
Pew Celpny lh. 5) ).-NE.0)STP=AIMAG(ALDEV(IEL,5)) 
Sel OGTOCEALT ) 
YM=(STL+FAC1)*1904+175.0 
Whi O10 
CALL MOVEA( XM, YM) 
Se ese 
iver SEP) 
Laine tel s ) TYPE 575 
575 Pot ee tiih OAD  STDE RELAY WILL NOT BE VERY 
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peetoLITIVE 


D7 0 


530 


Dic 


eG uner Allen ) 
mv = eenmbis CALTMAG AT DEVCIEL,2)))) 
Evie? 0. pe y, LTP 
BOMMATC DEV, ° ,21) 
DEE-1 TP=)5. 
Pitt CULE / 200.0) 
ice ) it D=) 
TVAL3=10.0 
Toa —= AT HIAGCALDEVCLEL, 6) ) 
bOms smell, 1TP, 1D 
Re=. 1 *( FLOATCEK )) 
Croom Ay CrbDEy RR. TVAL, TDSET, DEV) 
ci wencie 1 VAL )TVAL=TVAL3 
Ty Abs ov Al 
T=900-( ( TVAL+2.0)*180) 
R=(ALOG1O(RK)+FACI1+STL)*190+175.0 
GALL DRAWACT, FR) 
CONTINUE 
PAUSE 4 
Perm new hDEV(LEL,5))-~-NE-O)CALL DRAWA((900.0),R) 
C000 591 
OPSEIT=CUR I 
M=IFIX(80*OFSET) 
PVN TAG CALDEV( LEL, 2) ) 
Mee trerx  100* (ALOGIOCFALT) )) 
DO 587 K=M,KIM 
CVAL—=. => 201 
mae. > yOALL FUSE ( 1, ADEV, CVAL, OFSET,TVAL) 
WGA EV 6G E35 )CALUL 


Miers R( 1, IT FIX(ADEV),CVAL,OFSET,TVAL, 


WAIMAGCALDEVCLEL, 5) )) 
Tae 1) Alea 2. © 
Mah Abeer. IVA ANDOND.NE.O0O)GO TO 5389 
TVAL2Z—-TVAL 
ND=1] 
CALL 


DRAWA(900-( (TVAL)*180),(CVAL+FAC1 )*1904+175.0) 


Dio / 
yoo 


CONTINUE 
OF SET=CUR3 
CA NOVEACO.0,0.0) 
ND=0 
M=IFIX(80*OFSET) 
DO 588 K=M,KIM 
CVA ote «0 1 
Dewees anhL PuSh( 2, ANEV,CVAL,OFSET, TVAL) 
[OE ee eveerG hie) CALI 


@rPekKR (2,1 FIXCADEV),CVAL, OFSET,TVAL, 


PAM UAGCATDEVCIEL, 5) )) 


i _. Ce ee," Ge ee, ee, eee 


a 


. 
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EVAL TVALA2..:0 
Peo vAeeie ty aee AND.ND.NE.O)GO TO 591 
TVAL2=TVAL 


ND=1 

CIMGIE 
DRAWA(900-(CTVAL)*180), (CVAL+FAC1)*190+175.0) 
588 COUT. NUE 
SeKKKKKKHKKKKRKAKKKKRKKKEKERAKMMKKe KKK KKK KKK 
G NOW TO DRAW THE LINE SLIDE 
KKeEeKK KK KKK KK KKK KKK Ke KEKE KEKKEKKK KKK 
591 GrewenOvEA€ 0.050. 0) 


een mAb CALDEV(CLEL, 2))) 
PerewlDsGl.0)GO TO 595 
MeEV=EPF LX(-AlD) 
Sig heat DEV LEL, 3) ) 
Se —F ie) 1 Oe STL 
Powel DEV Clb > )).NE.0)STP=REALCAIDEVCIEL,5)) 
YM=(STL+FAC)*190+175.0 
er 
CALL MOVEA(C XM, YM) 
Siee=S1IP=10 
EER = brol ees ip ) 
Deere — Ns 
Pv teh. (Dley 200.0) 
ieee VE 0) 1 D= 1] 
eof O. 0 
Pirie nbALCALDEVCLEL, 6) ) 
Dome k=15,1TP, LD 
Rik = <r) 
Cm wigey cP DEVeRK. TVAL, TDSET, DEV) 
Perv i Gim lt VALS) EVAL=TVAL3 
TVAL3=TVAL 
ial 
NDRAWA(900-C CTVAL+2.0)%*180), CC ALOGLOCRK)+STL 
1+FAC)*1904175.0)) 
590 CONTINUE 
Ma@Gee Orbe y CIEL, 5)).NE.O0) CALL 
DRAWAC(900.0),(CCALOG1OCRK)+ 
ISTL+FAC)*1904175)) 
595 Mi@wirwne@wl Dey CLL, 2)).LT.0)GO TO 598 
OFohr=CUR 
Miiev—REALCAIDEVCLIEL, 2)) 
ND=0 
M=IFIX(80*OFSET) 
KIM=LFIXCLOO*ALOGIOCFALT) ) 
DO 597 K=M,KIM 
CV k =. ] 
ie Wives ON PUSECE, ADEV,CVAL,OFSET,TVAL) 
imG@evev.Gh. 395) CALI 
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@ereKR(1,1FIXCADEV) ,CVAL, OFSET,TVAL, 
Lean ADE VCLEL, 5) )) 
TVAL=TVAL+2.0 
Gn Serene AL SAND. ND.NE.0)GO TO 601 
TVAL2=TVAL 


oa 
Ata 
HeeewAC9O0O0-(TVAL*180) ,(CVAL+FAC)*190+175.0) 
597 CONTINUE 
601 OFSET=CUR2 
Gree ioOvEACO. 0.0.0) 
ipo 


fer 1X CS. 0%O0F SET) 
DO 599 K=N,KIM 
CVAL=K*.01 
ew ike O44) CMeIemh USEC? , ADEV, CVAL, OFSET, TVAL) 
Pen DEV wet. 35) CALL 
MremmrkR( 2,1FIXCADEV),CVAL,OFSET,TVAL, 
IREALCAIDEVCIEL, 5) )) 
TVAL=TVAL+2.0 
Mime b2e bh. TVAL.AND.ND«NE.O0)GO TO 598 
TVAL2=TVAL 


ND=1] 

CALL 
PeeewACIO0—(TVAL*180),(CVAL+FAC)*1904175.0) 
9 CONTINUE 
598 ieee G 0:0 
600 Pola CORT NISHEN??7° ) 

ND=0 


CALL MOVEA(0.0,0.0) 
CALL FLUSH 
ACGE PT 610, 1 


610 FORMAT(1) 
icine )co TO 2000 
GO TO 445 

2000 CLOSE(UNIT=30) 


CLOSE(UNIT=31) 
CLOSECUNIT=32) 
CHOSE CUN1T=3 3) 
Closet Uni T=34 ) 
Obie I NITECO,0) 


> OP 
END 
KAKA KE AARKKKR KKK KRKK AK KKK EK KKK KKK AKER KAKKKKEKKKAKE 
C Peeewb REPALTs: READS FCUR.~DAT TO OUTPUT FAULT 
SURRENTS 
C Nelo lege ho. ts 


pes k wR SOX eK oe oe KK KK KK HK KK oe KK ve HIE KKK oe de oe eK 
eo ENE RE PALMCIPALT, NEL, FALT, CUREL) 
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10 READ(33,15)1LEF,NN, TFAL, SFAL 
15 FORMAT(21, 2F) 
MiGlEE.EO.O)GO TO 20 
Pieper sie elrALL)GO TO 10 
MoweleenbL)yGO TO 10 
PAGE L= TF Al 
PALTZ —-SFAL 
EAcu—-AM ENT CFALT1, FALT2) 
PEOGUR=3-.0*CUREL 
PiGnevwet weer LOCUR FALT=AMAX1(FALT!,FALT2) 
FALT=AMAX1(FALT1, FALT2) 
20 REWIND 33 
Re URN 
END 
DEREEAEeC KARA AK AKA A KKK AKKHREKRKKRKKEKEAKKKAKRKKARK 
C SUBROUTINE FUSE;COMPUTES TIME VALUE FOR A GIVEN LOAD 
CURRENT 
PRAEKAKKKKEKKARKEKKKKKKRAKKKKKAKKKKEKAKKEKEAKKKRKKAKKKAKKERKS 
SUBROUTINE FUSE(ID, ADEV, CVAL, OFSET, TVAL) 
TVAL=0 .0 
vB Y= lp lx CADEY ) 
MAGE vec iets )GO TO 50 
X=CVAL-OFSET 
PeGuve Neel yGO. TO 40 
ie eae ee Oe? 4X )+0223~.7526*(X**2 ) )+(-751.279% 
eo 3 ) ) 
1401347 .-01*(CX**4 ) )4+0-1222.67%(X**5) )4+0440.813%C(X**6G) ) 
PVAL=TVAL=2..0 
GO TO 100 
40 K=CVAL—OFSET+. 1 
lien oa 7 O22 Je C121. 909*(X**2 ))4+(-213.420*(X* 
*3))+ 
1€201.771%*(X**4 ) )+(-96.9836%(X**5) )+018.6046* (X**6) ) 
TVAL=TVAL-2.0 
Gets le 4) TVAL=—1.95 
COM GO a0 0 
50 nC EOF SET 
PGi!) GO TO 70 
eer 7 OFC — 45 7262768*K )4+(1.21940* (X**2)) 
TVAL=TVAL-2.0 
GO TO 100 
70 X=CVAL~OF SET+. 1 
eee oa (=O 173224" X)+(3.231806%(X**2 ) )+(-.6057 
136*(X**3)) 
TVAL=TVAL=-2.0 
100 RETURN 
ERD 
KRAEAKKKKEKEARKKAKEKKAKRAKK RK KKKAKRHKKAKAKRHKKAKKKAKEKN 
c SUBROUTINE RELAY;COMPUTES TIME VALUE FOR A GIVEN 
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CURRENT 
C Tae eeOkRmANSOPECTFIC RELAY 
KKK KKK KRKKHK KKK KKKEKKKKRKKAKREKKKKKKRRKKRKEKES 
SOMME RELAY CNDBY, CURR, TVAL, TDSET, DEV) 
DitveNtolOn! DEV(20,0/4) 
Vie —Or0) 
hoe ne COUR R ) 
ite ere tel. 2) X=]. 2 
DOmwLO 7 1T=0 ,4 
TVALS=IVALIDEVCNDEV,1)*(X**1) 
10 CONTINUE 
IF(NDEV.~EQ.7) TVAL=TVAL-2.0 
ie iv Abst tT DSET 


RETURN 
END 
Kk ek kK ek aK RRA KKK KK AK KKK KK RAE KAKKAKKKAEKAK 
C (Pee  COGelLOG, Makes LOG LOG FORM FOR GRAPH OF RELAY 
CURVES 
HK KKK KKK KKK KKK KAR K KARR KK KK KER KEK KKK AKER KKKKKKKEKEKE 
C OUTLINE BORDER 
KRKK KKK KK KEKE KKK KEK KKK KKK KEK RE KEK KKK KK AK KKK KI 
C MUSietsh ALOG. LOG,@SYS:GRALIB 


wekkiokkekakeekekeeKKKKEKKKRKKKKKCKKRAKKAKKREKEEKKEKKEKRKEKKEKEKKSK 
SUBROUTINE LOGCIAOFF) 
CALL ERASE 
CMO veEAC 0. 0.2.0.0) 
Crm RAWACO.O., 750.0) 
CMR AWA 1000. 0,750.0) 
Cie AmACTO00.0,0.0) 
Cr weremenw NCO .0, 0. )) 
Gee Owe AC 900.0, 5.0) 
ew prea 900.0,750.0) 
DO 200 N=0,4 
Domo Mie. 10 
DN=10**N 
XM=iH+ DN 
Mee. 61235000) GO TO 225 
XX=CALOGIOCXM) )*180 
XTT=900-XX 
Cramipeiov HAC XTT, 175.0) 
Gemie DRAURCO.0,575.0) 


100 Corl due 
200 CONTINUE 
DOS CALL MOVEA(900.0,175.0) 


CeebieopiwAGe. 0,175.0) 
Dome ei. 2 

DOM Oe), 10 
DN=10**x 

XM=H* DN 





S010 
400 


460 


500 
p10 
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ho ALOG TOC XM) )4ro0 

= Xk 

CMe MOVE ACS 5.0, x7) 
Abba (805.0,0.0) 
CONTINUE 

CONTINUE 

Owe oO H=0, 3 

X¥N=1754+190*M 

CA MOVEAC 900.0, XN) 

Cobia NuMBERC1O,~ 12° ) 

Cone woweR CO.0, 10.0) 
J=IAOFF+M 

CAni NUMBERC J, LI ) 

CONTINUE 

Gwin GHRROT (C9 Q) 

Guyer MOVEACI60.0,400.0) 

Cure eANCHOSt CURRE ,5) 

Co mieeatelnos( NY © . 4) 

Cybi VN COS( NWN Are’ ,5) 

Cine ANGHUOSC’S”* ,1) 

Cupra noVveEn( 900.0,0.0) 
GCywblewieawAC900.0,750.0) 

DOP >D00 4=0,5 

xN=900-180*M 

Crm MOvEAC XN, 150.0) 
Peete.) ) CALL NUMBER(.01,°F3.2’ ) 
BeGuemOn 1) CALL HUMBER(.1,°F2.17) 
iim. 2 CALL NUMBER(1,°117) 
Geno. 5) CALL NUMBERC1I0, “12° ) 
Meatio.4 ) CALL NUMBERC100,°13° ) 
Machi. bO.5)GO TO 510 

CONTINUE 

SA bletovubnaA(20.0,140.0) 
CebeoaeeER (1000, ° 14°) 

Gb CHEROT CO) 

Cr ereeiovEat(200.0,/55.0) 

CA ANGHOS( TIME °,5) 
Cribberncuos(t [NSE ,5) 

Cre ANGCHOS(  CONDS’ ,5) 

CEP MOVE ACO.0,0.0) 

RETURN 

END 


RARKAKKRKKKKKKKKAREKRKEKRAEKKKKKKARKEKARKKAAKRKKS KKK KKK 


c 


CURRENT 


C 


Vo Oieikee cr PB RR= COMPUTES TIME VALUE FOR A GIVEN 


VALUE 


KKKEKKEKKRKEKKEKAKAKEAK KK KAKARKEKAKEKAAKAKKEAAKAKKEKEKS 


imROUnmint wCkKTRKR( 1D, MDEV,PVAL, OFSET,TVAL, AMVAL) 
CVAL=DVAL-OFSET 








150 


PeGiiawe. | )GO TO 50 

Cee 6 Vy hit. 1 122 

eee eer NiVAL. AND. AMVAL.NE.9)GO TO 90 

Mmaeeve ter .0')GO TO .80 

eo a Oot(-)OeaooZols(CVAL) )+023.65419* ( CVAL*®*2 ) )+ 
eG. 7922 

L*¥(CVAL*®*3))4+033.77472* ( CVAL**4 ) )+(-21.14536* (CVAL¥*5 
))+(6.24 

Poly foe CCVAL**G ) ) 

ive Al? . 0 

CO LO 100 
50, CvVAL—GVAL+. 2 32 

MCE v wy. il .0)GO TO°S0 

AMVAL=AMVAL+.075 

Pcie eG AMV Ao AND Se AMVAL.RE«0)CO TO 95 

Myo G5 067 224-4 O94 335 CVAL)+(0—.50690284*(CVAL**?2 ) 
)+€10.40 

12739* (CVAL**3) )+(-17.064956%* (CCVAL**4 ) )4+013.011543* CC 
VAL**5)) 

2+(-4.010724* (CVAL**6) ) 

TVAL=TVAL-2.0 


Gow TO 100 
80 TVAL=6.0 
Coto, 150 
90 TVAL=-2.0 
Coe ro! (100 
95 iMiri-ALOCTO(.025) 
100 CONTINUE 
150 RETURN 
END 


MKKEKRRAEKKKKKKKKAREKKKEKRKEKKEKAKKKKKEREKEKARKKKRKKEER NK 
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eT 


PL 


ARPBENDIX -C 


DESK -FILE FORMATS 








2 


Disk file RES. DAT 


a Double precision wire resistance coefficients 
ag Double precision wire resistance coefficients 
a Double precision wire reactance coefficients 
ae Double precision wire reactance coefficients 


biek file TEST.DAT 

element number/line bus/load bus/positive seq p.u. impedance/zero 

seq p.u. impedance/transformer neutral imp line bus/transformer neutral 
imp load bus 

00000 

Load data 

bus number/p.u. VA/per unit subtransient reactance 

0000 

Disk file XFER. DAT 

line bus/load bus/inrush/low withstand/high withstand/transformer size 
mCuUR. DAT 

line bus/load bus/current at line bus/current at load bus 

FCUR. DAT 

fault bus/element number/3 phase fault/1l phase fault 


eV . DAT 


device number/a /a,/a,/a,/a, 
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ADEV. DAT 


element number/element bus numbers/device numbers/line bus element 
low pickup, high pickup/load bus element low pickup, high pickup/ 
magnetic settings line, load/time dial setting line, load/CT ratios 
line,load/ 


ee 
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APPENDIX D 


SAMPLE PROGRAM DATA INPUT 








Ika), 


Before beginning the program sequence, RES.DAT and DEV.DAT 
contain the proper coefficients listed in Table 8.5. All source 
programs have been loaded and saved. PLOTD.F4 uses subroutines con- 


tained in the PDP1O system files so it must be loaded using: 


LOAD PLOTD.F4, @ SYS:GRALIB. 


The first program is started with run INPUT. Underlined words and 
Numbers are data input by the user. Words not underlined are prompts 


from the computer in the input sequence that follows. 


Do you want to read from a KYBRD? 
—Ch> tem empue from TlyY, | for read. 
<CR> 


Enter the problem base in MVA 
i 

Enter the problem base in KV 
15 


Bitwemmc bememt data as follows: Bus ( ) to ( ) 
(Length) (Size) (Transformer impedance (P.U.)) (Transformer size 
in MVA) (Transformer High Voltage in KV)(Transformer low 
voltage in KV) 

1 2 3000 6 -(element with no transformer) 

Pare Ome OO le 00> 1 1594.1 —(element with a transformer) 
lWpteeeranstormer ZN in P.U. refer to BUS 


Side, (BUS 2) (BUS 3). 
.001 .005 .001 .005 


(This is continued for each element until all elements of 
the type in RES.DAT have been entered) 


0 -(tells computer no more elements) 
At this point if another wire type is to be used the 


following procedure applies if not skip this section 
and proceed to 
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Input Bus 1 voltage. . .0 
Data is accepted and a list follows 


(Here the computer types all computed information on the 
elements entered so that they may be checked for error.) 


Now enter the load data in the following manner;. . 

Cams) (Cboadein KVA)G@Pwr Factor). . .or. . Waec roccr in HP) 

Ceci det@m iw. 205... bus) (Load in Amps) (Pwr Factor) 

What method will you use?. . .VA, HP, or Amps? 

=) -(writes element data computed to disk) 

End of execution 

(Load RES.DAT with the new wire type coefficients) 

RUN INPUT 

Do you want to read from a KYBRD? 

“Cn fOr input trom a TTY, L for read. 

1 

Old data is entered ready for new line data 

Enter element data as follows. . .(Program runs as 
imeehne first section) 


Picking up the program from completion of element data, 
the data entry follows. 


ImpursBuceee Voltage, . .15000 
Input Bus 2 Voltage. . .15000 
Input Buses Voltage. . .4100 








yoni s sere 8 Voltage. ~440 
Data is accepted and a list follows 
(Computer types computed element data). 


Number of buses is 8 

Now enter the load data in the following manner;. . 
CBucyioadeim KVA) (Pwr Factor). . <or. . Warton ioad: in Hie) 
(vende tom) Or-.  . . (bus) (Load in Amps) (Pwr factor) 
What method will you use? KVA, HP, or Amps? 


HP 

Enter. 

4 800 .80 

What method will you use?. . .KVA, HP, or Amps? 

He 

8 500 .90 
(This is continued until all loads are entered? 
ViAeeenictnoOdewill yoususe?. . .KVA, HP, or Amps? 

S 

End of Execution 

RUN LDFLO 

End of Execution 

RUN ZBUS 


Enter the zero sequence impedance for the loads at each bus. 
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8 .001 .005 -for each load where the zero sequence impedance 

: will not be assumed to equal the positive sequence 
: impedance 

QO -tells the computer there are no more entries 

Enter the mutually coupled elements in per unit in the following 
format, TEST.DAT has the element numbers. . ., (Element) (Element) 
(Zero Sequence Value) 

2 3 0 .001 (The zero must be entered because the coupling is a 
reactive value) 

0 -(no more or no coupling) 

End of execution 

RUN FALT 

Enter the fault impedance (Real) (Imaginary) 

00 

End of execution 

RUN COORD 

All load currents are ready, use device number list 

Select device for the load at bus 4 

4 


e 


Select the device for the load at Bus 8 

D 

Input temperature for ambient shift, 1 for less than 
IMG wwecweam. for 60 Dep C 

Cl erOmene Shit t 

<CR> 

Does the device have a magnetic element?(0)Yes 

0 

Select device for Bus 6 protecting element feeding load at Bus 8 
Z 

Select the Nestimenouse Co. relay, refer to Table 1 
Does the relay have an inst. element if yes 0 

0 


— 





(Until all backup devices have been set). 
Enter the device for the Bus 6 of element 5 
4 
lo necrmeicaewectinghiouse Co; “relay, refer to Table 1 
Does the relay have an inst. element if yes 0 
0 


— 


(Until all devices in the system have been set). 
Device listing ~(Here the computer lists the matrix ADEV). 
End of execution 
Rune PLO 
Devices have been set if you want to plot graph paper O if not 
type l 
0 








158 


This will plot curves based on the lowest voltage in the 
system 

imapubeenewcurment OLfSei for the graph, integer 1] 
Input element number 7 

qT 

Finished? 

limos 0 yes ) 

Devices have been set if you want to plot graph paper 
Ome not type 1 

i 

Input element number 

6 

Finished? 

© 

End of execution. 


If a device did not plot as desired its settings can be changed 
ii eADEV. DAT . 
































